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(2)  An  examination  of  the  spatial  range  of  the  molecular  binding  forces  indicates  that  a 
uniform  coating  one  or  two  molecular  layers  thick  is  sufficient  to  saturate  the  molecular 
forces  so  that  the  cohesion  is  determined  entirely  by  the  coating.  (3)  Binding  from  cap- 
illary forces  of  a liquid  bridge  linking  two  spheres,  if  it  occurs,  is  comparable  to,  or 
greater  than,  molecular  attractive  forces.  However,  such  liquid  bridges  may  not  be 
formed  under  the  experimental  conditions  of  interest.  (4)  Electrostatic  binding  appears 
to  be  negligible  except  under  special  circumstances  or  for  statistically  few  particles  in 
a large  collection.  (5)  Desirable  coating  materials  are  fatty  acids,  such  as  stearic  acid, 
and  fluorinated  hydrocarbons  and  polymers,  such  as  teflon.  Light  scattering  from  metaL 
lie  spheres  coated  with  various  nonabsorbing  and  absorbing  coatings  has  been  calculated 
with  dispersion  included  in  the  optical  constants  of  both  the  metallic  core  and  the  ab- 
sorbing coatings.  The  principal  results  of  this  study  arc:  (1)  The  scattering  from  an 
uncoated  metallic  sphere  is,  on  the  average,  within  three  percent  of  the  scattering  from 
a perfectly  conducting  sphere  of  the  same  radius.  (2)  For  a broad  range  of  coating  thick- 
nesses and  refractive  indices  greater  than  approximately  1.6,  the  scattering  from  the 
coated  sphere  is  similar  to  that  from  a perfectly  conducting  sphere  having  a radius  equal 
to  the  total  of  the  metallic  sphere  radius  plus  the  coating  thickness.  (3)  With  thin  coat- 
ings the  change  in  the  scattering  due  to  the  coating  is  linear  in  the  coating  thickness  and 
less  than  three  percent  for  coatings  having  a ratio  of  the  total  coated  sphere  radius  to 
the  core  sphere  radius  equal  to  1.004.  (4)  Absorption  bands  in  coating  materials  have 
been  investigated  by  the  use  of  Lorentzian  lines  of  varying  linewidths.  A Lorentzian  ab- 
sorption line  results  in  a localized  perturbation  in  the  scattering  curve  with  the  width  of 
the  perturbation  approximately  equal  to  the  absorption  linewidth  and  the  magnitude  of  the 
perturbation  determined  by  the  coating  thickness.  (5)  Scattering  from  an  actual  metal- 
oxide  coating  has  been  examined.  Large  magnitude  localized  deviations  in  the  scattering 
curve  are  observed  at  points  where  the  real  part  of  the  coating  refractive  index  is  small 
and  the  imaginary  part  is  large.  (6)  Scattering  from  small  dielectric  particles  is  less 
efficient  than  from  metallic  particles  of  the  same  size.  A localized  valley  in  the  scatter- 
ing is  found  at  a point  where  the  imaginary  part  of  the  refractive  index  changes  abruptly 
from  one  slowly  varying  value  to  another. 
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An  analysis  of  the  mechanisms  for  the  cohesion  of  small  metallic  particles 
suggests  that  the  binding  of  metallic  particles  can  be  reduced  by  coating  with 
one  of  several  different  dielectric  materials.  The  effects  of  dielectric  coatings 
on  the  light- scattering  properties  of  metallic  spheres  is  determined. 

Mechanisms  for  the  cohesion  of  particles  can  be  divided  into  three  catego-  ; 

ries:  (1)  If  two  particles  are  linked  by  a liquid  bridge  that  wets  both  particles, 
the  surface  tension  of  the  liquid  can  yield  an  attractive  force  between  the  parti- 
cles. (2)  Two  particles  in  close  proximity  experience  molecular  attraction  and  j j 

molecular  binding.  The  cold-welding  of  metals  at  points  of  intimate  contact  is  an 
example.  (3)  Electrically  charged  particles  are  attracted  to  oppositely  charged  | 

particles  and  to  image  charges  in  uncharged  particles.  | I 

The  principal  results  of  the  cohesion  study  are: 

• Tlie  molecular-attraction  mechanism  appears  to  be  the  dominant  mech-  I ! 

anism,  witJi  the  strengtii  of  the  binding  dependent  on  the  chemical -bonding  ' 

mechanism  appropriate  to  the  coating  material.  Pure  metals  are  expected 

to  have  the  largest  binding  forces,  and  nonpolar  molecular  solids  the  smal- 
lest. Intuitively,  covalent  and  ionic  materials,  such  as  metal  oxides,  are 
expected  to  have  binding  forces  that  are  sensitive  to  the  alignment  of  the 
lattices  on  the  two  particles  and,  hence,  to  have  intermediate  strength 
forces.  1 

• The  range  of  molecular  forces  has  been  examined.  It  is  estimated  that  a 
uniform  coating  one  or  two  molecular  layers  in  thickness  should  be  suffi- 
cient to  saturate  the  forces  so  that  the  cohesion  is  determined  entirely  by 

the  coating  rather  th.in  the  core.  ! 

• The  electrostatic  binding  mechanism  appears  to  be  negligible  except  under 
such  special  circumstances  as  those  in  a corona  discharge  or  for  statisti- 
cally few  particles  in  a large  collection  of  particles. 
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• Attractive  forces  from  the  liquid-bridge  mechanism,  if  it  occurs,  are 
greater  than  the  molecular- attraction  mechanism.  However,  liquid 
bridges  may  not  be  formed  under  the  experimental  conditions  of  inter- 
est since  negative  pressures  (tensile  forces)  within  the  liquid  are  pre- 
dicted for  a broad  range  of  particle  sizes  and  wetting  angles. 

Coating  materials  should  be  hygrophobic  and  should  have  nonpolar  molecular 
binding.  These  conditions  are  satisfied  by  several  fx)lymers,  of  which  such  fluo- 
rinated  hydrocarbons  as  teflon  are  notable  because  of  their  small  polarizabilities 
and  resulting  small  binding  forces.  Another  group  of  desirable  materials  is  the 
fatty-acid  group  including  stearic  acid.  These  acids  are  long,  linear  molecules 
having  one  reactive  end  and  an  inert  tail.  Under  favorable  reaction  conditions, 
these  acids  form  an  aligned  monomolecular  layer  on  metal  substrates,  with  the 
inert  nonpolar  tails  aligned  normal  to  the  surface. 

Li^t  scattering  from  metallic  spheres  coated  with  various  absorbing  and  non- 
absorbing coatings  has  been  examined  qualitatively  and  quantitatively.  The  present 
study  is  of  scattering  from  actual  metallic  spheres,  in  contrast  to  previous  studies 
which  considered  perfectly  conducting  spheres.  The  use  of  an  actual  metal  is  nec- 
essary because  the  effects  of  many  coatings  are  less  than  the  differences  between 
a metal  and  a perfect  conductor.  Means  of  overcoming  the  difficulties  that  arise 
in  the  Mie-theory  calculation  of  an  actual  metal  and  descriptions  of  the  calculation 
procedure  are  given  in  the  body  and  appendices  of  this  report.  The  principal  re- 
sults of  the  light-scattering  study  are: 

• The  scattering  efficiency  of  an  uncoated  metallic  sphere  is,  on  the  aver- 
age, within  three  percent  of  that  for  a perfectly  conducting  sphere  of  the 
same  radius.  In  the  wavelength  range  of  interest,  the  maximum  deviation 
is  less  than  ten  percent. 

• With  nondispersive  absorbing  and  transparent  coatings  having  refractive 
indices  greater  than  approximately  1. 6,  the  scattering  is  qualitatively 
similar  to  that  from  a perfectly  conducting  sphere  having  a radius  equal 
to  the  total  of  the  core  sphere  radius  plus  the  coating  thickness.  Witli 
smaller  magnitude  coating  refractive  indices,  tlie  scattering  is  better 
approximated  by  that  from  a perfectly  conducting  sphere  liaving  the  core 
radius. 
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• Quantitative  effects  of  nondispersive  coatings  were  evaluated  by 
calculating  the  root -mean -square  percent  (RMS%)  deviation  of  the 
coated  metallic  sphere  from  the  uncoated  metallic  sphere  over  a por- 
tion of  the  scattering  range.  With  thin  coatings,  the  RMS%  deviation  is 
linear  in  the  coating  thickness  at  a fixed  refractive  index.  For  coating 
refractive  indices  in  the  range  from  approximately  1.6  to  4.0,  the  RMS% 
deviation  is  insensitive  to  the  refractive  index  and  is  less  than  0. 3 % for 
a radius  ratio  of  the  coated  to  the  uncoated  sphere  of  1. 0004.  For  thicker 
coatings  having  radius  ratios  1.004  and  1.04,  the  RI^%  deviations  are 
less  than  3%  and  30%,  respectively. 

• Dispersive  absorption  bands  were  investigated  by  considering  a model 
coating  having  a single  Lorentzian  absorption  line  of  various  frequencies 
and  various  linewidths.  Ilie  absorption  lines  result  in  localized  deviations 
in  the  scattering  curve  centered  approximately  one  linewidth  toward  hi^er 
frequencies  from  the  absorption  line  and  having  a width  comparable  to  the 
absorption  linewidth.  The  shape  of  the  localized  deviation  is  primarily  de- 
pendent on  the  absorption  frequency  and  is  a localized  valley  near  the  first 
maximum  in  the  scattering  curve,  a localized  shoulder  near  the  midpoint 
between  the  first  maximum  and  first  minimum,  and  a localized  peak  near 
the  first  minimum.  For  an  absorption  coefficient  of  10^  cm'^  at  the  line 
center  and  a radius  ratio  of  1.04,  the  localized  deviations  have  a magni- 
tude of  approximately  10%  of  the  local  value  of  the  scattering. 

• Scattering  from  an  actual  metal-oxide  coating  was  examined.  The  localized 
deviations  in  the  scattering  curve  do  not  correlate  well  with  the  absorption 
bands  in  the  coating.  Rather,  the  most  severe  effects  occur  at  points  where 
the  real  part  of  the  refractive  index  is  small  and  the  imaginary  part  is  large. 

• Scattering  from  homogeneous  dielectric  particles  smaller  in  size  than  the 
primary  particles  was  examined.  On  the  average  the  scattering  is  less  effi- 
cient than  from  metallic  particles  of  the  same  size.  There  is  one  localized 
valley  in  the  scattering  curve  at  a point  where  the  real  part  of  the  refrac- 
tive index  is  slowly  varying  and  the  imaginary  part  changes  abruptly  from 
one  slowly  varying  value  to  another. 
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PREFACE 


This  final  report  covers  research  on  crystal  chemistry  studies  performed 
on  Contract  No.  F19628-75-C-0123  during  the  period  from  March  28,  1975 
through  June  30,  1976.  The  following  investigators  contributed  to  the  research: 
Dr.  M.  Sparks,  principal  investigator;  Dr.  C.  J.  Duthler,  research  scientist; 

Dr.  R.  Sherman,  research  associate;  and  Mr.  M,  Flannery,  research  associate. 


4 


TABLE  OF  CONTENTS 


Section  page 

SUMMARY  OF  RESULTS 1 

PREFACE 4 

I INTRODUCTION 13 

II  COHESION  MECHANISMS 14 

A.  Liquid  Bridge 14 

B.  Electrostatic 24 

C.  Molecular  Attraction 29 

D.  Summary  and  Conclusions  of  Cohesion  Study 37 

III  LIGHT  SCATTERING  FROM  COATED  SPHERES 42 

A.  Theory 42 

B.  Practical  Computational  Considerations 48 

C.  Tests  of  Program  Accuracy 55 

D.  Results 58 

E.  Summary  and  Conclusions  of  Scattering  Study 119 

REFERENCES 127 

Appendix  A DESCRIPTION  OF  THE  SUBROUTINES  AND  THEIR 

FUNCTION 129 

Appendix  B INPUT  DATA  CARD  DESCRIPTION 132 

Appendix  C LIST  OF  VARIABLES 136 

Appendix  D COMPUTER  CODE  LISTINGS 139 


5 


LIST  OF  ILLUSTRATIONS 

igure  Title  Page 


2. 1 Two  spheres  of  radius  a wet  by  a liquid  out  to  an  angle  0 from 

the  line  connecting  the  sphere  centers.  The  liquid  has  a con- 
tact angle  p with  the  sphere  surfaces  and  has  a radius  of 
curvature  r.  15 

2.2  Dependence  of  the  liquid  pressure  differential  on  wetting  angle.  18 

2.3  Dependence  of  the  attractive  force  components  on  wetting  angle.  20 

2.4  Dependence  of  the  critical  wetting  angle  on  sphere  radius.  21 

2. 5 Dependence  of  the  liquid-bridge  attractive  forces  on  sphere 

radius.  22 

2.6  A charge  of  magnitude  ne  at  a distance  of  2a  from  a conducting 

sphere  of  radius  a,  induces  an  image  charge  of  magnitude 

-ne/2  in  the  conducting  sphere.  27 

2.7  Dependence  of  the  molecular  attractive  force  on  coating 

thickness.  32 

2. 8 Separation  h of  the  sphere  surfaces  at  a distance  p from  the 

line  between  sphere  centers.  35 

2. 9 Comparison  of  cohesion  mechanisms.  38 

2. 10  Aligned  monolayer  of  fatty  acid.  41 

3. 1 Radiation  propagating  in  the  z direction  is  scattered  from  a 

spherically  symmetric  object  at  the  origm  to  the  observation 
point  (r,  e,p).  43 

3.2  Coated  sphere  having  core  radius,  a,  and  total  radius,  b.  47 

3.3  Magnitude  of  the  electric -multipole  scattering  coefficients  of 

a perfectly  conducting  sphere  as  a function  of  the  sphere  size 
parameter  ka.  49 

3.4  Magnitude  of  the  magnetic -multipole  scattering  coefficients  bj^  of 

a perfectly  conducting  sphere  as  a function  of  the  sphere  size 
parameter  ka.  50 

3.5  Comparison  of  the  scattering  from  a perfectly  conducting  sphere 

(symbols)  and  a metallic  sphere  (solid  curve).  60 


6 


LIST  OF  ILLUSTRATIONS  (Cont’d) 


Figure 


Title 


3.6  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols),  62 

3. 7 Scattering  from  a coated  metallic  sphere  normalized  to  the  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  63 

3.  8 Scattering  from  a coated  metallic  sphere  normalized  to  the  total 
cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols).  64 

3. 9 Scattering  from  a coated  metallic  sphere  normalized  to  the  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fecQy  conducting  sphere  of  radius  a (symbols).  65 

3. 10  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols).  66 

3. 11  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  67 

3. 12  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols).  68 

3. 13  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  69 

3. 14  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols).  70 

3. 15  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  7 1 

3. 16  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 
I cross-sectional  area  (solid  curve)  and  scattering  from  a per- 

" fcctly  conducting  sphere  of  radius  b (symbols).  72 


.(  1 


LIST  OF  ILLUSTRATIONS  (Cont’d) 


Figure  Title  Page 


3. 17  Scattering  from  a co"^ed  metallic  sphere  normalized  to  the  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  73 

3. 18  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols).  74 

3. 19  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  75 

3.20  Scattering  from  a coated  metallic  sphere  normalized  to  tlie  total 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols).  78 

3.21  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  79 

3.22  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols).  80 

3.23  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  81 

3.24  Scattering  from  a coated  metallic  sphere  normalized  to  tlie  total 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols).  82 

3.25  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  83 

3.26  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols).  84 

3. 27  Scattering  from  a coated  metallic  sphere  normalized  to  tlie  core 

cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  85 


8 


h 


LIS'l’  OF  ILLUSTRATIONS  (Cont’d) 


Scattering  from  a coated  metallic  sphere  normalized  to  the  total 
cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  b (symbols).  86 

Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a per- 
fectly conducting  sphere  of  radius  a (symbols).  87 

Root-mean-square  percent  deviation  of  the  scattering  of  a coated 
metallic  sphere  from  tr.at  of  the  uncoated  sphere  as  a function 
of  the  coating  refractive  index  for  various  coating  thicknesses.  90 

Contours  of  constant  root-mean-square  percent  deviation  of  the 
scattering  of  a coated  metallic  sphere  from  that  of  the  un- 
coated sphere  as  a function  of  the  coating  thickness  and  coating 
refractive  index.  92 

Real  part  C;{  and  imaginary  part  cj  of  the  dielectric  constant  of  a 
coating  having  a Lorentzian  absorption  line  with  a full  width  at 
half  maximum  of  50  cm"  ^ (a).  Scattering  from  a coated  metallic 
sphere  in  the  vicinity  of  the  absorption  line  (b).  94 

Real  part  and  imaginary  part  cj  of  the  dielectric  constant  of  a 
coating  having  a Lorentzian  absorption  line  with  a full  width  at 
half  maximum  of  20  cm"^  (a).  Scattering  from  a coated  metallic 
sphere  in  the  vicinity  of  the  absorption  line  (b).  95 

Real  part  and  imaginary  part  ej  of  the  dielectric  constant  of  a 
coating  having  a Lorentzian  absorption  line  with  a full  width  at 
half  maximum  of  10  cm"^(a).  Scattering  from  a coated  metallic 
sphere  in  the  vicinity  of  the  absorption  line  (b).  96 

Scattering  from  a coated  metallic  sphere  having  a 10  cm  ^ fwhm 
absorption  line  at  ka  = 1.01  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  97 

Scattering  from  a coated  metallic  sphere  having  a 10  cm  ^ fwhm 

• absorption  line  at  ka  = 1, 12  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  98 

Scattering  from  a coated  metallic  sphere  having  a 10  cm  ^ fwhm 
absorption  line  at  ka  = 1.26  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  99 


9 


LIST  OF  ILLUSTRATIONS  (Cont'd) 


Figure 

3.38 

3.  39 

3.40 

3.41 

3. 42 

3.43 

3.44 

3.45 

3.46 


Tide  Page 

Scattering  from  a coated  metallic  sphere  having  a 10  cm  ^ fwhm 
absorption  line  at  ka  = 1.43  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  100 

Scattering  from  a coated  metallic  sphere  having  a 10  cm  ^ fwhm 
absorption  line  at  ka  = 1.65  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  101 

Scattering  from  a corted  metallic  sphere  having  a 20  cm  ^ fwhm 
absorption  line  at  ka  = 1. 01  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  102 

Scattering  from  a coated  metallic  sphere  having  a 20  cm  ^ fwhm 
absorption  line  at  ka  = 1. 12  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  103 

Scattering  from  a coated  metallic  sphere  having  a 20  cm  ^ fwhm 
absorption  line  at  ka  = 1.26  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  104 

Scattering  from  a coated  metallic  sphere  having  a 20  cm  ^ fwhm 
absorption  line  at  ka  = 1.43  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  105 

Scattering  from  a coated  metallic  sphere  having  a 20  cm  ^ fwhm 
absorption  line  at  ka  = 1.65  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols^  106 

Scattering  from  a coated  metallic  sphere  having  a 50  cm  ^ fwhm 
absorption  line  at  ka  = 1.01  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  107 

Scattering  from  a coated  metallic  sphere  liaving  a 50  cm  ^ fwhm 
absorption  line  at  ka  = 1. 12  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  108 


10 


LIST  OF  ILLUSTRATIONS  (Cont'd) 


Figure  Title  Page 

3.47  Scattering  from  a coated  metallic  sphere  having  a 50  cm"^  fwhm 

absorption  line  at  ka  = 1.26  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  109 

3. 48  Scattering  from  a coated  metallic  sphere  having  a 50  cm”^  fwhm 

absorption  line  at  ka  = 1.43  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  110 

3.49  Scattering  from  a coated  metallic  sphere  having  a 50  cm'^  fwhm 

absorption  line  at  ka  = 1. 65  (solid  curve).  Results  are  normal- 
ized to  the  core  cross-sectional  area  and  compared  to  a per- 
fectly conducting  sphere  of  radius  b (symbols).  Ill 

3. 50  Scattering  from  a coated  metallic  sphere  having  a 10  cm’^  fwhm 

absorption  line  (solid  curve).  Results  are  normalized  to  the 
total  cross-sectional  area  and  compared  to  a perfectly  conduct- 
ing sphere  of  radius  b (symbols).  113 

3. 51  Scattering  from  a coated  metallic  sphere  having  a 10  cm'^  fwhm 

absorption  line  (solid  curve).  Results  are  normalized  to  the  core 
cross-sectional  area  and  compared  to  a perfectly  conducting 
sphere  of  radius  a (symbols).  114 

3.52  Real  part  n of  the  refractive  index  of  a metal-oxide  coating.  116 

3.53  Imaginary  part  k of  the  refractive  index  of  a metal-oxide  coating.  117 

3.54  Scattering  from  a metallic  sphere  coated  with  a metal-oxide 

coating.  H8 

3. 55  Real  part  n of  tlie  refractive  index  of  an  absorbing  dielectric 

material.  120 

3. 56  Imaginary  part  k of  the  refractive  index  of  an  absorbing  dielectric 

material.  12 1 

3.57  Scattering  from  a small,  homogeneous  sphere  of  absorbing  dielec - 

^ trie  material.  122 

B.  1 Typical  set  of  data  cards  for  the  coated  sphere  scattering 

program.  I35 


LIST  OF  TABLES 


[ I 

I 

Title  Page 


The  attractive  force  F and  components  of  tlie  attractive  force  F , 
from  the  pressure  differential,  and  F^,  from  the  direct  pull 
of  the  liquid,  for  several  sphere  radii,  a.  Tlie  critical  wetting 
angle  is  6q.  23 

Equilibrium  number  of  excess  or  deficient  electron  charges  (n), 
sphere  potential  $ , electric  field  E,  and  attractive  force  F for 
several  sphere  radii,  a.  F^a^  force  obtained  for  an  elec- 
tric field  of  3 xlO"^  V/cm  at  the  sphere  surface.  26 

Molecular  attractive  forces  for  spheres  of  various  radii,  a, 
coated  with  a nonpolar,  dielectric  coating.  34 

Description  of  program  variables.  137 

i 


1 

i 

i 

r 

f 


Table 

2.1 

I 

s 

2.2 

j 

[ 

i C.l 


12 


I.  INTRODUCTION 


4 

"1  ■ 

0«r  research  program  is  presently  concerned  with  light  scattering  from 
dispersed  metallic  spheres  that  have  been  coated  witli  a passive  layer  to  re- 
duce cohesion  and  to  prevent  excessive  material  degradation,  ht  Section  II, 
various  mechanisms  for  the  cohesion  of  spheres  are  investigated  and  the  con- 
ditions under  ^^dlich  particular  mechaiiisias  become  dominant  are  determined. 
Using  the  Lenard-Jones  potential  as  a mathematically  convenient  model  poten- 
tial for  the  molecular  attraction  of  spheres,  the  thickness  of  a molecular- solid 
coating  required  to  saturate  the  molecular  binding  from  strong  metallic  binding 
to  relatively  weak  van  der  Waals  binding  is  estimated. 

, the  theory  of  light  scattering  from  coated  spheres  is  inves- 
tigated. The  qualitative  relationship  between  actual  metallic  spheres  and  per- 
fectly conducting  spheres  is  explored  along  with  tlie  (jualitative  similarity  of 
coated  metallic  spheres  and  perfectly  conducting  spheres.  Quantitative  effects 
of  nondispersive,  absorbing,  and  nonabsorbing  coatings  of  various  thicknesses 
are  determined.  In  the  final  part  of  this  section,  Lorentzian  absorption  lines 
of  various  linewidths  are  used  as  models  for  absorption  and  dispersion  in  actual 
coatings.  A listing  of  the  computer  program  for  calculating  the  scattering  from 
coated  spheres  and  descriptions  of  the  subroutines  are  given  in  the  appendices. 


II.  COHESION  MECHANISMS 


There  are  three  mechanisms  that  can  contribute  to  the  cohesion  of  particles: 
a)  If  two  particles  are  linked  by  a liquid  bridge  that  wets  both  particles,  the  surface 
tension  of  the  liquid  results  in  an  attractive  force,  b)  An  electrically  charged  par- 
ticle experiences  attractive  forces  to  oppositely  charged  particles,  or  to  image 
charges  induced  in  other  uncharged  particles  and  other  uncharged  surfaces,  c)  Two 
particles  in  close  proximity  experience  molecular  attraction  and  molecular  binding. 
The  strength  and  range  of  the  molecular  attraction  depend  on  the  particle  material 
and  the  chemical-bonding  mechanism  appropriate  to  the  material:  i.e. , van  der 
Waals,  metallic,  covalent,  or  ionic, 

A.  Liquid  Bridge 

The  simple  theory  of  the  liquid-bridge  mechanism  for  particle  cohesion  pre- 
4-6 

sented  in  many  textbooks  seems  to  be  incorrect,  especially  for  very  small  size 
spheres.  A more  exact  theory  is  presented  in  the  following  and  the  shortcoming  of 
the  simple  theory  is  indicated.  Both  theories  yield  relatively  large  cohesive  forces 
implying  that  the  liquid-bridge  mechanism  would  be  the  dominant  mechanism  for 
particle  cohesion,  if  it  actually  occurs.  However  this  mechanism  seems  not  to 
a^ply  to  many  experimental  situations  because  the  state  of  the  adsorbed  liquid  is 
not  a thermodynamic  equilibrium  state  in  many  cases. 

Consider  tlie  contact  of  two  spheres,  each  having  radius  a,  as  is  sketched  in 
Fig.  2. 1.  An  adsorbed  liquid  wets  both  particles  out  to  an  angle  6 from  tlie  line 
connecting  the  sphere  centers.  The  conditions  for  the  liquid  to  wet  the  sphere  ma- 
terial are  discussed  in  Ref.  7.  Generally,  water  will  wet  metal  and  metal-oxide 
surfaces  but  will  not  wet  many  polymer  and  other  low  binding  energy  (hence,  low 
surface  energy)  materials. 

g 

The  attractive  force  between  spheres  can  be  written 

F = AP  trp^  + 2'rra'y  sui  6 sin  (0  + cp)  , (2. 1) 

2 

where  irp  is  the  cross-sectional  area  of  the  liquid,  y is  the  liquid  surface  tension, 
and  (p  is  the  contact  angle  that  the  liquid  makes  with  the  sphere  surface.  The  first 
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term  in  Eq.  (2. 1)  results  from  the  fact  that  the  pressure  within  the  liquid  is  lower 
than  the  exterior  by  an  amount  AP.  The  second  term  comes  from  the  direct  pull  of 
the  liquid  along  the  lines  of  contact  on  the  spheres.  With  a contact  angle  of  zero, 
Eq.  (2. 1)  becomes 

2 2 

F = APffp  + 2fraysin  6 . (2.2) 

The  pressure  differential  across  the  liquid  surface  is  inversely  proportional 
to  the  principal  radii  of  curvature  of  the  liquid  surface,  i.  e. , the  radius  of  curva- 
ture at  the  exterior  surface  and  the  radius  of  curvature  around  the  axis  of  symme- 

For  the  concave  exterior  surface  shown  in  Fig,  2, 1,  the  actual  shape  of  the 
liquid  surface  is  a nodoid.  At  large  wetting  angles  the  liquid  surface  becomes  a 
catenoid  or  unduloid  corresponding,  respectively,  to  a zero  or  a positive  pressure 
differential  across  the  liquid  surface.  The  solution  for  the  principal  radii  of  cur- 
vature of  the  actual  nodoidal  surface  involves  an  iterative  solution  of  equations  in- 
volving incomplete  elliptical  integrals.  For  simplicity  in  the  following,  the  surface 
is  approximated  by  a circular  arc  that  is  tangent  to  the  spheres  at  the  wetting  angle. 
With  this  approximation,  the  principal  radii  of  curvature  are  obtained  from  simple 
geometry  which  yields 

r = ( ^ ~ cos  9) 


_ / cos  9 + sin0  - 1 ' 


The  pressure  differential  across  the  liquid  surface  is 


Derivations  of  the  liquid-bridge  force  presented  in  several  textbooks'^ as- 
sume a circular-arc  liquid  surface,  neglect  the  second  term  in  Eqs.  (2. 1)  and 
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(2.2)  and  neglect  the  second  term  in  Eq.  (2.5).  These  approximations  are 
mathematically  valid  at  small  wetting  angles  and  yield 


= 2TTay 


as  B ->  0 


Although  the  attractive  force  in  Eq.  (2.6)  remains  finite  as  the  wetting  angle  ap- 
proaches zero,  the  pressure  differential  approaches  negative  infinity.  This  implies 
that  the  liquid  is  under  considerable  tension  at  small  wetting  angles.  This  does 
not  seem  physically  possible,  although  it  is  not  discussed  in  the  textbooks  present- 
ing this  derivation. 

At  finite  wetting  angles,  both  the  attractive  force  and  the  pressure  differential 
are  angle  dependent.  To  determine  the  dependence  of  the  pressure  differential  on 
the  wetting  angle,  a dimensionless  pressure  is  derived  by  substituting  the  radii  of 
curvature  from  Eqs.  (2.  3)  and  (2.4)  into  Eq.  (2.5)  and  factoring  a term  y/a  : 


(a/y)AP  = 


1 - cos  6 


cos  6 

COS0  + sin0  - 1 


This  dependence  of  the  pressure  differential  on  angle  is  plotted  in  Fig.  2.2.  At 

O 

small  angles  the  pressure  differential  goes  to  infinity  as  B . At  6 = 53. 1°  the 
pressure  differential  vanishes  because  the  two  principal  radii  of  curvature  are 
equal  and  opposite. 

Similarly,  the  two  force  terms  in  Eq.  (2.2)  can  be  written  in  dimensionless 
form  to  bring  out  the  angular  dependence  of  these  terms.  The  first  component  of 
the  attractive  force  p-p,  from  the  pressure  differential,  becomes 


F /2iray  = (cos9+  sinB-  1)  1 1 

P 2cos0  1-COS0  cos0+sin0-l 


The  other  term  F^,  from  the  direct  pull  of  the  liquid,  is 


F^/2iray  = sm  0 . 
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These  forces  are  plotted  as  functions  of  the  wetting  angle  in  Fig.  2.3.  At  small 
angles,  the  pressure  term  dominates  and  approaches  the  limiting  value  given  by 
Eq.  (2.6).  As  the  wetting  angle  increases,  contributions  from  become  non- 
negligible.  The  two  terms  are  equal  at  33.5°  and  the  pressure  term  vanishes 
at  53.1°. 

To  solve  for  the  attractive  force  ar  a function  of  sphere  size,  the  following 

approach  is  taken:  First,  a "critical  wettmg  angle”  is  determined  at  which 

^ ^ 

the  pressure  differential  equals  atmospheric  pressure  (1.01  X lO'^  dync/cm'^) 

2 

using  the  surface  tension  of  water  y - 70  dyne/cm  . Next  6q  is  substituted  into 
Eqs.  (2.3)  and  (2.4)  which  along  with  Eqs.  (2.5)  and  (2.2)  yield  the  attractive 
force.  The  critical  wetting  angle  6q  is  plotted  as  a function  of  the  sphere  radii 
in  Fig.  2.4.  The  two  components  of  the  attractive  force  and  the  total  attractive 
force  F = F^  + F^  are  plotted  in  Fig.  2.5.  Tabular  results  are  presented  in 
Table  2.1. 

The  above  derivation  makes  the  reasonable  and  convenient  approximations  of 
a zero  contact  angle  and  a circular  arc  surface.  In  addition,  the  sagging  of  the 
liquid  due  to  gravity  and  the  buoyant  force  of  the  liquid  are  neglected,  as  is  valid 
at  small  sphere  sizes. 

It  is  seen  from  Table  2. 1 and  Fig.  2.  4 that  the  critical  wetting  angles  for 
small  spheres  are  quite  large,  approaching  the  limiting  value  of  53. 1°.  in  this 
size  range,  the  attractive  force  presented  in  Fig.  2.5  is  dominated  by  tlie  direct 
surface  tension  term  F_^.  The  two  components  of  the  attractive  force  are  equal 
near  the  sphere  radius  of  2pm  (6^  = 33.5°)  and  the  pressure  differential  term 
Fp  dominates  at  sphere  radii  larger  than  this  value.  It  is  interesting  to  note  that 
the  total  force  F = F^  + F^  is  approximately  equal  to  the  naive  textbook  value 
given  in  Eq.  (2.6).  This  agreement  is  regarded  as  accidental. 

Should  the  particles  be  wet  to  angles  greater  than  0q  , the  cohesive  forces 
decrease  as  is  apparent  from  Fig.  2.2.  Larger  forces  than  those  presented  in 
Table  2. 1 and  Fig.  2.5  would  be  obtained  at  smaller  angles  if  large  negative 
pressures  (tensile  forces)  are  allowed  in  the  liquid.  There  is  no  discussion  of 
such  tensile  forces  in  the  literature  on  sphere  cohesion  except  for  a brief  state- 
ment in  Ref.  9a  where  it  is  stated  that  this  may  be  a limiting  effect.  A similar 
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Fable  2.1.  The  attractive  force  F and  components  of  the  attractive  force  , 
from  the  pressure  differential,  and  F^ , from  the  direct  pull  of  the  liquid,  for 
several  sphere  radii  a . The  critical  wetting  angle  is  6^  . 


a(|im) 

0g (degrees) 

Fp  (dyne) 

F^(dyne) 

F (dyne) 

0.05 

52. 1 

3.40x10'^ 

1.37x10"^ 

1.40x10"^ 

0.25 

48.4 

7.67  xlO'^ 

6. 15  xlO’^ 

6.92  xlO’^ 

0.50 

44.9 

2.70x10'^ 

l.lOxlO'^ 

1.37  xlO'^ 

2.5 

31.2 

3.80x10"^ 

2.95  xlO'^ 

6.75x10'^ 

5.0 

24.6 

1.02  xlO'^ 

3.81 X 10'^ 

1.40x10'^ 

25.0 

12.5 

7.74  xlO’^ 

5.15x10'^ 

8.26x10'^ 

50.0 

9.07 

1.71 

5.46x10'^ 

1.76 

. 


effect  occurs  with  very  small  capillary  tubes  wdicrc  negative  pressures  have 
been  observed. However  in  the  capillary  tubes,  the  tensile  state  of  the 
liquid  is  a metastable  state  and  is  lost  when  a bubble  is  nucleated.  Nevertheless 
such  metastablc  states  occur  regularly  m nature  as  is  evidenced  by  sap  in  the 
xylem  tubes  of  trees  that  are  over  10  m (34  ft. ) in  height.  Hence  if  wetting 
angles  less  than  the  critical  angle  6q  (x:cur,  the  liquid  will  be  under  tension  and 
in  a metastable  state. 

It  may  be  very  difficult  to  form  a liquid  bridge  between  particles  even  at 
angles  greater  than  0q  unless  one  starts  with  the  particles  completely  im- 
mersed in  the  liquid,  if  the  particles  are  initially  dry,  liquid  may  not  condense 
at  the  points  of  contact.  The  difficulty  in  condensation  here  is  analogous  to  the 
well-known  difficulty  in  the  homogeneous  nuclcation  of  small  water  drops  or 
small  bubbles.  It  may  be  possible  to  form  liquid  bridges  between  the  particles 
if  the  particles  are  initially  immersed  and  liquid  is  withdrawn  from  the  sample. 

In  this  case  sudden  evaporation  of  the  liquid  may  occur  when  the  critical  angle 
is  reached. 

B.  Electrostatic 

The  electrostatic  attraction  of  charged  spheres  is  complicated  by  the  facts 

that  the  mechanisms  for  charging  particles  are  not  well  understood  and  the  loca- 
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tions  of  the  charges  on,  or  within,  the  particles  generally  are  not  known.  ’ 

Only  with  conducting  spheres  that  are  free  of  contamination  can  predictions  be 
made  with  any  degree  of  confidence.  As  a known  example,  consider  the  contact 
of  spheres  of  two  different  metals  having  differing  electron  work  functions.  Elec- 
trons arc  transferred  from  one  sphere  to  the  otlier  until  the  Fermi  levels  arc 
aligned  (i.e. , as  in  a thermocouple).^'^  Excess  electrons  on  a conducting  sphere 
spatially  isolated  from  other  material  surfaces,  are  uniformly  distributed  over 
the  surface.  However,  die  presence  of  other  particles,  charged  or  not,  will  per- 
turb this  uniform  distribution. 

With  dielectric  particles  or  contaminated  metallic  particles,  charging  often 
involves  the  transfer  of  ions  absorbed  on  the  surfaces  as  well  as  the  transfer  of 
electrons.  The  ion-transfer  mcclianism  is  very  dependent  on  complexities  of  the 
surface  such  as  past  surface  treatments  and  the  relative  humidity  of  the  surround- 
ing atmosphere.  To  predict  the  charge  transfer  Ixitwccn  bulk  samples  of  various 
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materials,  researchers  have  attempted  to  set  up  triboelectric  series  of  materials 
analogous  to  electrochemical  potential  series.  However,  such  series  are  not  re- 
producible, with  series  set  up  by  different  researchers  being  very  dissimilar  and 
lacking  agreement  on  the  positions  of  materials  in  the  series.^^ 

Ignoring  the  mechanism  for  charging,  the  equilibrium  amount  of  charge  on  an 
isolated  sphere  can  be  derived  from  thermodynamic  considerations.^  A spatially 
isolated  sphere,  of  radius  a with  an  excess  or  deficiency  of  n electrons  uniformly 
distributed  over  its  surface,  has  an  electrostatic  potential  energy 

2 

e = (ne)  /2  a , (2.10) 

where  e is  the  charge  of  an  electron.  Assuming  a Boltzmann  distribution,  the 
probability  that  a sphere  has  n excess  or  deficient  electrons  is 

1 “ r -i 

f = exp  -(nc)^/2akT  / XI  exp  -(ne)^/2akT  , (2.11) 

J n = -00  [_  J 

where  k is  the  Boltzmaxm  constant  and  T is  the  temperature.  The  mean  magni- 
tude of  n is 

<n)  = (2akT/ite^)^'^^  . (2.12) 

Values  of  <n)  arc  listed  in  the  second  column  of  Table  2.2  for  several  sphere 
radii. 

Rather  than  attempting  to  specify  the  details  of  tlie  distribution  of  the  charge 
on  or  within  the  particle,  the  following  model  problem  is  considered:  With  a spa- 
tially isolated  and  uniformly  charged  sphere,  the  electric  field  outside  of  the 
sphere  appears  as  though  the  charge  is  concentrated  at  tlie  center  of  the  sphere. 
Hence  as  a first  approximation  for  a charged  sphere  in  contact  with  an  uncharged 
sphere,  the  excess  charge  is  considered  to  be  positioned  at  the  center  of  the 
charged  sphere.  If  the  spheres  both  have  radii  equal  to  a,  the  charge  ne  in  the 
first  sphere  induces  an  image  charge  of  magnitude  - ne/2  at  a distance  a/2  from 
the  center  of  the  second  sphere.  This  situation  is  sketched  in  Fig.  2. 6.  ( The  mag- 
nitude and  position  of  the  image  charge  is  appropriate  for  a conducting  sphere,  with 
dielectric  spheres  having  a smaller  magnitude  image  charge  at  a different  position. ) 
The  attractive  force  between  the  charges  is 
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Table  2.2.  Equilibrium  number  of  excess  or  deficient  electron  charges  <n), 
sphere  potential  electric  field  E , and  attractive  force  F for  several  sphere 
radii,  a.  is  the  force  obtained  for  an  electric  field  of  3 x V/cm  at  the 

sphere  surface. 


a(/im) 

<n) 

0>(mV) 

E(V/cm) 

F^^(lO'^^dyne) 

F (dyne) 

0.05 

0.  68 

19.6 

3910 

947 

5.57x10'^ 

0.25 

1.66 

9.55 

382 

226 

1.39x10"^ 

0.50 

2.37 

6.82 

136 

l^c 

5.57x10’^ 

2.5 

5.34 

3.07 

12.3 

23.4 

1.39xl0’^ 

5.0 

7.54 

2.17 

4.34 

11.6 

5.57x10"'^ 

25. 0 

16.9 

0.91 

0.39 

2.34 

1.39x10"^ 

50.0 

23.9 

0.  69 

0.14 

1.17 

5.57x10'^ 
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F = I ^ . (2.13) 

This  model  is  perhaps  most  valid  for  conducting  spheres  coated  with  a thin  dielec- 
tric layer.  For  conducting  spheres  where  an  insulating  layer  prevents  equilibration 
of  the  excess  charge,  the  result  in  Eq.  (2. 13)  is  the  first  term  in  a series  gener- 
ated by  alternately  placing  image  charges  in  the  two  spheres  to  satisfy  boundary 
conditions. 

The  force  given  by  Eq.  (2.8)  should  be  valid  as  an  order-of-magnitude  esti- 
mate for  many  seemingly  different  cases.  Notice  that  if  the  spheres  have  equal 
and  opposite  charges  positioned  at  the  sphere  centers,  the  factor  of  2/9  in  Eq. 

(2. 13)  is  replaced  witli  unity.  Larger  forces  than  those  indicated  by  Eq.  (2. 13) 
may  be  obtained  with  dielectric  spheres  having  oppositely  charged  surface  ions 
near  the  point  of  contact.  However,  the  ions  themselves  experience  large  attrac- 
tive forces  in  this  case  and  the  large  forces  may  be  ameliorated  by  the  ions  being 
torn  from  the  surfaces. 

Using  the  equilibrium  amount  of  charge  in  Eq.  (2. 13),  attractive  forces  are 
given  for  several  sphere  radii  in  the  third  column  of  Table  2.2.  These  attractive 
forces  are  many  orders  of  magnitude  smaller  tlian  either  those  given  in  the  previ- 
ous section  for  the  liquid-bridge  mechanism  or  those  derived  in  the  following  sec- 
tion for  molecular  attraction.  Hence  for  electrostatic  forces  to  be  important,  the 

magnitude  of  the  charge  must  be  much  greater  than  the  equilibrium  value.  ; 

i 

Under  favorable  circumstances  such  as  in  a corona  discharge,  spheres  can  j 

be  charged  far  in  excess  of  their  equilibrium  values  where  as  many  as  one  thou-  j 

sand  excess  electrons  arc  observed  on  micron-sized  particles.^’  To  search  ] 

for  a mechanism  to  limit  tiie  ultimate  possible  excess  charge,  the  electrical  poten-  ) 

tial  of  the  particle  and  the  magnitude  of  the  electrical  field  at  the  particle  surface  \ 

were  examined.  Values  of  these  quantities  are  listed  in  the  fourth  and  fifth  columns  | 

of  Table  2.2  for  spheres  having  equilibrium  numbers  of  electrons.  The  equilibrium  \ 

value  of  the  electrical  potential  is  in  the  millivolt  range  for  all  particle  sizes  and,  | 

hence,  is  not  expected  to  prevent  increases  in  the  charge  under  favorable  circum-  j 

stances.  On  the  other  hand,  thc>  electric  field  at  the  sphere  surface  is  quite  large  1 

for  small  spheres.  Attempts  to  increase  the  charge  on  small  spheres  beyond  the  j 

equilibrium  value  will  be  limited  by  air  breakdown  near  the  sphere  surface.  I 


r 
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To  estimate  the  maximum  possible  electrostatic  forces,  it  is  assumed  that 
the  sphere  charge  is  limited  by  a value  of  3 x 10“^  V/cm  for  the  electric  field  at 
the  sphere  surface.  Solving  for  the  charge  required  to  produce  this  field  and  sub- 
stituting the  result  into  Eq.  (2. 13)  yields  the  values  in  the  sixth  column  of  Table 
2.2  labeled  These  values  are  still  approximately  one  order  of  magnitude 

smaller  than  the  liquid-bridge  forces  obtained  in  the  previous  section.  However, 
at  the  largest  particle  sizes,  the  maximum  electrostatic  forces  become  compar- 
able to  molecular  forces. 

Such  highly  charged  particles  are  expected  to  be  the  exception  rather  than  the 
general  rule.  In  a large  sample  containing  many  particles,  binding  from  the  elec- 
trostatic mechanism  should  occur  only  for  relatively  few  particles, 

C.  Molecular  Attraction 


The  molecular  attraction  between  spheres  is  viewed  two  ways:^*^  First,  a 
model  potential  such  as  the  Lenard-Jones  potential  is  integrated  over  the  volume 
of  the  spheres.  Differentiation  yields  the  cohesive  force.  Secondly,  an  effec- 
tive area  of  contact  is  derived  from  geometrical  considerations.  This  effective 
area,  \\hen  multiplied  by  the  material  tensile  strength,  yields  the  force  of  attrac- 
tion. The  Lenard-Jones  potential  used  in  the  first  method  is  strictly  applicable 
only  to  nonpolar  molecular  solids  having  van  der  Waals  bonding.  However,  it  can 
be  used  as  a mathematically  convenient  model  potential  for  order-of-magnitude 
estimates  of  the  cohesive  force  between  other  types  of  materials.  Additionally, 
the  Lenard-Jones  potential  is  used  to  estimate  the  coating  thickness  required  to 
saturate  the  attraction  so  that  the  cohesion  is  determined  by  the  coating  rather 
than  the  core. 

Consider  two  spheres,  one  of  radius  a and  the  other  of  radius  b,  having  a 
distance  L between  sphere  centers.  The  Lenard-Jones  potential  between  volume 
elements  dVj  in  the  first  sphere  and  d^^  in  the  second  sphere  is 


d^U 


(2. 14) 


where  p is  the  distance  between  the  volume  elements  and  /3  is  a material -dependent 
constant.  Integration  over  the  volume  of  the  first  sphere  yields 
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d^u  = -4ffa^^dV2/3  , 


-4ira^^dV^/3  (r^ 


(2. 15) 


where  r is  the  distance  of  the  volume  element  d\^  from  the  center  of  the  first 
sphere.  Integration  over  the  volume  of  the  second  sphere  yields  the  potential 
between  the  two  spheres: 


( (L-b) 

I a^[(L-b)2-a2] 


(L+b)  1 

a^[(L+b)2-a^  2 


jln 


-L  J^n  ( ) - 

2 ' L-i-b+a  / 


+ 


L[(L-b)2-a2  L[(L+b)2-a2  [ . (2.16) 


Differentiating  with  respect  to  L and  collecting  terms  yields  the  attractive  force 
between  spheres: 


F = 


^2  2„ 

T ■’  ^ 


a b L 


[(L^  + b^-  g2)2  - 4b2  l2]2 


(2. 17) 


The  force  given  in  E,q.  (2. 17)  diverges  at  L - a + b because  we  have  neglected  to 
include  a repulsive  term  in  the  potential.  Rather  than  including  a repulsive  term, 
the  divergence  is  avoided  by  defining  sphere  contact  as  occurring  when  L = a + b 
+ 6,  where  6 is  a lattice  spacing. 

The  result  given  in  Eq.  (2. 17)  is  next  used  to  calculate  the  force  between  two 
coated  spheres  each  having  core  radius  a and  total  radius  b = a + 1,  where  t is 
the  coating  thickness.  For  coated  spheres  the  attractive  force  is  the  sum  of  the 
forces  between  cores,  between  cores  and  coatings,  and  between  coatings.  The 
force  between  cores  is 

'i.a  ’ T Al 

where  is  the  potential  parameter  for  the  core  material.  The  force  between 
each  core  and  the  coating  on  the  other  sphere  has  potential  parameter  This 
force  is  equal  to  the  force  lietwecn  the  core  and  a homogeneous,  solid  sphere  of 
the  coating  material  witli  radius  b minus  the  force  between  the  core  and  a sphere 
of  coating  material  having  radius  a : 


P — ^ ^ Q ^ 

a,b-  3 '’l2[(L2  + b2-a2)2-4b2L2j2  ' T * ^12  ' 

(2. 19) 

The  force  between  coatings  equals  the  force  between  homogeneous,  solid  spheres 
of  the  coating  material  of  radius  b minus  the  force  between  cores  of  coating  ma- 
terial and  minus  the  force  between  cores  and  coatings: 


3u3  T 
a b L 


64  2 


32  2 


64  2 


3k3t 
a b L 


TT  = _ « ft ^ aoL 

b,b  3 ^22^^3f^2,,^2j2  3 ^22  j ~2^'^2  . ,2,2 . 4,^l2  ]2 


^ 32  „2« 

^ 3 ^22  ,3rr2_.  2-,2 


(2.20) 


L-"[L^-4a^] 

Summing  Eqs.  (2. 18),  (2. 19),  and  (2.20)  yields  the  force  between  coated  spheres: 


T7  = 12  ^2  / ^ o ft  . o N a°  , 32  _2  „ b° 

- 3 " '^12^-  %2)  +-’'^22 


, 64  2 _ a^b^L 

3 12  22  [-(l2^j^2_  g2)2  . 4Jj2l2j2 


(2.21) 


To  examine  the  reduction  in  the  attractive  force  from  the  coating,  consider  a 
coatmg  having  potential  parameter  = 0*  j . The  relative  potential  parame- 
ter is  ^12=  0 11^22  = 0.3 16/3  To  obtain  the  effect  of  the  coating  thick- 

ness t,  we  take  b = a + t and  L = 2a  + 2t  + 6 witli  6 = 4 x 10"^  cm.  The 
attractive  force  between  two  coated  spheres  using  these  values  in  Eq.  (2.21)  is 
plotted  as  a function  of  coating  thickness  in  Fig.  2.7.  The  attractive  force  has  a 
minimum  value  for  coating  thicknesses  near  40  A.  For  thicknesses  greater  than 
this  value,  the  attractive  force  increases  due  to  the  increased  sizes  of  the  spheres. 
The  continuum  model  used  in  this  calculation  is  not  valid  for  coating  thicknesses 
less  dian  about  10  A.  Nevertheless,  the  rapid  decrease  of  the  attractive  force  for 
small  coating  tliickncsscs  obtained  from  the  continuum  model  indicates  that  only 
a few  molecular  layers  of  coating  are  necessary  to  saturate  the  forces  so  that  the 
force  is  determined  nearly  entirely  by  the  coating  rather  than  the  core. 
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To  obtain  numerical  values  for  metallic  spheres  coated  with  nonpolar 

dielectrics,  we  take  ^22  ~ where  Cl  is  the  atomic  volume  of  the  atoms 

in  the  coating  and  C.  is  the  constant  in  the  London  force  equation  for  the  force 
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between  two  spherically  symmetric  atoms.  ’ This  constant  is  evaluated  from 


= 3a  hp/4 


(2.22) 


where  a is  the  atomic  polarizability  and  hy  is  the  electronic  excitation  energy  of 
the  atoms.  As  typical  values  we  take  Cl  = {4  xlO’^cm)^,  a = lO’^^cm^,  and 
hp  = 1.60x10  ^^erg(lOcV).  These  values  yield 

^22  " 2.9x10"^^  erg  . (2.23) 


Assuming  that  the  coating  is  sufficiently  thick  to  saturate  the  attractive  force,  this 
value,  along  with  L = 2 b t 6,  where  6 = 4 x 10  cm  for  spheres  in  contact,  yields 
the  attractive  force 

F = ( 1.  49  dyne/cm  ) a . (2.24) 

Values  of  the  force  for  several  sphere  radii  are  presented  in  Table  2.  3. 

A second  way  of  viewing  the  molecular  attraction  between  two  spheres  in  con- 
tact is  from  a materials  strength  point  of  view.^’"^  In  this  method  the  spheres  are 
again  viewed  as  being  homogeneous  and  continuous,  rather  than  atomic.  The  hypo- 
thesis is  made  that  the  distinction  of  one  sphere  surface  from  the  other  is  lost  if 

the  surfaces  are  within  a distance  of  approximately  one  lattice  spacing,  i.e.,  about 
”8 

4 X 10  cm.  Referring  to  Fig.  2.  8,  it  is  seen  that  the  distance  h between  sphere 
surfaces  is  approximately  h s;  p^/ a,  wiiere  p is  the  distance  from  the  line  con- 
necting the  sphere  centers  and  a is  the  sphere  radius.  The  effective  area  of  con- 
tact between  spheres  is 

A = Ttp^  = 7rah  . (2.25) 

Multiplying  this  area  by  the  tensile  strength  of  the  material  yields  the  cohesive 
force 

F = o^^Trah  . (2.26) 
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Table  2.3.  Molecular  attractive  forces  for  spheres  of  various  radii,  a,  coated 
with  a nonpolar,  dielectric  coating. 


a ( ixm  ) 0. 05 

F(dyne)  7.4x10^ 


0.25  0.50 

3.7x10’^  7.4x10'^ 


2.5  5.0 

3.7xl0‘^  7.4x10'^ 


25.0  50.0 

3.7x10'^  7.4x10'^ 
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For  nonpolar  molecular  solids  the  two  methods  yield  equivalent  results. 

The  attractive  force,  derived  from  the  Lenard-Jones  potential  presented  in 
Eq,  (2.24),  implies  a tensile  strength  = 1.2  x 10^  dyne/cm^  (170  psi).  This 
value  is  approximately  what  one  would  expect  for  a molecular  solid.  The  tensile 
strength  of  metals  is  approximately  100  times  this  value  and  implies  an  approxi- 
mately 100  times  greater  attractive  force. 

Intuitively  the  cohesive  force  between  oxidized- metal  spheres  is  expected  to 
be  less  than  the  force  for  pure  metals.  With  metal  oxides  or  with  other  ionic  and 
covalent  solids,  estimates  of  the  attractive  force  are  complicated  by  uncertain- 
ties in  the  tensile  strength  of  these  materials.  Tensile  strengths  of  bulk  ionic  and 
covalent  solids  are  sensitive  to  small  cracks  wtiich  are  comparable  in  size  to  the 
smaUest  sphere  sizes  under  consideration,  if  cracks  are  not  important  with  small 
spheres  and  theoretical  tensile  strengths  are  used,  the  cohesive  force  is  as  large, 
or  greater  than,  the  cohesive  force  for  pure  metals.  Instead,  it  is  expected  that 
dielectric  coatings  are  sensitive  to  the  crystallographic  orientation  of  the  coatings 
on  the  two  spheres.  This  misalignment  of  the  lattices  may  result  in  a decreased 
cohesive  force  analogous  to  decreased  tensile  strengths  in  bulk,  dielectric  solids 
from  small  cracks. 

Additionally,  malleability  may  favor  the  cohesion  of  pure  metals  due  to  in- 
creases in  the  effective  area  of  contact  caused  by  deformation  around  the  point  of 
contact  under  the  action  of  molecular  attraction.^  When  two  elastic  spheres  are 
pressed  together  with  a force  F,  there  is  a flat  area  surrounding  the  point  of  con- 
tact having  radius 

p = 0.88  [Fa/E  , (2.27) 

where  E is  the  elasticity  of  the  sphere  material.  Using  the  attractive  force  given 
by  Eq.  (2.26)  and  using  E - lO^^dyne/cm^  and  = lO^dyne/cm^  as  typical 
values  for  a metal,  the  area  of  contact  is 

A = (2.84x10'^)  a^"'^  , (2^28) 

where  the  units  of  A arc  cm^  for  a expressed  in  cm.  With  a = 2.5  x lo"^  'm,  Eq. 
(2.28)  yields  4.  5 x 10  ^^cm^  which  is  comparable  to  the  effective  area  of  3.  1 
X 10  ^^cm^  given  by  Eq.  (2. 15).  These  values  indicate  that  the  effective  'rea  of 
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contact  for  a pure  metal,  or  for  a thin  oxidc-coated  metal,  may  be  enhanced 
compared  to  dielectrics  and  thick  oxide  coatings.  The  continuum  model  is  not 
thought  to  be  sufficiently  accurate  to  estimate  the  amount  of  enhancement. 


D.  Summary  and  Conclusions  of  Cohesion  Study 

Cohesive  forces  from  the  various  mechanisms  discussed  in  the  preceding 
sections  are  presented  graphically  in  Fig.  2.9.  The  attraction  from  the  liquid- 
bridge  mechanism,  if  it  occurs,  has  tlte  largest  magnitude.  However,  liquid 
bridges  may  not  be  fonned  between  particles  under  the  experimental  conditions 
of  interest,  as  is  discussed  in  Section  11  A.  Molecular  attractive  forces  are 
greatest  for  metallic  spheres  and  smallest  for  spheres  coated  with  a nonpolar 
dielectric  having  van  der  Waals  molecular  binding.  Coatings  having  covalent  or 
ionic  binding  are  expected  to  have  intermediate  force  values  with  the  range  of 
possible  molecular  attraction  forces  extending  between  (and  perhaps  somewhat 
beyond)  the  two  curves  labeled  molecular,  metal  and  molecular,  Lenard-Jones. 

j To  avoid  the  large  cohesive  forces  between  metallic  spheres,  the  spheres 

I should  be  coated  with  a nonpolar  dielectric  coating  having  van  der  Waals  binding. 

The  calculation  in  the  previous  section  indicates  tiiat  the  coating  need  only  be  a 
few  molecular  layers  tliick  to  saturate  the  forces  so  that  the  attraction  is  deter- 
mined by  the  coating  rather  than  tlie  core.  Tlie  liquid-binding  mechanism  can  be 
■ avoided  if  tlie  coating  is  hygrophobic. 

Tlie  concern  that  electrostatic  forces  may  become  troublesome  with  dielec- 
tric coatings  docs  not  appear  to  be  a problem.  If  tlie  particles  have  the  thermo- 
dynamic equilibrium  amount  of  charge,  the  electrostatic  forces  are  negligible, 
as  is  indicated  by  the  bottom  curve  of  Fig.  2.9,  The  amount  of  charge  required 
to  make  the  electrostatic  force  comparable  to  the  molecular  force  is  unreason- 
ably large,  with  large  electric  fields  occurring  near  the  particle  surface.  It  is 
thought  that  such  large  charges  will  only  occur  in  exceptional  circumstances 
and,  in  general,  will  be  present  only  on  a few  spheres  of  a large  sample. 

One  of  tlie  easiest  coatings  to  apply  to  the  particles  is  an  oxide  coating. 

Even  though  the  tensile  strength  of  oxides  can  be  as  large  as  metals,  the  sphere 
binding  force  may  be  smaller  for  oxidized  spheres  lx;cause  of  the  sensitivity  of 
oxides  to  crystallographic  orientation.  If  the  spheres  are  oxidized  wliile  in  contact. 
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an  oxide  bridge  of  large  cross  section  may  be  grown  between  the  spheres  with  a 
subsequently  large  binding  force.  It  will  then  be  necessary  to  mechanically  break 
the  spheres  apart.  Once  separated,  smaller  attractive  forces  are  expected  to  be 
obtained. 

Another  simple  method  of  reducing  cohesion  is  to  use  spacer  particles  smaller 
in  size  than  the  primary  particles  to  avoid  contact  of  the  primary  particles.  The 
spacer  particles  can  be  made  of  a material  such  as  carbon  that  has  weak  binding 
to  the  metal.  Additionally  because  of  the  smaller  size  of  the  spacer  particles,  the 
attractive  force  between  a spacer  and  primary  particles  is  smaller  than  between 
two  particles  of  the  same  size  so  that  the  cohesion  of  the  mixture  is  reduced. 
Spacer  particles  are  expected  to  have  little  effect  on  the  light  scattering  from 
the  dispersed  mixture  for  two  reasons;  First,  spacer  particles  are  only  a small 
fraction  of  the  total  mass.  Second,  the  size  of  the  spacer  particles  falls  in  a range 
where  the  scattering  is  inefficient,  i.e. , the  scattering  cross  section  is  smaller 
than  the  geometrical  cross  section.  Should  a spacer  particle  remain  bound  on  the 
surface  of  a primary  metallic  particle,  it  should  have  little  effect  on  the  scatter- 
ing of  the  primary  particle.  The  maximum  effect  on  the  scattering  will  occur  if 
the  spacer  is  positioned  at  a local  maximum  in  tlie  electric  field  on  the  surface  of 
the  primary  particle.  It  could  tlien  result  in  damping  of  one  of  the  multipole  mo- 
ments contributing  to  the  scattering,  but  the  effect  should  be  small  since  several 
multipole  moments  contribute. 

Coatings  that  have  desirable  properties  from  a molecular  point  of  view  are 
many  polymers  and  such  fatty  acids  as  stearic  acid.  Of  the  polymers,  teflon  is 
desirable  because  it  is  hygrophobic  and  because  it,  along  with  other  fluorinated 
hydrocarbons,  has  a small  polarizability  which  yields  small  attractive  forces. 
Methods  for  depositing  teflon  from  the  gas  phase  by  ultraviolet  photolysis  have 
been  published.^^'^^  A possible  problem  with  polymer  coatings  that  could  lead 
to  increasing  binding  is  incomplete  chain  termination  witli  loose  bonds  on  two 
particles  combining  and  forming  a bridge  between  particles. 

Low  cohesive  forces  should  be  obtained  with  metallic  spheres  coated  with  an 
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aligned  monomolecular  layer  of  fatty  acid.  ' ’ Stearic,  lauric,  and  other  fatty 
acids  are  used  as  additives  in  lubricating  oils  to  provide  boundary  lubrication 
under  high-pressure  conditions,  vacuum  conditions,  or  other  situations  vbere  the 
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lubricating  oil  may  be  lost.  The  fatty  acids  consist  of  a backbone  of  carbon 
atoms  with  a reactive  carboxyl  radical  on  one  end.  When  applied  from  a dilute 
solution  to  active  metals  such  as  copper  or  zinc,  the  carboxyl  radical  reacts 
with  the  substrate  metal  forming  an  aligned  monomolecular  layer  as  is  sketched 
in  Fig.  2. 10.  Such  a layer  will  not  be  formed  on  noble  metals  such  as  gold  and 
silver.  It  has  been  noted  that  care  is  needed  to  get  a reaction  on  aluminum,  with 
factors  such  as  the  degree  of  surface  strain  and  amount  of  oxide  present  on  the 
surface  mfluencing  the  reaction.  If  two  surfaces  coated  with  monolayers  of  fatty 
acid  are  brought  into  contact,  contact  is  made  between  the  inert,  nonpolar  tails 
facing  outward  from  the  coated  metals.  The  Lenard-Jones  potential  is  applicable 
for  the  interaction  between  the  tails  and  very  small  forces  of  attraction  should  be 
obtained.  Such  coatings  can  possibly  be  applied  to  small  metallic  spheres  by  re- 
action from  a volatile  solvent  such  as  benzine,  although  with  semi -inert  metals 
such  as  aluminum  more  care  may  be  needed  to  insure  a reaction  with  the  sub- 
strate. 


[ 

} 

I 
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Pig.  2.  10.  Aligned  monolayer  of  fatty  acid. 
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III.  LIGHT  SCATTERING  FROM  COATED  SPHERES 


A.  Theory 


Scattering  from  a spherically  symmetric  object  is  described  with  the  aid  of 

o o_oc 

Fig.  3. 1.  Incident  radiation  of  unit  intensity  and  polarized  in  the  x direc- 

tion propagates  in  the  z direction  to  the  particle  located  at  the  origin.  The  scat- 
tered intensities  at  the  position  (r,  0,<p)  for  the  two  polarizations  in  the  6 and 
Ip  directions  are 


L - \ -j  Is,  (6)  1^  cos^cp 

r 


^10 IT  I ^2^®^  I 

^ 4ff  r 


where  X is  the  wavelength.  The  amplitude  functions  Sj(6)  and  S2(0)  are 


S (8)  = L ^ 

■t  n=l  n 


2 n+  1 


a^7Tn(cos0)  + b^T^  (cos0) 


= „li  jfSr 


anT^(cos0)  + b^7Tj^(cos0) 


where  ir^  and  are  Legendre  functions  and  where  and  b^  are  Mie  scattering 
coefficients  which  are  determined  by  boundary  conditions  on  the  particle  surfaces. 

The  Legendre  functions  can  be  written  in  terms  of  associated  Legendre  poly- 
nomials of  order  n and  degree  unity: 

ff  (cos0)  = P^^^(cos0)/sin0  (3.5) 


T (COS0)  = dp5,^^(cos0)/d0  . 
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The  Legendre  functions  have  simple  values  for  scattering  in  the  forward  or 
backward  direction.  For  these  directions,  the  values  are 


7T^(0°)  = T^(0°)  = n(n+  l)/2 

in  the  forward  direction  and 

1T^(180°)  = -Tj^(180°)  = (-1)^  n(n  + l)/2 


(3.  8a) 


(3. 8b) 


in  the  backward  direction.  The  differential  scattering  cross  section  per  unit  solid 
angle  is 


o(6,(p)  = 4i7r  1(6, <p) 


(3.9) 


The  differential  scattering  efficiency  is  the  ratio  of  the  scattering  cross  section  to 
the  geometrical  cross  section  of  the  sphere.  Dividing  Eq.  (3. 9)  by  n & for  a sphere 
of  radius  a and  using  the  above  values  for  the  Legendre  functions  in  Eq.  (3. 8)  yields 
the  differential  scattering  efficiencies  (or  gains)  in  the  forward  and  backward  direc- 
tions: 


and 


G(0“)  = (4/a^)  I E n+  i (a  + b ) 1^ 

I n=l  2 n n'  I 

G(180°)  = (4/a^)  I E n + i (-l)'’(b-a)r 
I n=l  2 ' ' n n ' I 

where  a is  the  sphere  size  parameter 
OL  = 2ira/X 


(3. 10) 


(3.11) 


(3. 12) 
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Similarly,  simple  expressions  are  obtained  for  the  integrated  efficiencies 

for  scattering,  Q , and  extinction,  Q : 

SC  Si  OXu 


<^sca  - + lb„l\ 


= (2/oc^)  E (2n  + l)  Re(a  + b ) , (3.14) 

exi  n=l  “ ^ 

where  Re  denotes  the  real  part.  The  extinction  efficiency  gives  the  total  energy 
lost  from  the  incident  beam  to  scattering  and  absorption.  The  absorption  effi- 
ciency is  obtained  by  taking  the  difference  between  the  extinction  efficiency  and 
scattering  efficiency. 

It  remains  to  calculate  the  Mie  scattering  coefficients,  a^  and  b^.  These  co- 
efficients are  associated  with  the  various  multipole  moments  induced  in  the  sphere 
by  the  incident  radiation  with  representing  the  electric  dipole  moment,  at2 
representing  the  electric  quadrupole  moment,  bj  representing  the  magnetic  dipole 
moment,  etc. 

For  a homogeneous  sphere  having  a complex  refractive  index  m = n - i k , 
the  Mie  scattering  coefficients  arc 

i/j’ (ma)i^  (a)  - mb  ( ma ) (a ) 

a = — (3.15) 

" ^ ma ) Cj^(a  ) - ( ma ) (a ) 

and 

mi;;' (ma)b  (a)  - b (ma)^’ (a) 

b = 2 ^ 2 2 , (3.16) 

n nnpj^(ma)Cn(a)  - Vn 

where  primes  denote  derivatives  with  respect  to  the  arguments  of  the  functions. 
The  function  ^ is  a linear  combination  of  functions: 

c = il)  + iX  , (3.17) 


where  b„  and  are  Riccati-Bessel  functions  of  the  first  and  second  kind, 
n 
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In  the  limit  of  an  infinite  refractive  index  ( perfectly  conducting  sphere),  the 
above  expressions  for  tlie  Mie  scattering  coefficients  reduce  to 

3n  = (3.18) 

and 

bn  = ^n(a)/Cn(a)  . (3.19) 

Notice  that  these  equations,  also  the  various  scattering  efficiencies,  involve  only 
the  particle  size  parameter  a = 2 it  a / X.  Hence  parameters  for  perfectly  con- 
ducting spheres  provide  limiting  values  for  testing  results  for  spheres  having 
finite  optical  constants. 

The  Mie  scattering  coefficients  for  a coated  sphere  are  more  complex  than 

those  of  a homogeneous  sphere  because  of  the  necessity  of  satisfying  boundary 

conditions  at  tlie  coating-core  interface  as  well  as  at  the  outer  surface.  A coated 

sphere  having  core  radius  a and  total  radius  b is  depicted  in  Fig.  3.2.  Optical 

constants  in  the  core  and  coating  are  denoted  by  subscripts  1 and  2,  respectively. 

23 

Expressions  for  the  Mie  scattering  coefficients  given  in  Kerker  are  incorrect. 
The  corrected  coefficients  are  given  by  tlie  determinantal  expressions: 


m2<^;(mi0!) 

0 

X^im.^a) 

0 

li;(m2^) 

X;(ni2t^) 

0 

m2V;(t/) 

X^(m2l^) 

0 

m^X;(ra2«) 

ni2V  ;(m^a) 

0 

Vj^On^a) 

0 

0 

m,c;(.) 

0 

C,(.) 
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rr 


! i 


and 


b = 
n 


|mii^;(m^a) 

0 

m^Xn(>"2“> 

0 

v;(m2i^) 

\;(m2t^) 

0 

0 

m^X^(m.,0!) 

<i)  ^(m  jOf) 

0 

iQ-) 

0 

0 

c;(F) 

0 

m2Cn(>^) 

(3.21) 


where  a = 2ira/X  and  y = 2iTb/X. 

B.  Practical  Computational  Considerations 

The  calculation  of  tlie  scattering  efficiency  of  a sphere  involves  straightfor- 
ward computation  of  the  Mie  scattering  coefficients  from  Eqs.  (3.20)  and  (3.21) 
for  a coated  sphere  [ or  Eqs.  (3. 15)  and  (3. 16)  for  a homogeneous  sphere]  and 
substitution  of  these  coefficients  in  the  desired  scattering  efficiency  sum  in  Eqs. 

(3. 10),  (3. 11),  (3. 13),  or  (3. 14).  Factors  requiring  consideration  are:  the  number 
of  terms  required  in  the  Mie  scries,  an  efficient  and  accurate  means  of  calculating 
the  Riccati-Bcsscl  functions,  and  modifications  to  the  formulas  to  handle  overflows 
that  occur  with  the  large  values  of  the  refractive  index  for  metals. 

The  number  of  terms  required  in  the  Mie  series  is  dependent  primarily  on  the 
ratio  of  the  sphere  radius  to  the  incident  radiation  wavelength  and  to  a lesser  ex- 
tent on  the  sphere  refractive  index.  To  determine  the  number,  the  magnitudes  of 
the  first  eleven  scattermg  coefficients  for  a perfectly  conducting  sphere  of  radius  b 
were  calculated,  'fhe  logarithms  of  these  terms  are  plotted  as  a function  of  kb 
= 2'rrb/  X in  Figs.  3.3  and  3.4.  Prior  to  tlieir  initial  approach  to  unit  magnitude. 
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Fig.  3.3.  Magnitude  of  the  electric-multipole  scattering  coefficients  of  a 
perfectly  conducting  sphere  as  a function  of  the  sphere  size  parameter  ka. 
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the  scattering  coefficients  are  seen  to  increase  with  increasing  kb  and  to  decrease 
with  increasing  order.  lii  the  following  calculations,  the  Mic  series  has  been  trun- 
cated at  ten  terms.  For  sizes  and  wavelengths  in  the  range  of  Figs.  3.3  and  3.4, 
the  truncation  error  in  the  Mie  series  is  of  the  order  of  a^^  as  is  indi- 

cated by  the  dashed  curves. 

The  above  estimate  of  the  truncation  error  should  also  be  approximately  valid 
for  dielectric-coated  metallic  spheres  in  tliis  size  and  wavelength  region.  An  ex- 
amination of  the  scattering  coefficients  of  a few  dielectric -coated  metallic  spheres 
reveals  that,  even  with  thick  coatings  of  b/a  = 1.4,  the  major  differences  from  a 
perfectly  conducting  sphere  occur  in  the  region  near  magnitude  unity.  In  the  region 
where  the  scattering  coefficients  are  small,  which  determines  tlie  truncation  error, 
the  coefficients  are  quite  similar. 

Most  of  the  numerical  subtleties  in  the  program  are  contained  in  the  algo- 
rithms for  calculating  the  lliccati- Bessel  functions  and  their  derivatives.  Rather 
than  using  explicit  algorithms  for  calculating  all  orders  of  the  Riccati-Bc‘ssei  func- 
f tions  and  tlieir  derivatives,  two  of  the  orders  are  calculated  explicitly  with  recur- 

I sion  relations  being  used  to  generate  the  remaining  oraers  and  all  derivatives. 

' Routines  are  needed  only  for  the  psi  and  chi  functions  since  the  zeta  function  is  a 

linear  combination  of  tiiese  functions,  as  is  given  m Eq.  (3. 17). 

Recursion  relations  for  are  stalile  in  the  forward  direction,  so  and 
are  evaluated  using 

\^(z)  = sinz  + cosz/z  (3.22) 

and 

\^{z)  = -cosz  + 3\^(z)/z  , (3.23) 

except  for  arguments  having  large  imaginary  parts  as  is  discussed  later.  Use  of 
the  recursion  relation 

^n+1  = (2n  + l)X^/z  - X^,^  (3.24) 

yields  orders  two  through  ten. 
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Recursion  relations  for  are  unstable  in  the  forward  direction,  so  the 
backward  recursion  relation 


*^n-l  ~ (2n  + l);/)^/z  (3.25) 

is  used  with  starting  functions  and  to  generate  orders  one  through  eight. 
Derivatives  are  obtained  using 


’^n 


z ^n 


for  orders  one  through  nine  and 
*^n  = ^n-1  - 


(3.26) 


(3.27) 


for  order  ten.  These  same  derivative  equations  are  also  valid  for  when  X is 
substituted  for  ijj . 

Algorithms  involving  sines  and  cosines  multiplied  by  polynomials  in  z ^ can 
be  written  for  the  starting  functions  ij)g  and  , However  in  contrast  to  Eqs. 
(3.22)  and  (3.23)  for  Xj  and  X2 « the  analogous  algorithms  for  and  are 
inaccurate  at  small  arguments,  in  addition  to  having  difficulties  at  large  imaginary 
arguments  as  discussed  later.  For  arguments  greater  than  about  twice  the  sub- 
script, but  not  having  too  large  an  imaginary  part,  tiie  most  accurate  algorithms 
are: 


9 


and 


45  - 13860/z^  + 945945/z'^ 

- 1621 6200 /z^  + 344  5942  5/z^ 
- 1 - 990 /z^  i-  135135/z'^ 


sin  z/z 


- 4729725/z^  + 34459425/z^ 


cos  z , 


z I >20.0 


(3.28) 


IJJ 
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- ] - 1485/z^  + 315315^“^  - 18918900/z^ 

t 3J0134825/z^  - 654729075/z^^^ j sin  z 

- 55  - 25740/z^  + 2837835/z'^ 

- 91891800/z^'  + 654729()75/z^|  cos  z/z  , |z|>20.0  • (3.29) 
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At  small  /. , tho  terms  in  il\e  /.  ^ polynomials  Ivcome  very  large  and 
accuracy  is  lost  in  taking  the  difference  of  nearly  equal  large  numbers.  In  this 
region  it  is  more  accurate  to  exjiand  the  sines  and  cosines  in  power  series  which 
yield  the  over-all  power  series: 


n + 1 


.„(z) 


1*  3*  5-  • • (2n  + 1) 


1 - 


1 2 
2" 


1!  (2n  + 3) 


2!  (2n  -r  3)  (2n+  5) 


|z|  < 20.0  . 

(3.30) 


Tlie  power  series  in  Eq,  (3. 30)  is  not  used  for  all  magnitudes  z Lncludmg  those 
greater  than  20  because  it  is  inefficient  and  inaccurate  due  to  the  large  number  of 
terms  required  in  the  scries  with  large  arguments.  The  cross-over  point  in  accu- 
racy between  the  two  metliods  was  determined  by  computing  and  both  ways 
for  values  of  z from  jzl  = 10  to  [z  j = 30  and  locating  the  point  where  they  agreed 
to  greatest  number  of  decimals.  On  the  IBM 370  computer,  run  in  double  precision, 
the  results  of  the  two  methods  agree  to  six  places  in  the  region  around  the  cross- 
over point.  This  six-significant-figure  limitation  near  the  cross-over  point  is  the 
major  limiting  factor  in  the  accuracy  of  the  scattering  calculation  in  the  region. 

As  mentioned  above,  there  are  additional  difficulties  in  both  Riccati-Bessel 
functions  for  arguments  having  large  imaginary  parts,  as  occurs  for  those  con- 
taining the  refractive  index  of  a metal.  In  tliis  case  the  sine  and  cosine  factors 
in  the  Lliccati- Bessel  function  algoritluns  may  have  magnitudes  greater  than  the 
capacity  of  the  computer.  To  avoid  this  difficulty,  a large  positive  exponential 
term  from  the  imaginary  part  of  the  sine  and  cosine  arguments  is  factored  out. 

Consider  first  the  modifications  required  for  a metallic  core  sphere  having 
a complex  refractive  index  m^.  Using  the  convention  wliere  the  imaginary  part  of 
the  refractive  index  is  negative,  the  function  argument  involving  the  core  refrac- 
tive index  is  wi'itten 


m^a  = u - 1 V 


where  botli  u and  v arc  positive. 

, , f i(u-iv) 

cos  ( m^a ) = I - ' 


The  cosine  of  this  term  is 


(3.31) 


(3.32) 
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I’actorlng  our  the  term  having  the  divergent  exponential  part  yields 


i 


] 

t 


I 

1 

i 

1 


cos  ( m^a ) 


(v+  iu ) 
e 


, . -2v  -2  iu 
1 + e c 


Similarly  the  sine  of  this  argument  is  written 


(3.  33) 


sin  ( ) 


-2v 

e 


(3.34) 


Notice  that  the  same  divergent  factor  appears  in  the  first  square-bracketed 
term  of  both  Eqs.  (3. 33)  and  (3.  34).  Since  the  Riccati-Bessel  functions  are  linear 
in  sine  and  cosine,  this  same  term  can  be  factored  from  these  functions,  and  their 
derivatives,  leaving  the  second  square-bracketed  terms  in  place  of  the  sine  and  co- 
sine. Rt'fcrring  to  F,c|s.  (3.20)  and  (3.21)  for  tlie  Mie  scattering  coefficients,  the 
argument  m^a  appears  only  iji  the  tliird  columns  of  each  determinant.  By  the  rules 
of  matrix  algebra,  the  divergent  term  can  also  be  factored  from  the  determinants. 
When  this  is  done,  terms  factored  from  the  numerator  and  denominator  cancel. 
Thus,  the  divergent  term  is  entirely  compensated  by  replacing  the  sines  and  co- 
sines in  the  Riccati-Bessel  functions  of  tlie  core  sphere  by  the  second  square- 
bracketed  terms  in  Eqs.  (3.33)  and  (3.  34).  No  other  changes  are  necessary. 


Additionally,  it  is  desirable  to  expand  the  capabilities  of  the  computer  program 
to  handle  metal-clad  spheres  as  well  as  dielectric -coated  spheres.  In  this  case  it 
is  necessary  to  remove  divergences  occurring  in  Riccati-Bessel  functions  having 
arguments  0120:  and  m^u.  For  convenience,  only  a real  factor  c^/2  is  removed 
from  tile  algorithms,  for  example 


r " 

e V2 

r m -ill  -V  -] 

cos(m,jO!)  = 

.1  a a 

e f e e 

- * 

Functions  involving  the  coating  refractive  index  are  found  in  the  first  two  columns 
of  the  determinants  in  Eqs.  (3.20)  and  (3.21).  in  contrast  to  the  core  functions, 
there  is  no  common  factor  in  each  column  because  both  arguments  m20:  and 
appear  in  each  column.  However,  when  the  determinants  are  expanded,  there  are 
two  over-all  common  factors  1/2  exii(v^)  and  1/2  exp(vjj)  from  two  arguments. 
'Ilicse  common  factors  are  made  possible  by  the  zeros  in  tlie  third  and  fourth 
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columns.  As  with  the  core  functions,  the  divergent  factors  from  die  numerators 
and  denominators  of  the  scattering  coefficients  cancel,  and  it  is  only  necessary 
to  replace  the  shies  and  cosines  in  tlie  Kiccati-lkissel  functions  by  terms  similar 
to  the  second  square-bracketed  term  in  Eq.  (3.35). 

With  the  coathig  functions  it  is  also  necessary  to  remove  the  same  factor 
from  the  power  scries  algorithm  used  for  small  arguments  of  the  psi  function 
in  order  to  keep  functions  computed  from  this  algorithm  normalized  with  re- 
spect to  the  other  coating  functions.  Otherwise  erroneous  results  occur  in  the 
region  where  | m^o:|  < 20.0  or  lm.,p(  < 20.0.  Such  a difficulty  does  not  occur 
with  the  psi  function  of  the  core. 

Besides  tlie  difficulties  with  overflows  in  the  algorithms,  m some  cases  un- 
derflows occurred  in  products  of  functions  in  tlie  expanded  determinant.  These 
products  were  examined  luid  were  found  to  be  negligible  with  respect  to  other 
terms  in  the  expanded  determinants.  Hence,  the  underflow  difficulty  was  over- 
come, with  no  loss  in  accuracy,  by  testing  for  an  underflow  and  by  setting  under- 
flowing products  equal  to  zero. 

A listing  of  the  computer  program  is  given  in  the  appendices  along  with  a 
further  description  of  the  implementation  of  the  above  equations. 

C.  Tests  of  Program  Accuracy 

Four  series  of  tests  were  made  of  the  numerical  accuracy  of  the  program: 

1)  'Phe  Riccati-Bessel  function  subroutines  were  compared  with  published  tables. 

2)  A few  points  on  the  computer-calculated  scattering  cuiwes  were  compared  in 
detail  to  a calculation  on  a hand  calculator.  3)  It  was  verified  that  tlie  coated- 
spliere  formulas  reproduced  tlie  homogeneous  sphere  results  in  tlie  proper  limits. 
4)  Scattering  results  of  the  coated-sphere  calculation  were  compared  with  previ- 
ously published  curves, 

1.  To  test  the  accuracy  of  tlie  Riccati-Bessel  function  subroutines,  both 
algorithms  used  for  small  and  large  arguments  were  used  to  generate  and 
V|Q  for  a large  number  of  pure  real  and  pure  imaginary  arguments.  Recursion 
relations  were  used  to  generate  through  Since  tables  of  Riccati-Bessel 


functions  were  not  available,  these  functions  were  converted  to  spherical  Bessel 
functions  and  using 

^^(z)  = zjj^(z)  and  = zy^(z)  . (3.36) 

Values  of  spherical  Bessel  functions  are  available  in  Tables  10.1,  10.2,  and  10.5 
of  Abramowitz  and  Stegun.^^  Less  extensive  tests  were  made  of  the  functions 
since  only  a single  algorithm  is  needed  for  these  functions. 

For  large  z the  expressions  for  and  are  quite  complex  and  the  avail- 
able tables  are  rather  sparse  (only  z = 20,  20  i,  50,  50  i,  100,  100  i ).  Additional 
tests  were  made  of  these  fiuictions  l^y  generating  ^d  using  back- recursion 
relations  from  calculated  values  of  ipg  and  niesc  values  of  j/jj  and  ^2  were 
compared  to  other  values  calculated  directly  from  their  analytic  expressions. 
Values  of  and  ip2  calculated  by  both  methods  agreed  to  six  decimal  places  near 
1 z 1 = 20  with  improvement  to  ten  places  for  (z  j < 10  or  j z ] > 30.  This  is  the 
same  accuracy  that  is  expected  for  the  starting  functions  and  Hence  be- 
sides testing  the  starting  function,  these  results  indicate  that  no  accuracy  is  lost 
in  the  use  of  the  recursion  relations. 

2.  At  a few  selected  points  on  the  scattering  curves,  the  complete  calculation 
was  carefully  paralleled  step  by  step  by  a tedious  hand  calculation.  The  agreement 
of  these  calculations  gives  us  great  confidence  in  the  over-all  computational  pro- 
cedure. In  cases  where  products  of  Bessel  functions  produced  underflows  in  die 
computer,  this  calculation  revealed  that  there  were  always  additional  non-ncgli- 
gible  terms  that  summed  widi  the  underflowing  terms  so  that  there  is  no  loss  of 
accuracy  in  setting  the  underflowing  parts  to  zero. 

3.  As  an  additional  test  of  the  over-all  program,  it  was  verified  that  the 
coated-sphcrc  program  reproduced  the  homogeneous  sphere  results  in  the  ap- 
propriate limits.  Besides  providing  a test  of  the  formulas  for  the  scattering 
coefficients,  this  procedure  also  tests  die  Riccati-Bessel  function  subroutines 
since  the  coated  sphere  coefficients  involve  more  complex  products  than  die 
homogeneous  sphere  coefficients.  The  limit,  b finite  and  a = 0,  is  not  a proper 
homogeneous  sphere  limit  of  the  coated-sphere  formulas.  When  the  electromagne- 
tic boundary  value  problem  was  solved  at  r = a,  solutions  in  the  coating  involve 
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both  psi  and  chi  functions,  but  in  the  core  sphere  only  psi  functions  are  allowed 
since  chi  functions  are  singiQar  at  the  origin.  Hence  the  limit  a = 0 cannot  be 
taken  since  singular  chi  functions  would  be  introduced  into  the  region  containing 
the  origin. 

The  homogeneous  sphere  limit  was  taken  in  the  following  five  cases: 

• metal  core,  vacuum  coating,  b = a 

• metal  core,  arbitrary  coatings,  b = a 

• metal  core,  vacuum  coating,  b ^ a 

• metal  core,  same  metal  coating,  b ^ a 

• dielectric  core,  same  dielectric  coating,  b^  a . 

In  all  limits  the  results  agreed  within  the  six- significant-figure  accuracy  of  the 
program. 

4.  "Riere  are  no  previous  calculations  of  scattering  from  a dielectric-coated 
metallic  sphere  available  for  comparison  to  our  results.  Previous  mvestigators 
have  approximated  homogeneous  metallic  spheres  by  a perfectly  conducting  sphere. 
Scharfman  has  studied  dielectric -coated  perfectly  conducting  spheres  in  the  region 
near  ka  = 1.0.^^  In  the  course  of  our  computer  calculation,  scattering  from  per- 
fectly conducting  spheres  of  radii  a and  b was  calculated  for  comparison  to  the 
coated-sphere  results.  Our  calculations  for  a perfectly  conducting  sphere  agree 
with  our  own  hand  calculations.  Fig.  4.27  of  Kerker,^^  and  Scharfman.^^  Our 

97 

coated-sphere  results  presented  in  die  following  section  agree  witli  Scharfman 
with  the  differences  being  witliin  the  differences  of  a perfectly  conducting  sphere 
and  an  actual  metallic  sphere. 

For  comparison  to  a sphere  with  finite  optical  constants,  a dielectric  sphere 
coated  with  the  same  dielectric  was  compared  to  die  homogeneous  dielectric 
sphere  curve  in  Fig.  4.34  of  Kerker^'^  using  the  following  parameters: 

mf  = m2  = 1.29  - il.37  , b/a  = 1.4 

There  was  good  agreement  between  our  calculated  results  and  the  curve  in 
Kerker.^'^ 
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D.  Results 


Previous  investigators  have  approximated  scattering  from  metallic  spheres 
by  results  for  a perfect  conductor.  Scattering  from  coated  perfect  conductors  has 
been  considered  by  Scharfman  and  by  Rheinstein.  By  making  this  approxima- 
tion, the  difficulties  with  large  imaginary  arguments  of  the  Rice ati- Bessel  functions, 
as  occur  with  actual  metals,  are  avoided.  As  is  shown  in  the  following,  scattering 
from  a homogeneous,  metallic  sphere  is  well  approximated  by  scattering  from  a 
perfect  conductor.  However  for  examination  for  the  effects  of  dielectric  coatings, 
it  is  necessary  to  consider  actual  metals  since  changes  in  the  scattering  for  many 
coatings  are  less  than  the  differences  between  an  actual  metal  and  a perfect  con- 
ductor. 

The  results  of  four  studies  are  presented  in  this  section.  First,  scattering 
from  a homogeneous,  metallic  sphere  is  compared  to  that  from  a perfect  con- 
ductor. Second,  a qualitative  comparison  is  made  of  scattering  from  dielectric - 
coated  metals  and  dielectric -coated  perfect  conductors.  The  similarity  of  the 
scattering  of  both  types  of  coated  spheres  to  scattering  of  a homogeneous,  perfect 
conductor  is  determined.  Third,  differences  in  the  scattering  of  metallic  spheres 
due  to  various  nondispersive,  dielectric  coatings  is  investigated  quantitatively. 
Fourth,  the  effects  of  dispersion  in  the  dielectric  coatings  are  investigated  by  con- 
sidering model  coatings  having  a single  Lx) rentzian- shaped  absorption  line.  Scat- 
tering from  an  actual  metal  oxide  coating  and  from  small,  homogeneous  dielectric 
spheres  is  also  calculated. 

1.  Homogeneous,  metallic  spheres. 

The  infinite  refractive  index  of  a perfect  conductor  is  independent  of  the  wave- 
length. Hence,  scattering  curves  which  are  plotted  as  a function  of  the  particle  size 
parameter  a = ka  = 27Ta/X  can  be  interpreted  as  being  functions  of  changes  in 
either  the  sphere  radius  or  in  the  wavelength.  This  duality  does  not  occur  with  ac- 
tual metals  because  of  dispersion  in  the  optical  constants.  Metallic  spheres  con- 
sidered in  this  section  have  fixed  radii  and  a complex  refractive  index 

m2  (Xj  microns)  = -12. 5 + 3.  8 X - i ( 9.  9 + 5. 68X  ) . (3.37) 
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The  linear  dependence  of  the  refractive  index  on  wavelength  in  Eq.  (3.37)  is  a 
good  fit  to  handbook  values  over  tlie  wavelength  range  of  interest.  In  this  region 
the  reflectivity  of  tlie  metal  is  at  least  0.  978,  so  that  the  perfect  conductor  is  ex- 
pected to  be  a good  approximation. 

Scattering  from  such  a metal  is  compared  to  a perfect  conductor  of  the  same 
radius  b in  Fig.  3.5,  where  scattering  from  the  metal  is  denoted  by  the  solid 
curve  and  the  perfect-conductor  results  are  given  by  -f  symbols.  The  qualitative 
agreement  is  good  with  the  oscillations  in  tlie  two  curves  being  in  phase  through- 
out the  kb  range  from  zero  to  four  given  in  the  figure.  The  magnitudes  of  the 
scattering  agree  in  the  region  from  kb  equal  to  zero  to  the  first  peak.  After  the 
first  peak  the  amplitude  of  tlie  oscillations  in  the  metal  curve  are  greater  than 
those  in  the  perfect-conductor  cuive.  At  the  two  minima,  tlie  two  scattering 
curves  differ  by  approximately  ten  percent.  '^Tlie  mean  difference  in  the  region 
between  the  first  two  peaks  is  four  percent. 

There  appears  to  be  no  simple  method  for  estimating  the  quantitative  dif- 
ferences between  the  scattering  from  an  actual  metal  and  the  scattering  from  a 
perfect  conductor  without  doing  the  difficult  Mie-tiieory  calculation.  Notice  that 
the  approximately  ten  percent  difference  in  the  two  curves  at  the  minima  is  much 
greater  than  the  approximately  two  percent  difference  in  the  bulk  reflectivities. 

At  other  points,  the  two  cuives  cross  and  tliere  is  no  difference  in  the  scattering 
at  these  crossing  points. 

2.  Qualitative  results  for  coated  metallic  spheres. 

27 

Scharfman  has  investigated  the  scattering  efficiency  in  the  ka  region  around 
the  first  peak  of  perfectly  conducting  spheres  coated  with  nondispersive  dielectrics. 
For  high-index  coatings  widi  refractive  indices  greater  than  approximately  2.5,  his 
scattering  curves  have  a first  peak  at  nearly  the  same  position  and  of  nearly  the 
same  magnitude  as  a perfectly  conducting  sphere  of  tlie  same  total  radius,  b.  This 
result  is  somewhat  surprising  in  view  of  tlie  low  reflectivity  of  tlie  bulk  dielectric 
material.  Another  limit  is  the  limit  of  a vacuum  coating  (refractive  index  equal 
unity)  where  the  scattering  must  approach  the  results  for  an  luicoatcd  sphere. 

This  limit  is  approximately  valid  for  thin  coatings  having  refractive  uidices  less 
than  1. 6. 
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Metallic  splK'Tcs  coated  witli  nondispersive  absorbijig  ;ind  transparent 
dielectrics  were  examijicd  to  find  whether  similar  qualitative  rules  are  appli- 
cable in  this  case,  and  whether  the  rules  are  valid  at  other  maxima  and  minima 
besides  the  first  peak.  This  study  was  done  graphically  by  preparing  two  plots 
of  the  scattering  efficiency  of  tlie  coated  sphere  which  has  core  radius  a and 

total  radius  b:  In  the  first  plot  the  scattering  of  the  coated  metallic  sphere  is 

2 

normalized  to  the  total  geometrical  cross  section,  Ttb  , and  plotted  as  a solid 
curve  as  a function  of  kb.  For  comparison,  the  scattering  efficiency  of  a per- 
fectly conductmg  sphere  of  radius  b is  presented  as  a series  of  + symbols  on 
the  same  axes.  In  the  second  plot,  scattering  from  tlie  coated  metallic  sphere 
is  normalized  to  the  core  cross  section,  rra^,  and  plotted  versus  ka . Compar- 
ison in  this  case  is  made  to  a perfectly  conductmg  sphere  of  radius  a . 

For  rlie  investigation  of  nondispersive,  transparent  coatings,  a two-dimen- 
sional grid  of  coating  tliicknesses  and  refractive  mdices  was  used.  At  each  of  the 
three  radius  ratios,  b/a  = 1.04,  1.004,  1.0004,  the  scattering  was  calculated  for 

the  set  of  refractive  indices  m„  = 1.2,  1.35,  1.5,  1.75,  2.5,  3.25,  and  4.0.  Two 

^ 2 2 

plots  of  the  scattering  efficiency,  normalized  to  TTa  and  Ttb  , are  presented  in 
Figs.  3.6  to  3. 19  for  the  full  set  of  refractive  indices  and  the  single  radius  ratio 
b/a  = 1.04. 

Consider  first  the  scattering  curves  with  b/a  = 1.04  plotted  in  terms  of  the 
larger  radius,  i.  e. , kb  abscissa.  For  the  coating  refractive  index  m2  =4.0 
shown  in  Fig.  3.6,  the  curves  for  the  coated,  metallic  sphere  and  the  perfectly 
conducting  sphere  are  m -phase,  having  maxima  £md  minima  at  nearly  the  same 
positions.  Olie  curves  differ  in  amplitude  witli  the  magnitude  of  the  coated,  me- 
tallic sphere  ciirvi-  Ix'ing  gn-ater  than  the  pei'fectly  conducting  curve  at  tlie 
maxima  and  smaller  at  the  minima,  'llic  differences  in  amplitude  increase 
with  increasing  kb. 

As  the  coating  refractive  index  decreases,  the  amplitude  of  the  coated, 
metallic  sphere  scattering  curve  decreases  and  the  period  increases.  With  a 
coating  refractive  index  of  1.5  shown  in  Fig.  3. 14,  the  magnitude  of  the  first 
peak  of  the  coated  sphere  has  decreased  to  the  point  wliere  it  is  comparable  to 
that  of  the  perfect  conductor.  At  kb  values  beyond  the  first  peak,  the  coated 
sphere  curve  still  has  a larger  amplitude  than  the  perfect  conductor  curve. 
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HflVEVECTOfl  TiMCS  COflTEO  SPHEflE  RflOIUS,  RB.  ‘ 

COATED  METAL  SPHERE,  M2  = 4.0 

RB/RPI  = 1.0400 

I'ig.  3.6.  Scattering  from  a coated  metallic  sphere  normalized  to  the  total  cross 
sectional  area  (solid  curve)  and  scattering  from  a perfectly  conducting  sphere  of 
radius  h (symbols). 
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Fig.  3.7.  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  a (symbols). 
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WRVEVEcTflfl  Tines  coRTen  sPHene  rroius,  rb. 
COfiTED  METAL  SPHERE,  M2  = 3.25 


RB/Rfl  = 1 . 0400 

Fig.  3.  8.  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  b (symbols). 
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COATED  METAL  SPHERE.  M2  = 3.25 


RB/Rfl  = 1.0400 


Fig.  3.  9.  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  a (symbols). 
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HflVEVECTOR  TIMCS  COHTED  SPHERE  flflOIUS.  RB. 

COfiTEO  METAL  SPHERE,  MH  = 2.5 

RB/flfl  = 1.040D 


T.  10.  Scattering  from  a coated  metallic  sphere  normalized  to  the  total  cross 
sectional  area  (solid  ciirv(')  and  scattering  from  a perfectly  conducting  sphere 
of  radius  b (symbols). 
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Fig.  3.  11.  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conducting 
sphere  of  radius  a (symbols). 
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WflVEVECTOfi  TIMES  COATED  SPHERE  RADIUS,  RB. 

COFITED  METAL  SPHERE,  M2  = 1.75 

RB/Rfi  = 1.0400 


Fig.  3. 12.  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  b (symbols). 
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R8/Rfl  = 1.0400 

I'ip.  14.  Scattering  fi'om  n coated  metallic  spliere  normalized  to  the  total 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  b (symbols). 
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RB/Rfl  = 1.0400 
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Fig.  3. 15.  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  a (symbols). 
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Fig.  3.  16.  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  b (symbols). 
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WflVEVECTBfl  TIMES  CORE  SPHERE  RADIUS.  Rfl. 

COATED  metal  SPHERE.  M2  = 1.35 

RB/RA  = 1.0400 


Pip.  3.  17.  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere'  of  radius  a (symbols). 
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U.OU  L.CU  \.eu  C*UU  C.MU  C.9U 

HflVEVECTOfI  TIMES  CMflTED  SPHERE  RADIUS,  RB. 

COATED  METAL  SPHERE.  M2  = 1.20 


R8/RA  = 1.0400 

Fig.  3.  IS.  Scattering  from  a coated  metallic  sphere  normalized  to  the 
total  cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly 
conducting  sphere  of  radius  b (symbols). 
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COATED  METAL  SPHERE.  M2  = 1.20 

RB/RA  = 1.0400 


Fig.  3.  19.  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  a (symbols). 


A slight  shift  to  larger  kb  values  is  noticeable  in  the  coated  sphere  curve  peaks 
with  respect  to  the  peaks  of  tlie  perfect  conductor. 

At  the  still  smaller  coatmg  refractive  index  of  1.2  shown  in  Fig.  3. 18,  the 
shift  in  the  kb  values  of  tlie  peaks  is  much  more  noticeable.  'Ilie  magnitude  of 
the  first  peak  has  decreased  further  until  it  has  a smaller  magnitude  than  that 
of  the  perfect  conductor.  The  second  peak  now  has  a magnitude  comparable  to 
that  of  the  perfect  conductor. 

These  scattering  results  for  a refractive  index  of  1.2  are  plotted  in  terms 
of  the  core  radius  a in  Fig.  3. 19,  instead  of  in  terms  of  the  total  radius  b as 
was  done  in  previously  considered  curves.  The  positions  of  the  peaks  in  tlie 
coated  sphere  scattering  curve  m this  figure  agree  well  with  those  of  a per- 
fectly conducting  sphere  with  radius  a . When  die  amplitude  is  normalized  to 
die  core  radius  a as  it  is  in  this  figure,  die  coated  sphere  curve  has  a larger 
amplitude  than  the  perfect  conductor  cui-ve. 

Similarly,  curves  for  coatings  having  refractive  inc  ces  1.5,  or  greater, 
are  plotted  in  terms  of  die  core  radius  in  Figs.  3.7,  3.9,  3. 11,  3. 13,  and  3. 15. 
The  positions  of  the  peaks  in  the  coated  sphere  curves  and  in  the  curves  for  a 
perfectly  conducting  sphere  of  radius  a , shown  in  these  figures,  are  in  poorer 
agreement  than  they  were  for  a perfect  conductor  of  radius  b shown  in  the  previ- 
ous figures.  In  addition,  the  magnitude  of  die  peaks  in  die  coated  sphere  curves 
are  still  greater  compared  to  the  perfect  conductor  peaks  than  they  were  in  the 
previous  figures. 

Hence,  scattering  results  for  coated  spheres,  analogous  to  those  of 
Scharfman,  are  obtained  for  the  larger  range  of  kb  values  including  the 
first  three  peaks  in  the  scattering  curves:  For  large  values  of  the  coating  re- 
fractive index,  the  oscillations  in  the  scattering  curve  for  die  coated  sphere 
arc  approximately  in-phase  with  those  in  die  scattering  curve  of  a perfectly 
conducting  sphere  having  a radius  equal  to  the  total  radius  of  the  coated  sphere. 
For  coating  refractive  indices  less  than  approximately  1.6,  the  scattering  is 
better  described  by  a perfectly  conducting  sphere  having  a radius  equal  to  the 
core  radius.  In  both  cases,  die  amplitudes  of  the  oscillations  in  die  scattering 
curves  of  the  coated  spheres  arc  greater  than  those  of  the  perfectly  conducting 
spheres. 
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This  qualitative  rule  remains  valid  witli  tliiimer  coatings  where  similar,  but 
much  less  pronounced  effects  are  seen  in  tlie  scattering  plots.  To  illustrate  the 
magnitude  of  the  effect  with  thin  coatings,  scattering  plots  are  presented  for  the 
L radius  ratios  1.004  and  1.0004  and  for  only  two  coating  refractive  indices  of 

1.75  and  4.0  at  each  coating  thickness.  These  cuiwes  are  given  in  Figs.  3.20  to 
3.27.  Since  it  is  difficult  to  see  differences  in  the  curves  at  these  thicknesses, 
the  similarity  to  tlic  previous  case  of  b/a  = 1.04  was  verified  from  the  com- 
puter printouts. 

Including  even  a large  amount  of  nondispersive  absorption  in  the  coating  has 
very  little  effect  on  the  qualitative  behavior  of  tlie  scattering  curves.  Absorption 
lines  with  dispersion  are  discussed  in  Section  4.  To  investigate  the  effect  of  con- 
stant absorption  in  the  coating,  the  coating  refractive  index 

1112  " *”2  " ‘ ^ 

was  used  with  b/a  = 1.0012  and  witli  the  set  of  values  J m2  j = 1.2,  1.35,  1.5, 

1.75,  2.5,  3.25,  and  4.0.  The  scattering  curves  are  qualitatively  the  same  as  is 
expected  for  nonabsorbing  coatings  having  the  same  magnitude  refractive  index. 
Notice  that  the  imaginary  part  of  the  refractive  index  given  in  Eq.  (3.38)  implies 
a large  absorption  coefficient.  With  the  constant  imaginary  part  of  0. 1,  the  ab- 
sorption coefficient  varies  from  220  cnT^  at  kb  = 1 to  3000  cm'^  at  kb  = 4 . 

With  coating  thicknesses  and  refractive  indices  larger  than  those  used  in  the 
above  two-dimensional  grid,  the  scattering  begins  to  depart  from  the  simple  anal- 
ogies with  perfectly  conducting  spliercs.  llic  first  evidence  of  anamalous  results 
for  high- index,  thick  coatings  can  be  seen  in  Fig.  3.6  for  b/a  = 1.04  with  m^  = 4.0. 
In  this  figure  the  scattering  at  tlie  tliird  peak  is  as  large  as  tliat  at  the  second  peak 
in  contrast  to  tlie  decreasing  scattering  amplitude  witli  increasing  kb  of  the  per- 
fectly conducting  sphere.  If  tlie  radius  ratio  is  increased  still  further  to  b/a  = 1.4 
with  m2  = 1.5,  die  magnitude  of  the  third  peak  becomes  greater  than  both  the  first 
and  second  peaks.  This  case  is  plotted  in  Figs.  3.28  and  3.29  where  the  fourth 
peak  is  seen  to  have  a still  greater  magnitude.  Nevertheless,  even  with  diis  very 
thick  coating,  there  still  is  some  qualitative  agreement  widi  a perfectly  conduct- 
ing sphere  of  radius  b in  die  kb  region  of  the  first  two  peaks  in  Fig.  3.28. 
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Fig.  .3.20.  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  b (.symbols). 


CaflTED  METAL  SPHERE,  M2  = 4.0 

RB/RR  = 1.0040 

Fig.  3.21.  Scattering  from  a coated  metallic  sphere  normalized  to  tlie 
core  cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly 
conducting  sphere  of  radius  a (symbols). 
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3.23.  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  a (symbols). 
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KHVEVECTOfi  TIMES  CaflTEO  SPHERE  flflOJUS,  RB. 
COflTED  METAL  SPHERE,  M2  = 4.0 

FIB/RR  = 1,000*4 


Fig.  3.24.  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  b (symbols). 
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Fig.  3.25.  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  a (symbols). 
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CdflTEO  METAL  SPHERE.  M2  = 1.75 

RB/flA  = 1.0004 

Fig,  3.26.  Scattering  from  a coated  metallic  sphere  normalized  to  the  total 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  b (symbols). 
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COATlO  metal  sphere.  M2  = 1,75 

RB/RR  = 1.0004 


pMg.  3.27.  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  a (symbols). 
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HflVEVECTOR  TIMES  CdflTEO  SEMERE  RADIUS,  RB. 
COATED  METHL  SPHERE,  M2  = 1.50 


RB/Rfl  = 1.4000 

l-’ig.  3.28.  Scnm-ring  froni  n coiHcd  mi  tnllic  splierc  normalized  to  the  total 
cross-sectional  area  (solid  curve)  and  scattermg  from  a perfectly  conduct- 
ing sphere  of  radius  b (symljols). 


COATED  METAL  SPHERE.  M2  = 1.50 

RB/RA  = l.UOOO 


Fig.  3.29.  Scattering  from  a coated  metallic  sphere  normalized  to  the  core 
cross-sectional  area  (solid  curve)  and  scattering  from  a perfectly  conduct- 
ing sphere  of  radius  a (symbols). 
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At  kb  values  greater  tlian  tJiose  presented  in  Fig.  3.2«,  the  scattermg 

curve  has  a series  of  sharp  spikes  which  become  increasingly  narrow  and  in- 

28 

crease  in  magjiitude  as  kb  increases,  lliieinstein  lias  investigated  this 
phenomenon  for  perfectly  conducting  core  spheres  coated  with  a nondispersive 
dielectric.  In  his  analysis  tlie  anomalously  large  value  of  the  scattering  at  the 
resonance  spikes  is  attributed  to  constructive  interference  in  "creeping  waves” 
trapped  in  the  coating.  For  small  kb  values  witli  tiiin  coatings  of  low  refractive 
index,  the  "creeping  waves"  arc  considerably  attenuated  as  they  propagate  around 


the  circumference  of  the  sphere  and,  consequently,  have  little  effect  on  the  scat- 
tering. The  creeping  waves  begin  to  influence  the  scattering  when  the  waveguide 
conditions 
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TE  mode 
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(3.39) 

TM  mode 
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arc  met,  where  n is  an  integer,  if  die  coating  diickness  and  refractive  index  are 
held  constant,  once  kli  has  increased  above  tiic  value  required  by  (3.39)  to  turn 
on  a mode,  the  creeping  wave  imalysis  predicts  dial  the  Q of  the  resonances  in- 
creases with  increasing  kb.  Tliis  results  in  increasingly  sharp  spikes  in  the 

scattering  with  increasing  kb . Parameters  for  this  behavior  are  outside  the 

28 

present  range  of  interest  and  reference  is  made  to  Rheinstein  for  details. 

3.  liuantitative  results  for  nondispersive  coatings. 


Consider  scattering  plots  similar  to  those  presented  in  terms  of  die  core 
radius  in  the  previous  section,  where  now  scattering  from  die  uncoated,  metal- 
lic sphere,  rather  than  a perfectly  conducting  sphere,  is  given  by  symbols. 
Since  the  metallic  sphere  is  well  approximated  by  a perfectly  conducting  sphere 
of  the  same  radius,  die  curves  are  not  repeated.  Tlie  qualitative  behavior  re- 
mains the  same  as  iii  the  previous  section  widi  die  scattering  curve  for  the 
coated  sphere  having  large-r  maxima  and  lower  minima  than  the  curve  for  the 
uncoated,  metallic  sphere.  It  is  desired  to  descrilic  die  differences  bctivecn 
the  curves  quantitatively.  To  do  this,  one  complete  oscillation  m the  curve  is 
considered,  say  between  die  first  two  peaks.  This  region  from  ka  = 1.021  to 
2.337  is  divided  into  N points  and  the  root- mean -square  percent  deviation 
(RMS%)  is  calculated; 
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where  G^  and  g’^  arc  the  values  of  the  scattering  for  tiie  coated  and  uncoated  me- 
n n " 

tallic  spheres  at  the  ntli  point,  respectively.  Use  of  the  RMS  % deviation  avoids  dif- 
ficulties with  the  signs  of  tlie  deviations  when  the  curves  cross. 

Using  the  two-dimensional  grid  of  coating  refractive  indices  and  thicknesses 
given  in  the  previous  section,  curves  lor  the  RMS%  deviation  arc  given  in  Fig.  3.30 
for  each  thiclness  as  a function  of  the  coating  refractive  index.  In  the  limit  m2 
= 1.0,  all  of  the  curves  must  go  to  zero  deviation  since  tliis  limit  corresponds  to 
a vacuum  coating.  For  m.,  greater  thiui  appi'oximately  ttvo,  the  curves  become 
nearly  flat.  However,  tlie  values  in  tlie  flat  region  are  not  tJic  asymptotic  values 
for  an  infinite  refractive  index,  i.e. , the  limit  wliere  tiie  coated  sphere  becomes 
a perfectly  conducting  sphere  of  radius  b . 'Jlie  dashed  lines  at  tlie  right-hand 
side  of  the  figure  give  tlie  asymptotic  values  for  the  various  thicknesses.  The 
dotted  line  at  3.7%  is  tlie  asymptotic  limit  of  a coating  having  an  infinite  refrac- 
tive index  and  zero  thickness,  i.e.,  the  limiting  deviation  between  a metallic 
sphere  and  a perfectly  conducting  sphere  of  the  same  radius.  Notice  that  with 
b/a  = 1.04  and  m2  = 4.0,  the  RMS%  deviation  is  greater  tlian  its  asymptotic 
value  in  contrast  to  the  curves  for  thinner  coatings  which  are  below  their  asymp- 
totic limits,  litis  difference  is  a result  of  the  abnormally  large  cross  sections 
generated  by  the  waveguide  mechanism  in  spheres  witli  thick  coatings.  A similar 
enhancement  may  also  occur  for  j m2  1 » 4 in  the  curves  for  thinner  coatings. 

For  tliiii  coatings  with  radius  ratios  less  dian  aliout  1.004,  the  RMS%  devi- 
ation is  witlun  one  percent  of  being  linear  in  the  coating  thickness  at  a fixed 
coating  refractive  index.  This  linear  relationship  is  lost  at  larger  thicknesses. 

For  example,  linear  extrapolation  to  the  b/a  = 1.04  curve  from  the  b/a  = 1.004 
and  1.0004  curves  predicts  values  that  are  tlircc  percent  too  low  at  m2  = 1.2 
and  nineteen  percent  too  low  at  m2  = 4.0. 

The  ciirvx'  with  b/a  = 1.0012  in  Fig.  3.30  presents  the  RMS%  deviation  of 
nondispersive,  absorbing  coatings  having  complex  refractive  hidices  given  by 
Eq.  (3.  38).  Even  witli  the  large  absorption  coefficients  of  these  coatings,  the 
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RMS%  deviation  agrees  very  well  witli  values  for  nonabsorbing  coatings  having 
the  same  magnitude  refractive  index,  except  for  small  differences  at  small  mag- 
nitude indices.  The  dashed  curve  gives  values  for  a real  refractive  index  obtained 
by  extrapolation  from  the  b/a  = 1.004  and  b/a  = 1.0004  curves. 

In  Fig.  3.  31,  the  data  of  Fig.  3.30  have  been  cast  into  level  contours  of  the 
RMS%  deviations  as  a function  of  the  coating  tliickness,  b-a,  and  the  magnitude 
of  the  refractive  index.  As  the  coating  refractive  index  approaches  unity,  i.e. , 
a vacuum  coating,  the  curves  of  constant  deviation  approach  infinity  indicating 
that  an  arbitrarily  thick  vacuum  coating  produces  zero  deviation.  The  determina- 
tion of  asymptotic  values  of  tlie  error  contours  in  the  limit  of  an  infinite  refrac- 
tive index  would  require  another  set  of  computer  calculations  and  hence  has  not 
been  evaluated.  This  limit  corresponds  to  differences  in  tlie  scattering  between 
metallic  spheres  of  radius  a iind  perfectly  conducthig  spheres  of  radius  b. 

4.  Dispersive  coatings. 

Ojitical  constants  of  actual  solid  materials  have  absorption  bands  and  dis- 
persion, rather  than  die  constant  absorbing  and  nonalisorbing  refractive  indices 
considered  in  the  previous  sections.  Those  absorption  bands  can  be  conveniently 
represented  by  a distribution  of  Lorentzian  absorption  lines,  the  use  of  which 
guarantees  the  correct  dispersive  behavior  as  required  by  the  Kramers-Kronig 
relations.  The  Reststrahl  absorption  bands  of  many  ionic  solids  can  be  crudely 
approximated  by  a single  absorption  line  of  large  linewidth.  A better  model  in 
this  case  is  to  allow  die  relaxatior  frequency  of  the  Lorentzian  line  te  be  a func- 
tion of  the  radiation  frequency.  With  electronic  transitions  and  with  lattice  ab- 
sorption in  nonpolar  solids,  the  absorption  bands  consist  of  a superposition  of 
many  absorption  lines  of  narrow  linewidth. 

Hence  as  a model  calculation  applicable  to  many  materials  by  superposi- 
tion, the  effect  of  a single  Lorentzian  absorption  line  on  die  scattering  is  con 
sidered.  'Fhe  dielectric  constant  for  a Lorentzian  line  as  a function  of  frequt  v 
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where  Cq  and  are  the  static  and  high  frequency  dielectric  constants, 
respectively,  is  the  resonant  absorption  frequency,  and  F is  the  relaxation 
frequency.  For  a constant  relaxation  frequency,  F is  equal  to  the  full  width  at 
half  maximum  of  the  absorption  line.  The  scattering  from  an  actual  metal- 
oxide  coating  and  from  a dispersive  dielectric  sphere  is  considered  later  in 
this  section. 

The  effect  of  the  Imewidth  on  the  scattering  was  first  investigated  using 

three  moderately  thick  coatings  all  with  radius  ratio  1.04  and  with  a resonant 

absorption  frequency  given  by  the  ka  value  of  1.428.  The  constants  in  Eq.  (3.41) 

were  adjusted  for  the  three  linewidths  of  50  cm'\  20  cm'^,  and  10  cm’^  used,  so 

that  the  refractive  index  in  the  visible  was  1.5  and  the  absorption  strength  at  the 
4 - 1 

line  center  was  10  cm  . The  real  parts  (solid  curve)  and  imaginary  parts 
(dashed  curve)  of  these  dielectric  constants  are  plotted  in  Figs.  3. 32(a),  3. 33(a), 
and  3.34(a).  The  abscissas  in  these  figures  are  given  in  units  of  kb  for  the  radius 
ratio  b/a  = 1.04.  In  these  units,  the  line  centers  are  located  at  the  kb  value  of 
1.485.  In  parts  (b)  of  Figs.  3.32-3.34,  the  scattering  in  the  vicinity  of  the  ab- 
sorption line  is  plotted  as  a solid  curve  for  the  coated  metallic  sphere.  For 
reference,  scattering  from  a perfectly  conducting  sphere  of  radius  b is  shown 
as  a dashed  curve  on  the  same  axes. 

With  the  50  cm  ^ linewidth  absorption  shown  in  Fig.  3. 32,  the  scattering 
curve  for  the  coated  sphere  has  a gentle  S-shaped  deviation  from  a smooth 
curve  in  the  vicinity  of  the  absorption.  With  the  narrower  20  cm"^  linewidth 
shown  in  Fig.  3. 33,  the  scattering  curve  has  a narrow,  flat  plateau  at  the  ab- 
sorption line  frequency.  With  the  still  narrower  linewidth  of  10  cm"^  shown  in 
Fig.  3.34,  there  is  a small  magnitude  dip  in  tlie  absorption  curve  having  ap- 
proximately the  same  linewidth  as  the  absorption  line. 

To  investigate  the  effect  of  the  absorption -line  frequency  on  the  scattering, 

five  frequencies  are  considered  at  positions  evenly  spaced  in  wavelength  between 

the  first  maximum  and  the  first  minimum  on  the  scattering  curve.  As  above,  the 

scattering  is  calculated  for  the  tliree  linewidths  ( 10  cm'^,  20  cm'^,  and  50  cm"^ ) 

at  each  frequency  keeping  the  radius  ratio  constant  at  1.04  and  the  absorption  co- 

4 -1 

efficient  at  the  line  center  constant  at  10  cm  . These  scattering  results  are 
plotted  versus  ka  in  Figs.  3.35-3.49. 
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Fig.  3.32.  Real  part  and  imaginary  part  of  the  dielectric  constant  of  a coat- 
ing having  a Ix)rentzian  absorption  line  with  a full  width  at  half  maximum  of  50cm"^ 
(a).  Scattering  from  a coated  metallic  sphere  in  the  vicinity  of  the  absorption  ILne(b). 
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fwhm  = 20  cm 


Fig,  3,33,  Real  part  and  imaginary  part  Cj  of  the  dielectric  constant  of  a 
coating  having  a Lorentzian  absorption  line  with  a full  width  at  half  maximum  of 
20  cm  (a).  Scattering  from  a coated  metallic  sphere  in  the  vicinity  of  the  ab- 
sorption line  (b). 
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I'ig.  3.34.  Real  part  anfl  imaginary  part  cj  of  the  dielectric  constant  of  a 
coating  having  a Lorentzian  absorption  line  witli  a full  width  at  half  maximum 
of  10  cm  ^ (a).  Scattering  from  a coated  metallic  sphere  in  the  vicinity  of  the 
absorption  line  (b). 


Fig.  3.35.  Scattering  from  a coated  metallic  sphere  havuig  a 10  cm  ^ fwhm 
absorption  line  at  ka  = 1. 01  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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WRVEVECTOR  TIMES  CORE  SPHERE  RflOIUS.  Rfl. 
COATED  METAL  SPHERE.  LORENTZIAN  LINE. 

RB/RA  = 1.0400 


Fig.  3.36.  Scattering  from  a coated  metallic  sphere  having  a 10  cm  fwhrr 
absorption  line  at  ka  = 1. 12  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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Fig.  3.37.  Scattering  from  a coated  metallic  sphere  having  a 10  cm'^  fwhm 
absorption  line  at  ka  = 1.26  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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Fig.  3.38.  Scattering  from  a coated  metallic  sphere  having  a 10  cm  ^ fwhm 
absorption  line  at  ka  = 1.43  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 


WflVEVECTOfl  TIMES  CORE  SPHERE  RADIUS.  RR. 
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F’ig.  3.39.  Scattering  from  a coated  metallic  sphere  having  a 10  cm  fwhn 
absorption  line  at  ka  = 1.65  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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Fig.  3.40.  Scattering  from  a coated  metallic  sphere  having  a 20  cm'^  fwhm 
absorption  line  at  ka  = 1.01  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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WflVEVECTOR  TIMES  CORE  SPHERE  RflPIUS,  Rfl. 
COATED  METAL  SPHERE,  L0RENT2IAN  LINE. 

RB/RA  = 1.0400 


Fig.  3.41.  Scattering  from  a coated  metallic  sphere  having  a 20  cm  fwhm 
absorption  line  at  ka  = 1. 12  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 


Fig.  3.43.  Scattering  from  a coated  metallic  sphere  having  a 20  cm  ^ fwhm 
absorption  line  at  ka  = 1.43  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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Fig.  3.44.  Scattering  from  a coated  metallic  sphere  having  a 20  cm  ^ fwhm 
absorption  line  at  ka  = 1.65  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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Fig.  3.45.  Scattering  from  a coated  metallic  sphere  having  a 50  cm  ^ f\vhm 
absorption  line  at  ka  = 1.01  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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Fig.  3.46.  Scattering  from  a coated  metallic  sphere  having  a 50  cm  fwhm 
absorption  line  at  ka  = 1. 12  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 


Fig.  3.47.  Scattering  from  a coated  metallic  sphere  having  a 50  cm  ^ fwhm 
absorption  line  at  ka  = 1.26  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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Fig.  3.48.  Scattering  from  a coated  metallic  sphere  having  a 50  cm  fwhm 
absorption  line  at  ka  = 1.43  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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Fig.  3.49.  Scattering  from  a coated  metallic  sphere  having  a 50  cm“^  fwhm 
absorption  line  at  ka  = 1.65  (solid  curve).  Results  are  normalized  to  the 
core  cross-sectional  area  and  compared  to  a perfectly  conducting  sphere  of 
radius  b (symbols). 
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The  effect  of  the  absorption-line  frequency  on  the  scattering  is  similar  for 
all  the  linewidths.  First  consider  die  10  cm  ^ ILnewidth  results  shown  in  Figs. 
3.35-3.39.  The  localized  perturbation  from  the  absorption  line  is  a localized 
valley  near  the  first  peak,  becomes  a shoulder  near  the  midpoint,  and  is  a 
localized  peak  near  the  first  minimum.  Qualitatively  similar  results  are  given 
for  the  20  cm  ^ linewidth  in  Figs.  3.40-3.44  and  for  the  50  cm”^  linewidth  in 
Figs.  3.45-3.49,  The  widths  of  tliesc  localized  deviations  are  approximately 
equal  to  the  absorption  linewidtiis.  Roughly  speaking,  the  magnitudes  of  the 
localized  deviations  are  approximately  proportional  to  the  unperturbed  scatter- 
ing with  the  localized  valley  at  the  first  peak  having  a greater  amplitude  than  the 
localized  peak  near  the  first  minimum.  Both  localized  deviations  are  approxi- 
mately ten  percent  of  tlie  unperturbed  magnitude. 

A careful  comparison  reveals  diat  the  frequencies  of  tlie  centers  of  the 
localized  deviations  and  resonant  frequencies  of  the  Lorentzian  lines  do  not 
coincide.  Instead,  the  localized  deviations  are  displaced  approximately  one 
linewidth  from  the  line  center  toward  larger  ka  values  (higher  frequencies). 

It  is  noted  in  passing  that  a plane  surface  of  bulk  material  having  a single  ab- 
sorption line  is  most  reflecting  at  frequencies  slightly  greater  than  the  absorp- 
tion frequency.  Apparently  because  of  the  Mie -scattering  oscillations  in  the 
scattering  curve,  a small  particle  may  be  eitlier  locally  a better,  or  worse, 
scatterer  at  this  frequency. 

Although  the  most  readily  apparent  effect  of  an  absorption  line  is  the  local- 
ized deviation  centered  near  the  absorption  frequency,  dispersion  in  the  real 
part  of  the  refractive  index  results  in  less  dramatic  changes  over  a larger 
range.  To  make  this  effect  more  apparent,  scattering  curves  for  a 10  cm'^ 
linewidth  absorption  line  in  tlie  very  thick  coating  having  the  radius  ratio  1. 4 
are  shown  in  Figs.  3.50  and  3.51.  When  these  curves  are  carefully  compared 
to  Figs.  3.28  and  3.29  for  the  same  thickness  nondispersive  coating  having  the 
same  refractive  index  in  the  visible,  it  is  seen  diat  the  coated-sphere  curves 
cross  the  perfectly  conducting  sphere  curves  at  slightly  different  positions. 

The  above  results  demonstrate  that  strongly  absorbing  Lorentzian  absorp- 
tion lines  in  sufficiently  thick  coatings  have  an  observable  effect  on  the  scattering. 
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HfiVEYECTOR  TIMES  CORE  SPHERE  RflOIUS,  Rfl. 

COflTEO  METAL  SPHERE.  L0RENTZIRN  LINE. 

RB/Rfl  = 1.4000 


Fig.  3.51,  Scattering  from  a coated  metallic  sphere  having  a 10  cm  fwhm 
.nbsorption  line  (solid  curve).  Results  arc  normalized  to  the  core  cross- 
sectional  area  and  compared  to  a perfectly  conducting  sphere  of  radius  a 
(symbols). 


F"""; ^ 


with  the  deviations  from  a smooth  curve  being  highly  visible  since  they  are 
localized.  With  the  relatively  thick  coatings  having  I'adius  ratios  1.04  and  1.4, 
the  magnitudes  of  the  localized  deviations  are  comparable  to  the  differences 
between  nondispersivc  coatings  and  uncoated  spheres.  Effects  from  thinner 
coatings  and  less  absorbing  lines  are  expected  to  be  less  severe. 

The  real  and  imaginary  parts  of  the  complex  refractive  index  m = n - ik 

of  an  actual  metal  oxide  are  plotted  in  units  of  ka  in  Figs.  3.52  and  3.53,  re- 
29 

spectively.  Two  approaches  can  be  taken  to  put  these  optical  constants  in  a 
form  for  use  by  the  computer  in  calculating  the  scattering  from  a metallic 
particle  coated  with  this  material:  First,  the  imaginary  part  of  the  refractive 
index  can  be  fit  by  a superjiosition  of  several  Lorentzian  lines.  The  correct 
dispersive  behavior  of  tlie  Lorentzian  lines  should  automatically  reproduce  the 
real  part  of  die  refractive  index.  This  set  of  Ixircntzian  lines  can  tlien  be  put  in 
analytic  form  in  a subroutine  of  the  computer  program.  Secondly,  one  can  simply 
tabulate  the  refractive  index  at  the  ka  points  where  the  scattering  is  calculated. 

Since  many  Lorentzian  lines  are  needed  to  fit  the  optical  constants  in  Figs.  3.52 
and  3.53,  the  second  approach  was  taken  in  calculating  tlic  following  results. 

The  scattering  from  a metal-oxide  coated,  metallic  sphere  is  shown  in 
Fig.  3.54  as  a function  of  ka  for  a radius  ratio  of  b/a  = 1.04.  On  the  average, 
the  scattering  is  similar  to  that  from  a sphere  coated  witli  a nondispersivc  coat- 
ing of  the  same  thickness.  However,  there  are  four  narrow,  localized  valleys  in 
the  scattering  curve:  two  of  small  amplitude  at  the  ka  values  of  0.76  and  1.00 
and  two  of  large  amplitude  at  the  ka  values  of  1.39  and  1.43.  Tlie  positions  and 
magnitudes  of  diese  localized  valleys  do  not  correlate  well  with  the  positions 
and  magnitudes  of  the  absorption  peaks.  Time  did  not  permit  a detailed  study 
of  this  effect.  It  appears  that  the  localized  effects  of  dispersion  on  the  scatter- 
ing may  be  most  severe  when  the  real  part  of  the  refractive  index  is  small  and 
the  imaginary  part  is  large.  A possible  explanation  of  the  large  magnitude, 
localized  valleys  is  that  the  small  real  part  of  the  refractive  index  allows  much 
of  the  incident  radiation  to  penetrate  the  coatmg.  This  radiation  is  then  absorbed 
by  the  large  imaginary  part  of  the  refractive  index.  Any  future  studies  of  this  ] 

effect  should  also  consider  tlic  possibilities  of  subtle  resonance  effects  in  the  i 

Mie  scattering  or  an  effective,  localized  horizontal  displacement  of  the  scattering  I 
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WRVEVECTOR  TIMES  CORE  SPHERE  RflOIUS,  Rfl. 

COATED  METAL  SPHERE,  ABSORPTION  BANDS 
IN  COATING.  RB/RA  = 1.0400 


Fig.  3,54.  Scattering  from  a metallic  sphere  coated  with  a metal-oxide 


at  such  points.  A small,  localized  horizontal  displacement  m a rapidly  varying 
portion  of  tlic  scattering  curve  can  result  in  a large  vertical  displacement. 

Recall  that  a similar  effect  was  observed  above  for  a single  Lorentzian  line 
where  it  was  observed  that  die  center  of  the  localized  deviation  in  the  scattering 
curve  was  slightly  displaced  toward  higher  frequencies.  In  diis  region  approxi- 
mately one  linewiddi  above  the  resonant  frequency,  the  real  part  of  die  refrac- 
tive index  passes  through  a minimum  wliile  die  imaginary  part  is  slowly  varying 
and  relatively  small.'^^ 

The  real  and  imaginary  parts  of  the  refractive  index  of  an  absorbing  dielec- 

29 

trie  material  arc  plotted  in  Figs.  3.55  and  3.56.  Particles  of  this  material, 
smaller  in  size  than  the  primary'  particles,  may  be  useful  as  spacer  particles 
to  prevent  contact  and  cohesion  of  the  primary  particles.  Tlie  scattering  from 
the  homogeneous  dielectric  spheres  is  plotted  as  a function  of  ka  in  Fig.  3.57. 
Tlie  ka  range  in  this  figure  is  smaller  dian  previous  figures,  extending  only  to 
the  first  peak,  because  the  wavelength  range  was  kept  constant  and  a smaller 
particle  radius  was  used.  In  diis  region  the  scattering  from  the  dielectric 
spheres  is  less  efficient  dian  from  metallic  spheres  of  die  same  size.  There 
is  one  localized  step  in  the  otiierwise  slowly  varying  scattering  curve  at  the 
ka  value  of  0.32.  At  diis  point,  die  real  part  of  the  refractive  index  is  slowly 
varying  but  the  imaginary  part  drops  abruptly  from  one  slowly  varying  value 
to  another. 


E.  Summary  and  Conclusions  of  Scattering  Study 

Scattering  from  spheres  is  typically  presented  as  a plot  of  the  scattering 
efficiency  as  a function  of  ka,  w^ere  k = 2tt/X  is  the  propagation  constant  of 
the  radiation,  X is  the  radiation  wavelength,  and  a is  die  sphere  radius.  In  pre- 
vious sections,  scattering  from  coated  spheres  has  been  presented  both  in  terms 
of  the  core  radius  a and  total  radius  of  the  core  plus  the  coating  diickness  b . 
The  plot  of  die  scattering  for  a metallic  sphere  is  an  oscillating  curve  having 
peaks  near  the  ka  values  of  1.0,  2.3,  3.6,  etc.,  with  the  peaks  decreasing  in 
amplitude  as  ka  increases.  There  is  no  simple  niediod  for  calculating  the  scat- 
tering when  ka  is  of  the  order  of  unity,  hi  diis  region  it  is  necessary  to  use  tile 
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Fig.  3.57.  Scattering  from  a small,  homogeneous  sphere  of  absorbing 
dielectric  material. 
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complicated  Mie-theory-serics  solution.  This  is  in  contrast  to  the  simple 
approximations  of  geometrical  qitics  ;md  Rayleigh  scattering  (first  term  in 
Mie  series)  that  can  1k'  made  m tiie  limits  of  very  large  and  very  small  ka 
values,  respectively. 

Mie-scattering  results  for  perfectly  conducting  spheres  can  be  used  as  lim- 
iting approximations  for  scattering  from  coated  and  uncoated  metallic  spheres. 
However,  for  the  determination  of  the  effects  of  coatings  on  the  scattering  of  a 
metallic  sphere,  this  limit  cannot  be  used  because  the  effects  of  the  coatings 
are  often  smaller  than  the  differences  in  die  scattering  between  actual  metallic 
spheres  and  perfectly  conducting  spheres.  The  Mie -theory  calculations  per- 
formed in  this  study  may  be  the  first  such  calculations  for  an  actual  metal. 

These  results  for  coated  and  uncoated  metallic  spheres  have  been  compared 
qualitatively  and  quantitatively  to  scattering  from  perfectly  conducting  spheres. 

With  an  uncoated  metallic  sphere,  the  positions  of  the  maxima  and  minima 
in  the  oscillations  of  the  scattering  curve  agree  well  with  those  of  a perfectly 
conducting  sphere.  The  amplitude  of  the  oscillations  is  greater  for  an  actual 
metal  than  for  a perfect  conductor,  with  the  differences  becoming  greater  with 
increasing  ka.  The  maximum  difference  occurs  at  the  first  two  minima  where 
the  scattering  curve  of  the  metallic  sphere  is  approximately  ten  percent  below 
that  of  a perfectly  conducting  sphere.  The  mean  difference  in  the  region  between 
the  first  two  peaks  is  four  percent. 

For  coated  spheres  having  coating  refractive  indices  greater  than  approxi- 
mately 1.6,  the  scattering  is  qualitatively  similar  to  that  from  a perfectly  con- 
ducting sphere  having  a radius  equal  to  the  total  of  the  metallic  core  radius  plus 
the  coating  thickness.  With  coating  refractive  indices  less  than  this  value,  scat- 
tering from  the  coated  sphere  is  better  approximated  by  that  from  a perfectly 
reflecting  sphere  having  the  core  radius.  These  respective  approximations  are 
obvious  in  the  limits  of  a perfectly  reflecting  coating  (infinite  refractive  index) 
and  of  a vacuum  coating  (refractive  index  of  unity).  The  small  value  of  the  re- 
fractive index  at  the  cross-over  point  is  surprising  in  view  of  the  small,  bulk 
reflectivity  of  such  a coating.  In  both  comparisons  the  scattering  curves  for 
the  coated,  metallic  sphere  and  the  appropriate  radius  perfectly  conducting 
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sphere  are  approximately  in-phase,  having  maxima  and  minima  at  nearly  the 
same  positions.  'Flic  amplitudes  of  the  coated,  metallic  sphere  curves  arc 
greater  than  those  of  the  corresponding  perfectly  reflecting  sphere  curves  with 
the  difference  increasing  widi  increasing  ka . 

With  coating  thicknesses  and  refractive  indices  greater  than  those  of  cur- 
rent interest,  waveguide  modes  can  be  set  up  in  the  dielectric  coatings.  In 
coatings  where  this  occurs,  these  modes  result  in  large  magnitude,  narrow 
spikes  in  the  scattering  curve  at  large  ka  values. 

To  investigate  the  effects  of  the  coatings  quantitatively,  the  root-mean- 
square  percent  deviation  (RMS%)  between  the  coated  sphere  scattering  and  that 
from  the  uncoated  metallic  sphere  core  was  calculated  for  the  ka  region  between 
the  first  two  peaks  in  the  scattering  curve.  With  thin  coatings  that  have  a ratio 
of  the  coated  sphere  radius  to  the  core  radius  less  than  approximately  1.04,  the 
RMS%  deviation  is  approximately  linear  in  the  coating  thickness  at  a fixed  re- 
fractive index.  For  coatings  having  refractive  indices  less  than  4.0,  the  RMS% 
deviation  is  less  than  0.3%  for  a radius  ratio  of  1.0004,  less  than  3%  for  the 
radius  ratio  1.004,  and  less  than  30%  for  the  radius  ratio  1.04. 

With  nondispersive  absorbing  coatings,  the  scattering  results  are  quali- 
tatively and  quantitatively  the  same  as  those  for  nondispersive  transparent 
coatings  having  tlie  same  magnitude  refractive  indices,  except  for  small  quan- 
titative differences  at  small  magnitude  indices. 

Actual  absorbing  coating  materials  have  dispersive  absorption  bands.  The 
effect  of  dispersive  absorption  on  the  scattering  was  first  studied  using  a model 
coating  having  a single  Lorentzian-shaped  absorption  line  of  varying  linewidths 
and  varying  frequencies.  The  use  of  a Lorentzian  line  guarantees  the  correct 
dispersive  behavior  of  the  refractive  index.  In  this  study  the  coating  thickness 

was  held  constant  at  a radius  ratio  of  1.04  and  the  absorption  coefficient  at  the 

4 -1  -1 

line  center  was  held  constant  at  10  cm  . Three  absorption  linewidths  (10  cm  , 

20  cm  \ and  50  cm  were  examined  at  five  frequencies  equally  spaced  in  wave- 
length between  the  first  maximum  and  the  first  minimum  on  tlie  scattering  curve, 
'fhe  absorption  lines  resulted  in  localized  deviations  in  the  scattering  curve  lo- 
cated near  the  resonant  frequency  of  tlie  Lorentzian  line  and  having  approximately 
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the  same  width  as  the  absorption  line.  A close  examination  reveals  that  the 
localized  deviation  is  displaced  in  frequency  approximately  one  lincwidtli  toward 
the  high  frequency  side  of  tlie  absorption  line  at  the  point  where  tlic  real  part  of 
the  refractive  index  is  a minimum. 

The  shape  of  the  localized  deviation  is  primarily  dependent  on  the  frequency 
of  the  absorption  line  and  is  dependent  to  a lesser  extent  on  the  Imewidth.  The 
results  for  the  10  cm  ^ linewidth  are  representative:  At  the  first  maximum  in 
the  scattering  curve  the  absorption  line  results  in  a localized  valley  of  approxi- 
mately 10  cm  ^ linewidth.  Near  the  midpoint  between  the  first  maximum  and 
the  first  minimum,  the  localized  deviation  becomes  a shoulder  in  the  scattering 
curve.  A localized  peak  is  observed  for  frequencies  near  the  first  minimum. 
The  magnitudes  of  the  localized  deviations  are  approximately  proportional  to 

the  local  value  of  the  scattering  curve.  With  the  radius  ratio  of  1.04  and  the  ab- 

4 -1 

sorption  coefficient  of  10  cm  used,  both  the  local  valley  at  the  first  maximum 
and  the  local  peak  at  the  first  minimum  have  a magnitude  of  approximately  ten 
percent  of  the  value  of  the  scattering  curve  at  these  positions  for  aU  absorption 
linewidths. 

The  localized  deviations  from  absorption  bands  in  an  actual  metal-oxide 
coating  were  also  investigated.  It  is  observed  that  magnitudes  of  the  localized 
deviations  in  the  scattering  curve  do  not  correlate  well  with  the  peaks  in  the 
absorption  coefficient  of  tite  coating.  Rather,  the  most  severe  localized  effects 
appear  to  occur  when  the  real  part  of  the  refractive  index  is  small  and  the  im- 
aginary part  is  large. 

Scattering  from  homogeneous  dielectric  particles  smaller  in  size  than  the 
primary  particles  was  also  examined.  On  the  average,  the  scattering  is  less 
efficient  than  from  the  same  size  metallic  particle.  There  is  one  localized 
valley  in  the  scattering  curve  at  a point  where  the  real  part  of  the  refractive 
index  is  small  and  slowly  varying  and  the  imaginary  part  changes  abruptly  from 
one  slowly  varying  value  to  another. 

With  the  1.04  radius  ratio  and  the  absorjition  coefficient  of  10*^  cm  ^ used 
in  the  study  of  the  Lxjrentzian  absorption  lines,  the  magnitudes  of  the  localized 
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deviations  are  comparable  to  those  from  nondispersivc  coatings  of  the  same 
thickness.  However,  in  contrast  to  nondispersivc  coatings  where  the  deviations 
are  spread  over  large  ka  ranges,  the  deviations  caused  by  absorption  lines  are 
more  visible  because  tliey  are  localized.  These  effects  are  expected  to  be  less 
severe  with  thinner  coatings  and  smaller  absorption  coefficients  and  more  severe 
with  thicker  coatings  and  larger  absorption  coefficients.  The  localized  effects  of 
the  absorption  lines  may  bo  ameliorated  somewhat  in  practical  situations  where 
there  is  a distribution  of  sphere  sizes.  In  this  case  some  smoothing  is  expected 
from  the  fact  that  the  localized  deviation  from  a fixed  wavelength  absorption  line 
may  be  a local  valley  for  one  size  and  a local  peak  for  another  size  yielding  par- 
tial cancellation. 
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APPENDIX  A.  DESCRIPTION  OF  THE  SUBROUTINES  AND  THEIR  FUNCTION 


( 


Listings  of  the  computer  codes  may  be  found  in  Appendix  D. 

SCATER:  Main  program.  The  main  program  has  four  nested  loops.  The 
outermost  loop  is  executed  NCASES  times.  It  is  normally  used  to  run  sequences 
of  spheres  of  different  radii,  but  may  also  be  used  to  control  the  plotting  of  each 
case  or  to  vary  the  range  over  which  the  RMS  % deviation  is  computed.  The  next 
loop  is  executed  MVM  times  and  is  used  to  initialize  the  refractive  index  formu- 
las. The  third  loop  varies  ka  from  zero  to  4.0.  The  refractive  indices  are  eval- 
uated at  each  wavelength  by  calls  to  COREMl  and  COATl,  and  then  all  the  various 
arguments  of  the  Riccati -Bessel  functions  are  formed.  The  next  section  evaluates 
the  psi  and  chi  Riccati-Bcssel  functions  by  calls  to  INITAL  which  in  turn  calls 
PSI,  PS1910,  and  FOR  WAD. 

The  innermost  loop  completes  the  computation  of  the  Riccati -Bessel  functions 
by  calls  to  RECURR  and  sums  the  scattering  coefficients  of  orders  1-10  to  get  the 
scattering  amplitudes.  Upon  leaving  the  innermost  loop  the  scattering  is  computed 
and  the  wavelength  range  has  been  covered,  the  RMS  % deviation  is  computed  by 
calls  to  SUMSQ,  and  the  results  are  plotted  by  calls  to  LFARGE  and  PLOTTT. 

PSI910:  generates  4)g(z)  and  from  the  sine-cosine  polynomials  (3.28) 

and  (3.29)  for  large  arguments. 

PSI:  generates  ;^)g(z)  and  !/)jq(z)  from  the  power  series  (3. 30)  for  small 
arguments. 

FORWAD:  is  used  to  generate  by  backwards  recursion  from  and 

or  to  generate  by  forward  recursion  from  Xj  and  X2* 

RECURR:  computes  the  zeta  functions  and  all  function  derivatives  of  any 
given  order  V,  between  one  and  ten,  and  tlien  computes  the  Vth  order  scattering 
coefficients  for  the  core  and  coated  spheres.  Common  area  REX  passes  argu- 
ments and  function  values  between  RECURR  and  SCATER. 

SUMSQ;  calculates  the  RMS%  deviations  between  the  coated-sphere  gain 
curve  and  the  perfectly  conducting  or  core-sphere  gain  curves,  over  a specified 
range  of  ka  and  kb : 
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i = KF 

RMS%  deviation  = 100  I X/ 

\ i = KI 


(g1(0  - G2(i))/G2(i) 


CQATM 2 : entries  COATl  and  COREMl.  GOATM2  is  a user-supplied  routine 
that  provides  the  dispersion  relations  for  the  core  and  coating  refractive  indices. 
For  each  pair  of  radii  (or  each  value  of  MV  in  the  range  1 ^ MV^MVM),  SCATER 
malci's  a call  to  COATM2  before  each  run  through  its  range  of  wavelengths.  This 
call  should  be  used  to  initialize  the  dispersion  relations  of  the  core  and  coating 
for  that  particular  run.  Subsequently  SCATER  maJees  successive  calls  to  COATl 
and  COREMl  for  the  complex,  double  precision  values  of  M2  and  Ml  at  each  wave- 
length. 

Four  routines  already  exist  and  can  serve  as  models.  (See  Appendix  D for 
listings.)  One,  COAT,  provides  for  linear  dispersion  in  die  real  and  imaginary 
parts  of  the  refractive  indices,  while  another,  COAL,  provides  for  a Lorentzian 
absorption  line  in  the  coating  material  widi  linear  dispersion  in  the  core.  Nu- 
merical tables  for  realistic  complex  refractive  indices  are  used  in  COAS  for  a 
metal-oxide  coating,  and  in  COAC  for  a solid  dielectric  sphere. 

LFARGE:  entry  PLOTTT.  Tliere  is  one  call  to  LFARGE  for  each  set  of 
axes.  Calls  to  LFARGE(NCK)  give: 

NCK  = 1;  initializes  plotter  by  a cal  to  PINIT. 

NCK  = 1 or  2;  reads  cards  containing  x and  y axis  lengths,  variable  ranges, 
plot  scale  factors  and  axis  and  graph  titles.  Then  axes  and  titles  are  plotted. 

A call  to  PLOTTT  is  made  for  every  data  curve  to  be  plotted.  Any  numlDcr  of 
curves  may  be  placed  on  one  set  of  axes.  The  last  call  to  PLOTTT  for  a given  set 
of  axes  also  controls  the  advance  to  the  next  plot  if  NTG  = 2,  and  closes  the  plotter 
data  set  on  the  last  call  of  a complete  run  by  setting  NTG  = NENDPL  = 2 . 

Note  that  the  plotter  routine  may  change  the  data  arrays  given  to  it.  Presently, 
if  y-axis  array  elements  have  magnitudes  outside  the  range  of  the  y axis  drawn  by 
the  previous  call  to  LFARGE,  the  offending  elements  are  truncated  to  the  nearest 
limit.  Further,  tlie  plotter  routines  require  storage  of  the  minimum  value  plotted 
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on  the  axis  and  the  number  of  units  per  inch  to  be  stored  near  the  end  of  the 
arrays.  Normally,  if  all  elements  are  plotted,  diose  data  are  stored  in  the  last 
two  array  elements  which  are  reserved  for  that  purpose.  Additional  explanation 
may  be  found  in  the  comments  in  LFARGE  in  tlie  listing  of  ^pendix  D. 

The  routines  called  by  LFARGE  are  standard  CALCOMP  software  and  de- 
scriptions may  be  found  in  the  manufacturer's  ))lotter  manuals.  There  are  three 
exceptions: 

PINIT  (XX):  initializes  the  plotter  output  on  Fortran  device  number  XX, 
and  writes  a user  ID  on  the  plot. 

DFACT:  sets  both  plot  scale  factors.’ 

ENPLT:  closes  plotter  output  file  and  plots  user  ID. 
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APPENDIX  B.  INPUT  DATA  CARD  DESCRIPTION 

Tills  appendix  describes  tlu>  input  data  cards  required  by  the  computer 
program  SCATER  and  its  subroutines.  Included  are  card  ordering,  variables 
read,  and  Fortran  formats,  as  well  as  a description  of  tlie  function  of  the  vari- 
ables in  the  program  and  a sample  data  card  set. 

On  some  of  the  data  cards,  blank  formats,  such  as  "lOX,"  are  specified  to 
allow  the  user  to  identify  the  variable  on  the  data  card  by  its  name,  a mnemonic 
or  its  format. 


CARD 
Cards  1-5: 


VARIABLES 

COMENT 


FORMAT 

(20A4) 


These  cards,  which  can  contaijj. alphanumeric  characters,  provide  a five- 
line  descriptive  title  for  each  run.  They  may  be  left  blank  but  cannot  be  omitted. 

Car^;  NCASES  (lOX,  14) 

Tliis  card  contains  the  nuir±)er  of  different  pairs  of  radii,  (a, b),  to  be  studied 
in  the  run. 


The  user  must  provide  one  set  of  all  the  following  cards  for  each  pair  of 
radii,  NCASES  sets  in  all. 


Card  7 : 


RA,  RB 


2(10X,D24. 17) 


Reads  core  and  coated  sphere  radii. 


Card  8: 


MVM 


NPLOT 


KI,KF 


MVM.  NPLOT,  Kl,  KF  4(10X,  14) 

is  the  number  of  different  refractive  indices  or  materials  to  be 
run  for  this  pair  of  radii. 

controls  die  flow  of  data  to  die  plotter:  plots  if  NPLOT  = f,  no 
plot  if  NPLOT  = I. 

is  a pair  of  integers  used  by  subroutine  SUMSQ.  SUMSQ  finds 
the  RMS%  deviation  of  the  gains  over  a ratige  of  ka  whose  upper 
and  lower  lx)unds  are  APLOT(KI)  and  APLOT(KF),  respectively. 
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Tlie  user  must  supply  one  set  of  all  the  following  cards  for  each  material- 
M\'M  sets  for  each  pair  of  radii. 

Card  9:  RMl,  SRMl,  EMI,  SEMI  (7X.D13.6) 

Card  10:  RM2,  SRM2,  EM2,  SEM2  (7X,  D13. 6) 

These  cards  are  read  by  subroutine  COATM2  to  provide  the  initialization 
constants  for  complex  slope -intercept  dispersion  relations  for  the  core  and 
coating.  If  the  user  supplied  a different  subroutine  for  calculating  the  refrac- 
tive indices,  the  number  and  format  of  these  cards  could  change. 

• •• 

The  following  cards  arc  read  by  subroutine  LFARGE  to  provide  data  to 
the  plotter.  They  are  provided  only  if  NPLOT  = 1.  There  are  two  sets  of 
these  cards  for  each  material  or  refractive  index.  The  first  set  is  for  a plot 
of  and  vs  ka , while  the  second  is  for  a plot  of  and 

vs  kb. 


Card  II:  LABE  LX  (15A4) 

60  character  alphanumeric  title  for  tlie  x axis. 

Card  12:  LABELY  (15A4) 

60  character  alphanumeric  title  for  the  y axis. 

Card  13:  LADELG.  LABELH  2(10A4) 

Two  line,  80  character  alphanumeric  title  for  the  graph. 

Card  14;  XML\,  XMAX,  SIDEX,  XFCT,  YMIN,  (8G10) 

YMAX,  SIDEY,  YFCT 

The  first  four  variables  refer  to  the  x axis  and  the  second  four  to  the 
y axis.  We  describe  the  x-axis  variables  only  since  the  y-axis  variables  have 
the  equivalent  meanings. 


YMIN,  XMAX 


arc  the  minimum  and  maximum  values  of  the  x variable 
written  on  the  x axis. 
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SIDEX 


is  the  length  of  the  x axis  in  inches. 

is  a scaling  factor  that  maltiplies  the  x dimension  of 
graphs,  titles,  and  axes. 


XFCT 


It  is  suggested  that  if  the  over-all  dimensions  of  the  graph  axis  are  signi- 
ficantly less  than  8"  X 8”,  this  should  be  accomplished  by  making  XFCT  and 
YFCT  less  than  unity.  XFCT  and  YFCT  reduces  the  sizes  of  titles,  axes, 
and  graphs  together  thus  keeping  them  proportionate. 

Also  note  that  in  most  cases  and  have  a maximum  value  of 

2 

about  4.0  when  plotted  vs  kb  and  normalized  to  b . But  will  be  larger 

by  the  ratio  (b/a)^  when  plotted  vs  ka  and  normalized  to  a^.  This  should  be 
taken  into  account  when  choosing  the  YMAX  values  if  (b/a)  is  significantly  dif- 
ferent from  one. 

Figure  B.  1 is  a listing  of  a tNqhcal  set  of  data  cards  for  the  coated-sphere 
scattering  program  when  tlie  coating  has  a Lorentzian  absorption  line.  In  the 
two  columns  on  the  right-hand  side  of  the  figure,  the  cards  are  given  a sequen- 
tial card  number  and  a card-type  number  which  refers  to  the  card  descriptions 
given  above.  The  repetition  of  the  card  types  is  evident  from  the  figure.  Briefly, 
since  NCASES  = 2 on  card  # 6,  we  have  two  complete  data  groups  beginning  at 
cards  #7  and  #39.  On  card  #8,  MVM  = 3,  so  three  different  material  subcases 
follow.  In  particular,  these  three  sets  of  data  cards  were  used  to  generate  the 
graphs  in  Figs.  3.35-3.42  for  metal  cores  clad  with  materials  having  different 

Lorentzian  absorption  lines.  The  second  group  beginning  at  card  #39  has  a = 2.5 

2 

and  b = 3.5,  so  (b/a)  ~ 2 and  the  coated-sphere  gain  curve  will  be  about  twice 
as  large  as  the  perfectly  conducting  sphere  curve  when  plotted  vs  ka . TTierefore 
on  card  #46  we  have  douliled  the  length  of  die  y axis  in  the  manner  suggested 
above. 
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Fig.  B.  1.  lyjjical  set  of  data  cards  for  the  coated-sphere  scattering  prograin. 
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APPENDIX  C.  LIST  OF  VARIABLES 

The  variables  in  all  program  segments  except  SCATER  ajid  RECURR  arc- 
well-defined  by  the  subroutine  dc'scriptions  (A^ipc-ndix  A),  llu'  input  data  card 
descriptions  (Appendix  13),  and  comment  cards  in  the  program  listings  (Appen- 
dix D).  This  appendix  contains  a descriptive  list  of  the  important  variables  used 
in  SCATER  and  RECURR  other  than  the  input  variables  described  in  Appendix  B. 

The  names  of  all  the  Riccati -Bessel  functions  and  their  derivatives  are  con- 
structed from  the  following  mnemonic: 

1)  The  first  letter  of  the  name  indicates  the  type  of  Rlccati-Bessel 
function  according  to: 

ip^S.  x^X,  4-^  z. 

2)  The  remainder  of  the  name  is  the  computer  name  of  the  argument: 

A,  N,  MIA,  M2A,  M2N  . 

3)  The  Ip  and  X functions  are  stored  in  arrays  whose  subscript  equals 
the  order  of  the  Bessel  function. 

In  all  other  cases  the  order  is  the  current  value  of  the  integer  vari- 
able V. 

4)  The  name  is  preceded  by  a "D"  if  it  is  tlie  derivative  of  the  function. 

In  Table  C.  1 all  arrays  are  underlined  and  all  arrays  used  by  the  plotter 
end  in  "PET." 
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Variable 

Table  C.l.  Description  of  program  variables. 
Description 

KIP 

= 

number  of  frequency  (or  wavelength)  values 

KEP 

= 

dimension  of  plotter  arrays 

PI2 

= 

2 IT  j 

I 

= 

1 

complex  imaginary  unit 

ZERO 

= 

complex  zero 

RA 

= 

core  sphere  radius,  a 

RB 

= 

coated  sphere  radius,  b 

L 

= 

wavelength,  X 

Ml 

= 

core  refractive  index,  m^^ 

M2 

= 

coating  refractive  index,  m2 

A,  DA 

= 

a = ka  I 

N.DN 

= 

, second  variable  is  complex 
i;  = kb  ) 

APLT 

= 

ka 

BPLT 

= 

kb 

MIA 

= 

mj^a  = mj^ka 

M2A 

= 

m2«  = m2ka 

M2N 

= 

m2V  = m2kb 

TA 

= 

complex  sine  of  Bessel's  function  argument 

TB 

= 

complex  cosine  of  Bessel's  function  argument 

V 

= 

order  of  the  current  scattering  coefficient  ' 

VPl 

= 

V+  1 

VMl 

= 

V - 1 , 

Table  C.  1.  Description  of  program  variables  (Continued). 


Variable 
AN,  BN 

ALN,  BLN 

G 

GI 

GAL 

GPLT 

GIPLT 

GALLON 


Description 

= sequentially,  these  become  the  scattering  coefficients  for 
the  coated  sphere  and  then  the  perfectly  conducting  sphere, 
of  order  V 

= core  scattering  coefficients  of  order  V 
= scattering  amplitude  of  coated  sphere 
= scattering  amplitude  of  perfectly  conducting  sphere 
= scattering  amplitude  of  core  sphere 

2 

gain  of  coated  sphere  normalized  to  a 

2 

gain  of  perfectly  conducting  sphere  normalized  to  b 

2 

gain  of  core  sphere  normalized  to  b 
= gain  of  coated  sphere  normalized  to  b^ 
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APPENDIX  D.  COMPUTER  CODE  LISTINGS 


L- 


CH 

a\ 


o — (\jro<ju^'ON't>'7'o-*rg'^'<fir»'CiNao(7'o  — rjn«rir<CN(riOO'-fvj^<tir'£;K(na'0  — run<jir)'iJNCPC‘o  — 
oooooooooo*^  — — — — 

OnOCOC*OOOOOOOOOOOOOOooOOOOOOC>OooOOOOOOO« 

X’Z  Z ^ Z ZZZ  TZi:  Z 7,  ^ ZZ  7 /,/  7.  Z 7 /.z.  ^ 7.  V ^ ^ 7.  ^7.7  y 7 Z.Z  Z Z'Z  Z.^'/  iL. 

<<<<<<<<<<<<<<<<<<<<<<<<<«f<<<<f<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
Z27ZX222Z2ZS7  72.2.>2S^7S>^3L2:r:^Z^:^2^32.2.2.?:2:222  2722ZZl^.  22.SZZ>Z:^2:2 


H 

<1 

a 


X X 

U. 

• • < 

o 

< 

»-  >»  z 

• 

• -J  < 

z 

• 

N 

iAZ  • 

o 

v> 

• 

U, 

• 

D 

• 

Z fU 

in 

O 

X 

• 

• a:  o 

z 

K 

w 

II 

s 

^Uicg 

lU 

z 

V 

in 

IV 

• * K w 

2 

s 

D 

««  • 

Z O K K 

M 

a 

w 

V fU 

< -Q  J 

o 

« 

a 

A 

O 

V ♦ 

a o ja 

N 

» <1 

< 

» « 

o » a ID 

• 

liJ 

rvj 

2 

• w 

or 

OCD 

OCD  O 

z » 

X 

o 

o 

M 

• 

• 

uJQr 

a *a  • 

-J 

K 

ru 

J 

• 

M • 

cr 

X 

»- 

c *r  o«^ 

C50 

* 

u. 

• 

• 

UJ 

n 

'to 

wu.ru 

» — o 

m 

o 

cr 

m ji  in 

K 

n 

Ok 

o o o 

Z w o -« 

* < 

o 

w 

II 

1 2 

.7 

• 

u* 

Z 

-J  < - w 

fU 

o 

ZK  O 

O 

• 

- o 

a 

< w <r 

a 

s 

z 

w 

•> 

eg 

•UJ 

ax  < In 

• 2 Z X 

UJ 

a 

D 

M 

h-  -l»- 

•> 

Z 

Xfsj 

luu:  J 

Z.  l/iCM  •> 

X 

UJ 

o 

< Cl  < 

X 

• 

in- 

V)  a 

(D  Z -* 

X 

in 

< 

o z 

o 

• 

• J 

I-  a >-  CD 

• i/)o 

o 

• 

in 

a 

o cr 

X 

V < • 

< z < • 

z - 

z 

a 

lU 

< 

(\J 

K 2 or 

uocr  *- 

< w 

< 

M 

a 

M- 

H ••  > 

cu 

K cr  o 

in  u a 

• < 

X 

a 

z 

Z 

• 

» 

— OK 

• < o 

z in 

in 

• 

a 

UJ 

UJ«X 

N 

• 

z z u 

M 

• • 

UJ 

N 

D 

a 

a II  UJ 

• 

IL 

•“« 

UJ 

> >•  O-' 

< 

D • 

O 

in 

Hi 

UJwQ 

< 

V 

z 

L.  • > 

< • Z >- 

* • 

-Jf'J  o 

IL 

u. 

u.  z 

J 

• 

• 

Z X 

» ^ - J 

z 

< ■*. 

o 

u. 

U.  K — 

D 

II 

• 

O UJ 

I « 0. 

fu  o 

>a  o 

K 

M 

* 

-o 

2 

in 

X 

IDDO- 

2 Z 

>-• 

in 

o 

Q JUJ 

cr 

z> 

N 

» » < 

Z 

• w . S4 

> X »£> 

-J 

X- 

a > 

o 

•-< 

X • » 

< a J • 

< Z -^O 

IJ 

CD 

o 

^ u. 

» » 

u.  - 

u. 

o 

• 

onj 

J I O'^ 
a D u.  rsi 


Z f\J  X o 

> 2 


z M m V a 
uj  cr  ^ »- 


H X 
OO 


• C£ 

-J  » * O 


a 

a 

Ol 

(T 


OO  w — 

Z U 05  7 Z O 

»-v-#oo—  • — 

ZO<»~'--OZ  *sz 
n o lu  — H fvj  * > X io 

0 • cr  </>  j — ^ u] 

UJ  • ua  z • ^ a 

if)  •►-Qrfn<#\ 

- -i  ct 

1 *un^  -Jtnu'Z 


o a^'Otn*4-z  zvtu  uj 
— ^(tC-'OUJCDOliJ  U !<»</>  l/)X 
'-’UJ<f'J*«'200>  Z<<-‘<-- 

f-fT  * a <\»o  nou<fCJS 
7 7 a«x50fr'-''-u<— u(vz^zv 
lUDuiv  •r5»~  •••*•• 

7 Z»-fV<OfruIOO<<?  >rXN 
# » O h w cr  V)  < ) o o f.i  I o o o r 
— -‘UuJuaN.UJC*'  o — 

O O X -J  O ft  C w ^ ^ ^ w o 

K n rr  2 iiJ  o ¥-~  K K M 


c 

cr  <\< 

• o 
^ » 


X 

UJ 

o 

z 


r>t 

• o 

y • w 

Od  O z 

< 

a 

C < 

UJ 

UJ 

> 

Ql 

UJ  X UJ  rg 

• < << 

C UJ  --  o 

o 

UJ 

a ^ 

uJK  cr 

o 

cr 

• ••^c:5q:  I 

I 

— 

Z| 

or  • X 

> - rg 

UJ  X r-  •♦■  u - 

CM 

CD 

z ^ 

K 'DO  u. 

z 

UJ 

X V UJ  K 

K 

K 

UJU.  K H 

X 2 X z 

cr  OO  • 

♦ 

X 

• V 

0 2.  2 UJ 

z 

X O Vk  o 

u 

I 

X • w' 

• Z X M 

u.  • ro  o -u  — 

D 

J - 

JO  cr 

> z 

^ w <►-  Z 

u. 

< 

' 

vl 

o a ui  z c. 

rum 

m in  • o 1 •-*■ 

X 

Z 

>N 

a zer 

UJ 

2 

V *ZO  UJ 

o 

cr 

1 

oin  # — o 

y o 

-‘U.>-aDo<  *o 

w 

• 

2 — 

Z*  OK 

> 

» • 

V • -*UO  -J 

IL 

oi 

O u • 

• • 

oo<-^  *inooK 

J-  win 

in  • 

• • 

• • z 

•- 

(D  • 

V (n  < • UJ  UJ 

z 

UJ 

■^LU  »-  ttX^XO*Nl>  OJa-^ 

Jxoa  •o*--  x-<  • _j 


-^iS'-'U.actsos- 


- w n 


UJ  ITJ  • -X  ►•  < «t  V-  < in  M ^ in  h-  n »•-  o 


XZ  •<  ►-•JUZZN 
C£X4-*V  U.>U.  ’-'•V  >0^»2ak'V<«-'^X 

• ou.  •CfUJftUJ  • — (t—*-'-' 

— -«i/)o  DO*azo  -uifro>'-*iiur<^_j< 
0'“'«fOr.J  K O X > < »- o z o ii  •- u UJ  o >- V M 

o<Mu  — — aaiiin<— >vQcru<«<'-‘ 
uJ>NO  >-••<•  •7J«  o OC£> 
II  K-D>2*-0  K-»*l.>»KO  OUV.‘UJVj« 


^-'fKX(/)X<^»n  *o  ►-  <NJ*- 


-JZ 

II 

OJ 

& > 
« -J 


ocr«j»-joju'2 

^ z 

.J  l/l  UJ  *f  N — r>  O 2 

z?na  Z2»nn 

n > in  y • 

m < J z > 

o^»x»«r>K<j^ 

ffiu.a  «r 

-1  »-  3 

CNJ  rg 

< « J »-  .J  — «r  (I 

*-  n < IT 

rj  < 7 

•f  ii  < « nr 

a •-  _j  — ti 

— (0  < w K 

J D -J  u _j  II  a _i 

U K 2 K U) 

UJ  2 o 

2 2 

UJ  K -f  «i  < II  UJ  c.) fT  n til  o ri  o o n u a o 

O O K n a o 

•»ouuj<j<ha<  Md 

Z < « in  IT 

a -u 

inKirft  O _JU  — tru 

fi  u.  u u.  Q 

U.  z ->  fl  U.  Z Cl  u. 

CJ  Q — u Q U. 

Q a-/ou<o<<o2® 

uu  < 

^ rg  n 

— w < 

-J  7 

^ eg  n <f  m < z 

> 

— m J 

X u o 

r.  m 

•-  rg 

> — 

UJ 

K 

o 

o O 

o 

r>  O 

o 

u 

o 

o o 

o 

O O 

o 

< 

m4  9~* 

•>« 

*•« 

•_< 

r,»U  U 

uu  u 

U o 

c 

c 

u 

140 


L 


1 


— ^ — ^ — -t(vrgrMfvj(NjM(N^fVi 


o^f^r^-<>;r>'r)^-ccO'o^f\jn«u'''Cf-crc^o  — rsin-fj^o 
oooooooooo  — — — — ^ — — o^r\jfvjrgryr\;ry 


K K K ►-  V-  ►- »-  K > 


• ►-  ►-  V-  I 


» 

S3 

o 

» 

O 

u. 

• 

a 

03 

UJ 

• 

o 

< 

z 

w 

o 

•-• 

K 

« 

• 

X 

z 

U/ 

in 

J 

u 

a 

o 

A 

X 

o 

0} 

o 

• 

H 

u 

z»^ 

• 

n • 

D • 

< 

• ^ 

^ N 

H 

I > 

H N 

• 

T • 

»— 

U — 

7 

7.  V 

— * 

Z 1 

O 7 

• 1 

►- 

D V 

— 7 

CAD 

u 

IL  H 

• 

to  • 

z 

J 

O- 

V)  a 

o 

J • 

7 HO 

mm 

V)  1 

< j nj 

7 - - 

nr)  — « 

• o- 

z 

— Uj  — 

7 0 0^ 

• • o — 

y • - 

H UJ 

\ in 

oo  • 

XT.-' 

in 

z cr  o 

— o in 

X • in 

o 

UJ  H 

^ ^ in 

* • w in  — 

X L.  in 

* 

zu. 

w • o in 

>:  7 in  m • 

• » in 

<■* 

0-0 

X 7 . • 

7 7 mm  7 

O u;  i/> 

•-• 

a O ro 

X 7 O -- 

7 i m m 7 

• • 

UJ 

X » <03-* 

X a M u* 

• • m » 7 

O _l  u — 

z 

UJ  D'^  H H J 

X • • w 

•»  a •• 

o < — o 

H ►-  ® 

►-•-•*  K tn 

Z ? jx^  O 

*-*•-<<  X u 

II  U1  X j.  O 

uuj  a i 

u z z > X 

l/)  — U * 5 _J  fTj 

• ^ Q K 

X r>  <J  VI  t 2 7.  # 

»•-  o — z M Stu  — 

u ir  j III  -j  — o • 

(TiOn  7 UJUJ  II  II 

ziA*-oa  ui-i- 
o 


-‘HU.  * 0 a H > 

j < a o X K \ 

uj  T lu  uj  uj  o cn 

uj  ►-  Cl  t cc  fy  o H 

0:  z • »»  •«  ♦ 

c UJ  J X X — H 

X V)  J • Ui  U.I  • 11 

UJ  j j o:  a-  o ^ 

O O — UJ  ♦ -♦  UJ  — 

II  H ■ III  «i  CD  * 

a a I ■ I-  n X 

H y H *>-  II  II  X 

U*  LJ  U-  < U)  u.  X 

a lii  X 

1 

U <\j 


♦ xhv)  O'oaxu) 

ox^l/)  —XujV^ 

• X •!/)  w^-ujac/) 

nx— vAo^rr^^A 

•♦•  — yw  tn</)#o'-‘' 

CD 0 70  tnv»o  «o 

K<7— 

I 'it  w O ^ tvj  M « 

II  nr  V)  — N *-*-  H ft 

oao*  — U.Z 

^ o D a. 

X _j  ij  _iH  j j[T)\rt  JO 
X<o.  <0<*tV)V)<UiZ 

xo*-uouoim/>uffUJ 


J43 


I 


OI 

z\ 

I 

OI 

ai 

<1 

U! 


-•  cj  r>  in  >£;  r-  c 0“  o — ^ t o N a:  o o 


o o o <■>  o rj  o o o o < 
o -«  f j n ^ ir*  c c o 
o o o f ••  O o o o 


* |^j  r)  t oO  o cG  a o ■'■•  <M  f » sj 

. ^ .4  .«  fc  'Xl  f.J  ^i 


O O O o f j o o r,  c-  o o O c»  f » O o o C O 


a)  *nir>(/>  in  t/J  in  i/>  tr>  lo  (/;  vun  m tn  m to  u>  </)  in  «/)  m 'n  lo 
QQanaQ(ir*aanQCiarLaua.aac.Gciaa 


o 


O 

o 

A 

n 

II 

• » 

o 

A m 

H 

UJ 

f i 

» *-« 

O 

<1 

y cr  — 

o 

OI 

n ui  o 

1 

•-•in  rr 

Hi 

'-‘I-  < 

N 

Di 

✓ u z 

01 

• Z OJ 

1 

r» 

ZD  Nj 

Q 

• u.  m • 

—> 

O 

1 

Tsi  7.  7 

o 

• 

NJ 

c>  • m L-  ® o ' ) 

o 

o 

1 

<■>  z - l Ji-  u. 

*■ 

fO 

o 

fU 

C*  D_JH  D Z • 

c 

*^H 

orl 

OU  VJ  JUl^«: 

o 

o 

» -J 

w cn  u *'  :•  * 

• 

H 

•-*  o 

to  C3  > DO 

o 

w «« 

{ 

tOincj  i'»cr 

• 

c; 

_Ji 

n 0 CD  n z rr  < 

o 

— f\J  t?  ^ 

(rt  M o / -)  o > Njfroi.'' 

(1 1 o ui  ii  ii.  o ♦ ••  N - • • « < «r  a 

n ^ 

y\  lo  ►-  T)  w u»  • o M y 0 \ It  fv n in  + 4 

f H • » a / / . I K ® o ‘ c'  ♦ r j ) n n / z' 

M i :zi.)  a tn  tr  i cr  in  » u ♦ -i  < ) ~>  /■ 

I >•  I ) ri  zz"  ti  > j ♦ •*-  «r  o . . «4  ® I # M M o * oi  M If  ti  tr 

|j  I f ) fr  nr  •!  Tf  n h m . « » m > :•  oi  t . « rvj  • . • 1 1 u ii  ■> 

« T : j > H w > II  i-  iM  . II  II  no  04  o ^ / 1 (X 

I II  ! < J - • If  *>  ( » II  •'  u oj  01 (")  It  rr  It  ‘ ) > ft 

i'<  / tn  II  II  / M t i / < Ii  o ■ ) *1.  • « :•  z'  o -c  n > •♦  • « u n a ui 


4\ 

1 ) 

D 

► 4 

• • 

Z > U M 

o 

o c* 

1 i 

» 

f\i  n 

Uti 

t.> 

*■  S 

«r 

.4 

r<t 

144 


!j 


f 


1 


OCOOOOOOOOOOCjOOOoOo 
o^run-tf  iTJ'0h-<ca'O--c\jr)'<yui'£ih-a:; 
OOOOQ,OOcOo  — 

ooooooooooooooooc>oo 

'L2.-Z.  Z~^-Z.7  Z Z‘Z  Z-ZV.  y 7 Z y.  7 T. 

aCLQaaQaciaadaaaaaaaQ. 


CO  ® 
0^0 
* CO  i( 

o o o 

• <f  • o 
in  o in  + 
f\j  • c 
<f  m cn  (I 
O'  nj  <M0 
m fO  ro 
<}  a - -x  CO 
>t  m o ^ 
ro  -vj  m 
4.  to  CO 
j3  m + ro 

a + -o  ■*• 

•It  O C f\j 
o O « o 

• « o « 

Cj  o • o 

o #o  • 

fvjinoo 
•O  C\J  O'  Nf 
CO  ^ 

fvj  cr  in  ' 


<n 

o <M  O'  ri  o 

o 

^ S 5)  1 / 

• 

1 <f  ^ O'^ 

o 

<t  1 1 • (0 

o 

o <f  in  o 

<1 

• 

h c oin  * 

<s 

O * f . o 

•«— 

o 

**  O c « • 

cn 

u 

• 

in  • • c»  in 

1 

♦ 

<\J 

•j  m in  r- 

Kl 

cn 

o 

m o 

ri 

» 

• 

U-)  ^ rn  1 O' 

ai 

N/ 

<t  m in  f\.' 

:;i 

o 

ox  r)  — N 

^1 

> 

• 

4 ^ n oj  <t 

1 

• 

(/» 

nj  h ♦ o m 

VI 

O X 

• 

O <V  <M  « 'O 

1 

o -- 

h u 

# o o cc  ♦ 

c o 

7 

o » h o -n 

ai 

o ^ 

UJ  N» 

• o o H CJ 

-tf’i 

y o 

7 

c>  • #1  o « 

Ul 

•-<  >-* 

D > 

'D  o m • o 

1 

viin 

O 

aj  O'  cruo  a 

.M 

a Q 

tr  X 

fn  o <T  ^ o 

CJl 

II 

< V 

-«  1 wO  o 

HI 

O Jj 

.oc 

1 o 1 cc- 

hi 

lij  7 

K • • ^ O 

O a O CM 

Z\ 

(/» *-• 

• I'J  — ^ •+  -■*  (M  <M 

t o -•  a O' 

o\ 

• H •o.j'sj'Mxnoooc3in-'*^<ra) 

^»l 

X 3 O' 

— Qw^MC««  *• 

a 't  « — ' in  -•« 

1 

^ w u 'll  'H  O'  1 

(il 

ocry)(no*~'On-«oooo«  iini  i 

Dl 

(T  n)  Q 

a u co  u > H II  II  II 

M '/)  1 1 

\ 

11  D 

O — CM  O <r-  II  M 1 

OI 

7 m II 

II  II  II  II  u o c o O o X V 1 

71 

•“i 

•-•  X 

Y 

z 

s 

c 

hi 

h 

inl 

o 

n 

•t 

J 

XI  u u u 

OI 

71 

I 

CM 

crl 

<l 


• <7  ir.  <.  t~  o'  ?*o 

« ^ ^ ^ ^ rj 


* o o o r«  Tj  o <■*  o o o 
o '-*  cj  m <1  »Ti  -o  N a,  a a '-•  (Vi  iO  o cj  o- 

c.  o C)  O o C-  o o o O . 4 -4  w 

o o o O o ' ' < ' o c o o o f./ ' > o r*  o o o 

3C  ^ .-c  <<*<<•<<<  i jc 

or  cr  (r  a r/  nr  fi  rr  fr  fr  cr  rr  tr  r:  G.  O'  tr  cr  a (j: 

o o o o c f j o o n o ( * D r j n n n o o o 

lLlLll.UliLU.llU.  UL.U.U.UILUU.I1LLU. 


a 

L> 


<1 

(/)  • li 
fTj  0 0 

<1 

— fi  7 

*— 

— 

OI 

n • 'on  — 

CM 

1 

7 r 7 r 

1 

4- 

Hi 

*.4  7 r*  II  *-<  U. 

»-« 

“> 

D* 

-CJ*-* 

■*- 

w 

Cl 

7 » - (/)  .-4  li. 

0 

0 

71 

• loa  7 

1 

1 

►-I 

k;  3 r>  »•  111 

0 

N 

NJ 

\ 

►DU  </)  7 > 

f\i 

V 

V 

VI 

C>  0 U tj  till*-,  d 

— > 

1 

0 iiJ  Cf  D > U H 

0 

-4 

C)i 

0 £’)  a 0 j 7 '^H 

1 

* 

Ql 

0 <1.  V)  C u UJ 

• 

»-4 

**> 

<1 

^ < 1/)  0 > .^4  > 

<-* 

n 

w 

orrtorrii  -id^o 

c 

0 

1 

u n ./  •!  n 7 fvi  7 u 

N 

ft 

.n 

U li  «f  •<  0 0 

(.1 

—s 

C.'l 

II  5 y 7 • < 0 f?  < 

..4 

7 '“4 

7 

n 

‘X! 

i}  til  0 u • jy- 

• 

► 1 

•4 

i 

III 7 0 <c 4- 1 » no  »*' 

.. 

n 

Hi 

n 

~> 

71 

l/)  -4  U .T)  10  7 7 Ill 

» 

I>I 

M r 

II 

ft 

IJI 

• H 7 X -t  1 

• • 

►4  <M 

<0 

4-4 

4-4 

rvj 

Ul 

X D ••  ••  111 » ri 

lilO  '0 

■7  n 

* 

1 

>.  0 0 4.4  > ^ 

J 

7 0 0 II 

ir  fr  0 

*■  4 

II  cr 

U 1 

(Jc  w II  -4  7 on 

0 

..4  — 

:■>  '1 

n 

1 

*) 

OI 

0 0 II  II  0 ' — u > 

w -t 

» 4 

1 < 

n 1 

n 

1 

M ‘ > — 0 0 

nil  « 

rj  — 

III  V 

( ^ 

II 

‘-‘•U  7 

111 

7 t/»  f)  n II  4 4 II 

u 

•■4  O’ 

00 

ir  U 

1 ) 

~)  0 i t 

‘Ll 

/I 

l.»  7 7 

•( 

‘«l 

• 4 7 0 M 

u 

0 

CIj  0 

0 

» 1 

H 

«-4 

n 

(/fl 

Li 

<1 

ra  u u u U LI  i.) 

u 

u 

146 


r 


OI 

'Z\ 


a» 

<] 

m 


V - r.  oOoo<^ooriooOooooooooooooc'Oooooooooooo 

f\j  r»  <*  iH  vo  ^ a;  cr  o — (\i  r»  sj  uo  o N cr-  a o — rr  IT)  c N X c o — (\j  ’O  O'  o — f>J  r>  s;  to  C N X O'  o — rv  N- 


o C*  c o o O O o o O C.I  c o O C ■ o o c c 


ryfvjfNj  rj  r\,  rj  rvj  (V  fvj  cnj  nf’' fO  f*r^.  nr. nr/. ^ c <J  <t  <f  <»<x  < >Jinir>u)ir>inuiintn 
. - C O c r.  O O C- O O <5  O o f Cj  c o o r- o o <*>  o o o o o O o o O o C- o o O O o O o 
JDPJ*_?D:^D^Dl’3DP.).“):>JD:)3DDDI>DDDD3Dr>DDDr>DppDp 


Ui  UJ  Ui  Ui  U>  Ijl  Ui  UJ  UJ  U U>  Lij  u L)M  UiUJ  lilUt  >jJ  Vi>  11  LJ  Uj  tiJ  UUi.'^i  U1  ui  U1  U ' UJ  Ui  Ut  111  U'  Ul  Ut  UJ  U1  UUll  Ui  Ui  UlUJ  UiUI  Ui  U1  UJ  Ui  U 

a a cr  Q cr  c/  (r  cr  ct  cr  .'z  0 0 i"'  ‘I  fi  a fT  (f  cr  c/  a cr  L'  ft  ti  (r  cr  O'  ft  fr  Cl  rr  fr  c a ir  (r  tr  cr  cr  a 0.  cr  .T  o 0 or  Q cr  c:  a cr  c a 


y 

u 


yj 

7 

U. 


o 

z 


to 

7. 

o 

K 

U 

z 

D 

U. 


Q 

U 

ft 

U 

u. 


X 

UJ 

a 

> 

o 


CL)  C'  - -- 
• — ^ O 

^ — o 

_J  <I  — w 
< ^ w < 

. ? Z X 

z to  r\j 


QJ 


7 

< 


> ^ 
lO  o 


< 

to 


y 

to 

d 

ft 


to 

Q 


< <i 
NJ  rs. 

00 

ft  « 

> > 
<1.  < 

> y 
to  to 
« I 

1 -* 
I y 

< to 
— n 

♦ y ft 

u;  to  n 

0 • O k; 

UJ  nO  O fJ  ft  « 

1.  U V . 0 < *- 

Q < rr  i‘J  > 

r ft  ft  X ^ w 


J 

UJ 

> 

< 

^ 

- * 



> 

X 

> 7 

7 

X > 

cC  > 

X 

n f*-  c s.  N 

UJ 

H 

• 

r 

a Q a 

a n 

— * 

-< 

to 

tO 

ftw  r.1  V 

to^ 

to 

^ O u.  u.  CO 

to 

r 

'< 

7 

> > > 

> > 

y y y 

y 

y 

< o *5  n o r.  < 

r 1 

r.  n o Ui 

< tn  < 

< 

<1 

to 

< 

1 

fv  r> 

w w 

> > > 

> 

> 

CM 

ft 

OJ  ft 

ft 

i o; 

ft  CM 

ft  « 

ft  K 

1 1 U I to  to 

UJ 

u 

T) 

y 

z 

-f  «r 

< 7 

w w 

w- 

> 

*— 

> 

y 

Ul  lA  a r* 

n 

o 

ift 

• 

» 

• 

N.. 

e-  f»j  r; 

CM  CxJ 

< «T  X 

/ 

to  > 

X > 

> 

y»  to 

> X 

>>>GU.U  0 5 

ft 

rj 

1 

1 

cr 

< 

y' 

o 

y y V 

y y 

— CM  CM 

01  OJ 

« 

— 

« w 

W ft 

w ft 

w 

tJ 

< ca  w 

or 

» 

rj 

7 

f 

• 

— 

to  to  f/. 

> X 

y y "• 

y 

: 

CM 

o;  / 

7 

X 

.*  -■• 

X 

» ft  > 

> 

►- 1 

1^- 

UJ 

D 

1 1 1 

... 

1 I 

to  (.*:  to 

— 

X 

X 

; 

CM 

y CM  rj 

o;  y 

CM  y 

CM  eg  rj  r- 

»-  - • cr  — ' 

< • 

o 

• 

» 

Q 

«i  .f  z 

n 

< X 

ft  ft  ft 

» ♦ 

ft 

4- 

1 

y 

1 :■ 

y 

y 1 

y 1 

y y 

y C) 

u u.  ( ' ** 

ai 

u > rr 

1 

< 

<T 

> 

^ ro  . 

> 

fj 

y 

«t  < X 

< 

X 

x; 

- )0 

X 

t/) 

X 

to  > 

tOL 

«<  23 

— 

7\ 

u 

u 

• 

— 

ro 

*— 

— y y : 

.... 

-• 

y > 

> 

-*  CM  C-J  > 

•«- 

C'J  OJ 

> 

ft 

> ft 

ft 

•2  > 

1 > 

1 1 

1 

fc-  »-•  ^ 

y 

— a 

u 

y 

y 

/ 

z a V \ % 

z 

N S 

_ 

... 

: > y 

— ■ 

r~. 

y 

y 

OJ 

ro 

o.. 

CM  w 

...  ... 

— 

to 

\ N to  tO  to 

1 

0 UJ  X 

• 

• 

Z X 

.Wtf. 

tf« 

^ — 

X 

... 

7 

; 

y 

X 

\ 

< 

y 

*t : 

y 

y < 

> *1 

> 

»- 

I--  ft 

1 

M 

O > OH 

UJ 

o. 

to 

r, 

1 

->  > > 

1 

w 

> > 

to  > 

>. 

> 

.— 

X 

1 

*-  1 

1 

1 - 

. . 

7 

si  c 7 - 

< 

O < fJ* 

1 

• 

• 

y 

> 

r> 

z 

< ^ ^ 

X 

< 

ft 

w- 

> > > 

ft 

w 

> 

> 

> 

X 

y 

> 

y 

7 y 

y 7 

7 f..' 

u u UJ  y 

— • 

cJ 

o nr  H CJ 

V 

S 

> 

• 

C 

N. 

tr.  < -f  / 

\ 

X 

< 7 

X 

<s 

•W.  ....  _ 

X 

w 

w- 

-- 

f' 

to 

•> 

to  > 

lO  tO  to  to  NJ 

isf  — 

< r 1 

y 

trl 

u cr  ij 

LJ 

» 

— 

<-• 

h- 

1 — CM  ri 

— 

1 

n t>j 

\ to  < < X 

s 

X 

«c 

/ 

A 

ft 

ft 

ft 

ft  ft 

■ft  ft 

ft  U 

ft  ft  u < 

to 

D .ii- 

UJ  f/ 

'(l 

y 

o 

> 

<■  y y y 

> 

•t 

y y 

— 

ft 

-*  CM  OJ 

.... 

ft 

Oi 

CM 

X 

> 

7 

/ 

7 

^r'  ^ ^ 

o 

ul 

u t>  tJ 

J 

Ul 

« 

0. 

» 

-« 

D 

— * 

V to  lO  to 

■w 

N.  X X 

> 

< 

y y y 

•f 

> 

ft 

ft 

... 

to 

if. 

CM 

n; 

U 

u u u y 

ft 

i 

ujuito 

ft 

1 

-/ 

> 

•— 

O 

/ 

ft  ft  ft 

/ 

ft  ft 

w 

\.  to  U/  to 

s 

X 

X 

► . 

«T 

*y 

ft 

« 

•ft 

ft- 

u. 

< cu.  to 

J 

rr  cr  UJ  to 

n 

4i 

► 

N 7 

CO  > 

X 

> 

— 

—• 

ft  ft  ft 

/ 

ft 

ft 

ft 

X' 

V' 

«« 

7 < 

X X 

XU) 

I 1 L-'O 

> 

a 

M H UJ 

to  lO 

*; ' 

> 

r tr. 

— 

ft 

--  a a Q 

ft 

a n 

to  > 

X 

A 

A 

A 

> 

> 

fl 

n O' 

CM 

0.’  CM’  (M  CM  ri  r.  r 

<M  0»  C ft 

rrJ 

O UJ  ~)  > 

Ul 

f£ 

< 

UJ 

o 

<>>> 

< 

> > 

♦» 

w 

I 1 1 

ft 

W. 

1 

1 

> 

X 

•— 

ft 

y 

J. 

y 

7 y 

: y 

:•  u L.  u u — 

< 

HI 

u;  i'  n — 

UJ 

• <j  cr 

• 

• 

a 

to  II  M M 

n 

X 

i>  II 

> < 

II  It  11 

> < 

II 

11 

— ■ 

to  ft 

to  X 

>■ 

X to  X X t''  X tr. 

•T  X > 

to 

Zl 

to  - 5 >- 

DOT  ft 

ft 

• 

> 

« 

> 

ft 

1 

to 

1 

X 

X 

n 

<— 

ft 

r.  CM  n c j CM  o;  n n r> 

r r r.  to  ft  ft  o w 

•ft 

'■.I 

• H Cj 

c 

rr 

o r>o 

-* 

H 

-- 

II 

— . 

n, 

... 

n 

« 

II 

ft 

to 

II 

> < ft 

y 

ft  y 

y 

y ft 

ft  ft 

ft  ft 

ty  7 

- X ^ 

u 

1 :')U> 

r- 

U U X U17 

— 

UJ 

ft 

a 

a 

1 

> 

II 

w. 

to 

— 

V 

« to  X cO  Oj  CM  Ol  01  M O' 

u u *-  to 

y 

7 

1 

a »*• 

( 

ft  OCJ 

w 

• 

> 

;*  -J  < z 

> 

< 

C) 

L» 

1 

< X 

1 

•I 

•t 

c 

y 

: o r n :■ 

y y. 

y y 

: rr 

r.  a — 

«- 

n 

ul 

o cr  o cr 

t i 

r 

■ 

} u> 

> 

< 

/ 

> 

u 

h ^ n <\j 

I) 

r>  f .j  V- 

II 

H 

— rj  CM 

i; 

fy  CM 

tr, 

II 

II 

II 

II  II 

II 

II  II 

II  II 

H II 

II  u. 

<—  UJ  II 

M 

D 

C)J 

T:  f.  rr  u 1 

7 

n 

•r  ►- 

y 

fVj  fj 

w 

— 

< y y y 

/ 

•r 

y \ 

•T 

y y y 

/ 

y 

X 

it 

.. 

OJ  n >i 

to  7.  ■ - 

oi  n 

<f  iOnTk 

II  1'  »-•  X 

J 

V-  O 

» 

M ;if « 

( 1 

> 

UJ  /* 

Q 

y 

> 

u 

a u. 

lO  l/j  l/J  OO 

H 

X 

> V 

» ) 

y 

tot/- 

to  t/j  (0 

X 

X 

X 

X 

-r' 

fx 

«< 

ftj  u 

u 

li  u 

u 

u u 

u u 

L U 

u. » 

X X 1 J 

_J  U’  .' 

Zl 


•t 


tj  tJ  O 


lO 

n 


h 

n 


u 

to 


u 

to 


147 


3 


1 


2t 


cri 

<1 

ui 


— — — ——  ——  — — — C^Cv.Cg 


0 f\  n ^ L*^  o S ^ c-  O — (V  T'  <?  iT  r ^ o ^ 

OOOCCOOOOO  — — ^ — — — — — 
OOOOOOqOOOOOOOOCOOC-Ooo 
in  LT,  t/)  tn  to  i/^  CO  C-'.  tTi  (.'“■(/)  l/J  CO  (/5  cn  cot/)  lO  </)  ;/)</> 

1 5 T > ? 5 5 5 1-  > ? >-  r > :•  > > > - y >• 


to  CO  to  to  I/)  CO  to  to  to  to  CO  CO  yo  CO  to  to  m to  CO  t/)  cT-'  to 
or  z 

UJ 


» 

C. 

u 

(VJ 

to 

3 

O 

O 

UJ 

7 

o 

lT 

C 

7 

o 

U 

< • 

u 

Lu 

> 

^ ft 

o 

3 

O tit 

UJ 

H 

X 

H 

U 

- a 

<1 

X 

to  to 

o 

X 

U) 

u 

<r 

> o 

7 

Ul 

1'  z 

3 

— 

Q 

— 

D • - 

s 

V,'  fc- 

r 

'1) 

O 

X 

X U 

7 

• 

7 

• 

2 .J 

rv 

< 

re 

— ^ c 

— 

— 

< V 2 

a 

UJ 

Cj 

UJ 

O—  • 3 

tu 

» 

Li 

< lit  u 

H 

• 

H 

• 

CO  rr 

< 

<" 

'<  LJ  >• 

O 

Q 

c 

-hCX  j 
7 nan 

u 

H 

u 

H 

* '/,  to  u 

z 

- 

7. 

• 

Xu,'—  u; 

L. 

» 

X 

• 

< L'  7 U 'X 

u. 

X 

tL 

r ■*:  — u rj 

3t 

• 

■< 

• 

X H C UJ 

H 

A. 

H 

(V 

<t 

* LJ  X u C 

— X’ 

— 

— X 

70'- 

c- 

X c 

X 3 

Ul 

<1 

— < n II 

to,  • 

• 

< 

rj 

2 7 X 

4.  C 

2 Z 

- 

j 7 

• 

1 

i U OjCVi  — 

Co  c 

X c 

X 3 

— 

• — i u c 

— 

^ — 

n 

V — 

II 

I?! 

<r  H 3 — 

« 

< k 

< f- 

C.I 

♦ «f  O'  ••  » <3’ 

— V ^ 

• < 

* <i  • 

3 

21 

— tL  X • 

— — 

— — 

jy 

— — 

C 

— ; 

u > fi  a 

^ — 

z > 

« 

7 > I 

» 

1 

♦ UJ  L>  X 

W 4- 

*-•  X 

X 

►-*  JJ  — 

X 

SI 

n C 7 T 7 — 

s-  7 

> o 

res 

1 

O'.:)  c r. 

V c: 

V 

UJ 

X 't 

UJ 

n 

l.  '0  CO  u 

< s r 

w .NT 

u 

w • 

to 

Li 

o c to  3 

r-  — X 

<s  • 

z 

< re  7 

<\ 

to  ri  2 c 3 

*'  — 1 

* c' 

< 

• -0- 

< 

ul 

> 2 cr  u.  u<  t*  1 — 

► V 

>f>  r T 

C.' 

to.  ;•  u O' 

1 

3 3 '7  I--  ►-  - 

7 7 

• CI  — 

. r>  * 

.ji 

to  c'*  co  ;*  «r  ^ — 

— ''j ; » 

- S liJ 

u. 

f 1 

n u 'T  ^ r 0 

' ? V 

<•  • 0 

X 

'■>  • 

r 

Cl 

r x>-  rt  uu 

y 1 -a 

f -•  • 

H 

O II 

► 

Hi 

y.  7 3 7 

n ^ s • 

- I OJ 

— X 

y\ 

co-o  IJ  n u G 

- — vo  - 

„«  — 

n 

'*UJ 

■» 

Cl 

• H ^ T •■■• 

^ — w 

— UJ 

It.’ 

— fT  UJ 

Ul 

T 3 G H - - CO 

o — ► 

V ' 

> 

»-  UJ 

> 7 

1 

H '1  G t?  O 7 C 

r.  o r<  n H -I  - 

C f* 

.-•11 

( i D 

U 1 

GC:  ff  UJ  • 

ro  — tM- 

> 

> ^ 

7 ^ r. 

* 3 

i.A 

n : ■ <■)  UJ  > o 

-to 

•-•  a II 

H 

— ii  ir 

H C 

1 

H 3 t/)  M 

U II  II  c 

O'  3 

uj  a tj 

Ui7 

C.'l 

.’■•/)-  o - r-j  c i/i  c t'".  c r.  u 0. 

a X 

cr  UJ 

yt 

t-'  3 t'  II 

fj 

(V 

-1 

- to  7.  7. 

c 

— K. 

re  T 

HI 

H 

o 

O 

o 

toi 

L> 

rt) 

o 

o 

-1 

< 

— 

-rt 

ru  u u u u 

148 


— — — ^ — — -.(\,rg<\ifg 


OOOOOOOOOOOOOOOOOOOOOOO 

o-'C^n-ijiTvC^^-c^o  — cv.r:<juT\0Ni:c'O'^?v. 
oooooocooo  — — — — — — — — — '^rvjCs.fvj 
OOOOr>OOOoOOOOOC-OOOOf^OC>C/ 

<<<<<»?<<<<  *t  <<?<<<<<< 

ooC)cjooouoocc:aac)CjCocoucG 


* » 

* » 

X 

X 

I m 

I U) 

•_■  » 

— » 

s o 

''V  ^ 

m * 

m • 

— r) 

— r> 

_)  — 

J-- 

• 

S UJ 

i U’ 

in  • 

« • 

P * 

UJ>- 

— • 

r\j  • 

ua 

fV 

S’ 

S 

— UJ 

j; 

Ul  II 

U.’  II 

ou 

ID 

U) 

tn 

/f  ft 

IT; 

— 

r\j 

— UJ 

» 

♦ 5 

* 3 

H 

rg 

UJ 

UJ 

UJ7 

? 

— * 

f\.'  • 

> •“ 

— UJ 

Z[  • 

3 - 

X • 

UJ  X 

U'  X 

K-O 

1 <V 

IT. 

in 

— 

0-. 

uz 

L 3 

* • 

• • 

.J 

_J 

< 

wq: 

_JC 

* 

:< 

a 

O O') 

•<  • 

t • 

«■ 

« 

U.  UJ 

• 

* X? 

« X 

rg 

>- 

u.  a 

« CM 

— — 

f\J  — 

3 

3 

• X >: 

S U- 

3 ID 

ID 

U 

_j 

X UJ  (V 

X •' 

cr  • 

'Ji 

in 

<1 

e m 

ID  _J  • 

in* 

m • 

♦ 

4- 

»-l 

z 

O • G 

f\J 

— 

<1 

Z X 3 > 

* II 

•f  II 

3 

y 

a 

ho 

— XJ  O UJ 

U 

ID 

1 

< or 

G U 

— — 

fV  (\J 

• 

• 

c u. 

UJ  Z • • 

3 3 

3 S 

_J 

J 

d; 

u 

> - — - 

xcr  * 

ft.  ft 

3 

Cll 

z 

— 3 

w to 

— « 

iil 

X c 

X UJ  1 UJ 

- I 

• 

— (\) 

-J  — 

H 

z - 

u > u - 

II  • — 

II  . — 

J 3 

y 

— < 

X N 

X \ 

3 ft 

• ft 

oz  X d:co 

g^LjU  u-Uisr^s  *01  -vT  cr-*.  bj-*^ 

cojcr  ^-rvj  ft  — 

o 5 c '/)  a ft  •-'  • • • f’*  • r»  <>-  u 3 

-J  h-  U — u.  - r.  X r>  - n'. O o:  u a 

<<  0*0u<ii5— — ^-U'  — UJ  U—  U.W 

oov  r/ajcrovii)-  u»*  u-x 

L.UUJC*-—  ..  w_j 

M H < -.Oh  IMU!- X s • . ^ fl  > a 

c u;  D u I < cr  a o u u 3'  uj>- 
u r.  z UJ  c ("  c — r:  M ii  < u ft  u 

— rvj  cn  no 

•HO  Jill  — ww,--  f\J  5 UM  Cm 

X DU  > H n u H » - : ut  ^ jj  ;•  cr  7 ft  >*■ 

f-O  x<  — <<cto  ftiAtrVHTXUJcr 

c)^fvu><.rr  m?:-*-  ••Dft<DftnD 
cnfp3j  »-wOv»~t>v-0 

I!  _)  H ft  II  < : aj  <■: u:  7 lj  uj  / o u*  7 

7 i/^  c > o c-  - • O'  u.  u ;r  UJ  u a iiJ u.  cr  bj 

u — njr’<‘iT' 

•-  \jO  ^ \Lh  — M 

» oo 

C -D 


O O o r.  o ''i  O O O O o o O O C;  o o r;  O CJ  o o o o o o o 

o — ^ c <0  'J'  o — fvj  r»  c iT  N ^ cr  O — Cvj  ' -t  i/' 

o o o o o <v  c>  o c — f'j  fj  fj  f.j  rj  cv. 


O o f i <u  o O c*  C > <■>  o o o c o 


<<<<<<f<<<<<<j<.t<%'<r<^-r  •!<<<:<<< 

OOOOCCOUCCQCjrjCCitJCOOCLJCCJCjOQQ 


1 


•C  X 

\ l/l 

a - 

> 0 

• • 

< • 

-J 

X 

u n 

X ui 

i.  — 

s 

^ » 

— UJ 

a 

SnO 

* 

z 

• • 

+ 

< 

r> 

0 

J - 

^ n 

« 

• • 

:<  uj 

(/>  H 

« » 

^ 0 

c 

til  U 

•™  > 

,<  _j 

• 

U Ui 

a a > 

• <5 

*— 

— V 

> i U'  II 

.J  - 

* 

O Lf 

< <r  (Tj 

< • 

0 

i'  “1 

a L'j  — 

X 

• 

— H 

• • > 

N ifl 

0 

z 

o o u 

0 • 

•jj  •" 

-J  _J--- 

• -0 

4- 

> 

< > • 

• 

►-«  o • 

J • • Ul  > 

1 n 

H / Ul 

1 f:>  O tfy 

J 

o < ^ 

U.  t/3  10  . . 

(j  u ; 

•«  -J 

< 

--a  a .1  vc 

.J  • 

« < 

C£  UJ  _J 

C Ul  U;  .<  • 

-:  - 

J * 

u.  a 

^ • » fl  n 

« 

^ — 

liJO  z 

i y z 

C-  II 

w > 

Ct  .JO 

X — >.  LJ 

_l 

UJ 

X 

w .X  Cf  a • 

< 0 

a iTj 

<\ 

<.0  H l; 

. 

^ la 

• ^ 

Hi 

^ 5 G O 

. u — — 

\ a 

^ — 

*3\ 

— o a z 

X ;■  7 3-  •*•  11 

C'  UJ 

1 >■ 

Ol 

H cr  a *•* 

1 J c:  UJ  UJ 

• - 

CM  Ul 

1 

< u oo 

C.  U fi  IT,  — — 

0 

Hi 

O T liJiT**  •>> 

w X 

UJ  .J 

3.1 

U /^  < > 

— a a • 

w O') 

X 

L\ 

o 

N z Z ^ IG 

N • 

— 

C « 

0^7  H 

uj  1 Ul  u)  • X c a 

0 

^ ^ _j  — 

Z y.  <r  o > u • • jj  • - 

y;  • ^ 

J 

-J  --  - 

• 1 

< a — < 

- • - - X \ 

a n\ 

.■5 

< \.  y • a 

Vl 

o y u;  M a 

►-  Ul  Z — uo  N. 

lij  \ 

B 

- 0 0 a ta 

1 

O-'U'C.  t-  L 

1 a cr  T'  • o 

11  >4 

0 

aon_i  U14- 

r.i 

o f j 'r  «/>  'i> 

<«i’u*)Xuo»  a • 

* - • 

.1 

»-  a Ul  a — 

QJ 

o :?  O -1  a 

^ . • p-\ 

- r- 

< 

< UJ  a 13  z 

<1 

H ► u r:)  ^ 

H •-  r,  X 

rvj 

w 

u iL  a 

Ul 

< <f  Cl  • 

IJ  B 3.  -.7.  — — UJ 

» II  UJ  V 

( J / w u_  ^ 

1 

Cj  V t/"i  fj'  j I. 

(T  -J  a f - IT  V IJ  • 

B r 

0 

y ► w y 

U tJ  < H — 

4 V • V \ ^ 

r. 

--  a fr  - J 

'•} 

M H HI  •C  1*3  I- 

UJ  {..  — X r-  V.  - 



•«• 

B >.  n 

U! 

o Ui  r;  »r 

-r  (1  OH  c M 

(r  nr  u 

h y LO  L!  y 

Hi 

U / n (A  UI  r , 

O 0 — 0 T |< 

0 

<»*n  ecu 

y\ 

t/i  — :-  -J  -I . I u • — c - - 

II  - a 

• 

n C'  t.J  0 n 

c-l 

*y  (j  7 'll 

w-  w w .<  > 

> 

11  .1  u 11 

J >u»o  > u 

II  LJ  H 1 :•  Ul 

[ *. 

II  .« 

y (c  y (17 

1 

H ( . 'f  <f 

• *M|  1/3 

> ■ • 0<^(T 

> ( >^  n > liitr 

Ij  1 

M fr  r 1 / « <r 

I ri  7 ► 

• \ » 

L3  J 

n 1 •!  -,  'T  fr  7) 

< Jl 

n ( » > ►-*  » 

!(•  f>  -in  ► • fr 

\ 

llIH 

1 lO  1 1 » H tj  1 

1 

M . M III  ► II  *f 

r : 1..  1 : (V  ( J 

:•  Ul 

n t * 11/  / u Ul 

K'A 

(/,  -t  n •!  f 1 

--  u ij  0 1. 

^ r; 

Ul  u U.  (I  UlU.  CL 

/I 

U ( ( f -J  WU  — CM  r.  .|  IT, 

• • U t.3  h 

ll.  -•  cu 

» 

0 c. 

vJ 

0 Ci 

< 

f,t  f^JUUUUi.J 

150 


OI 

3 

a 


ol 

! 

tt| 

y 


— - — fw(Nj(\jfv.'(\f\iCvcw(\*(vir)nnr‘>fnr)ror)r)rx^^^<^c^«<tf)ir)U’/iDtf)ipL')(/)if) 


ct 

o • 

• IL  o 

V)  z 

4UK  •- 

ua  K 

hii 

f^UO 
Z C£  U 
•-•Ui 

»-  u. 

ILIZO 

> •- 

^ Ct 

►-OZ 

< < CD 

a 

u.  UJ  Z 

uja  o 
ao- 
_/»- 
ou;a 
z a 

« I o 
K o to 

or  m 
o u. 
u 

ZliJ 
Q OH 
Z ^UJ 
— < (/)  J 

o z a a 

O'-'UJ  UJ  > 

otvjcr  a *3 

O z O 

< < o a 

00(0  o 
(jvjuja  u. 

II  I-  < 

0 in  O UiO 

UJ  z Q z o 

,/)«  7 ^ y 

• I-  n JO 

1 OU  CJ 
»“  C)  X J 
Oa  Ojlil 

ino a jui 

II  D ►-  0 J 
7 Mi  <->  C 
O OM< 
^ UU>~ 


• o 

X Z 
UJ  - 
O •!- 
Z X < 
i-iUJO 

ou 
UJ  z 
> « cr 
« o 

HliJU. 

u > 

< “•  X 
CC  M u 
u.  uo 

UJ  < Z 


•OUJU. 

I Z (Y  o 

I-  *- 

O H UJ  K 
7 <(f(t 
UJ  O Q < 
-lUUO 
UJ 

> II  M J 

4 4 

* n:  rc  UJ 
►-  UJCt 

II  -r  cr 
o o M 

. JUU 

^ u.  u.  < 


< 

o 

u 

a 

• Q 

ou.  ' 

z c 

■-•  X 

^ UJ  - 

• <o- 
ooz  <« 
z u*--  < 

H U-  LJ  « 

4 U > L 

o -c 

u X I- 
UJ  O'- 
U.OC 

oza 

-U.  L 

X UJ*> 

uju.  q:  V 

oo  • 

Z Ui 
-I-  X 

cri- 
u.  < 
oa  u. 
o 

I-  >- 
orai- 
<<a 
a z < 

JO 
< <>- 
UJ  a (t 

a"-< 

z 

ft  a: 

ooo 

u.  u.  < 

a 

10  to  — 

UJUI 
— -UJ 
IJO  I 
Z Z I- 
UjUJ 
o r>  to 
O 13  1 
UlUJ  Z 
ft  ft  — 
u.  u.  a 

11  II  II 

©uo 


• N 

o no 

— * o 

♦ -'■o 

X o n U) 

UJ  o » 
*a. 

0 . Q n C4 

r » 

0*^C4  • 

I.JO  — in 

w n <r 
NJ  U • N 

1 

o*^cr  • 

♦ o o © 
UJ  O'  • 
w wfO'n 
CO  © • »\ 

♦ © 
J «-*  N • 
4 o o ^ 
UJ  o • » 
II  — 1^»0 

< • n 

z • • 

►-  o n 

• 

w UJ  -r  m 
Z ui 
. j III  < — • 
Cl  > • 

a rf  <r 

•-CC* 


O'-oo  — -•ooo'-n© 

I I I I 4 I I I I I I 

UJ  UJ  UJ  UJ  UJ  U.' UJ  UJ  Ui  UJ  UJ  UJ 
o®ncjCT!--  — <»  liO  — 

cya-ooiniOTin—o  — -* 
ZmJ  iJ4n^^f^JNO^•— o 

f\jCM(/  — — 

0 — — o — o ooo  — no 

1 I 4 I I I I 1 I I I I 

UJ  UJ  UJ  UJ  UJ  UJ  U UJ  iU  © UJ  UJ 

ro  ro  K o s ru  -s  o m u“i  o ft 
fMiricrotoooNnooin 

<JOOCOC0<?O(\lfMOtf) 

— f\jnft‘*t<jr\jNcnfu 

— oj  — o — o o o o — n 


♦ I I 


I I I 


<3  n 
n <3j 


0 — 0 
ft  < © 

nn  o 


o N © © 

iT)®  r-  fwN 

o • « • 

• *o  o n 

o — c\j  ow 


t\jrj  — 

4 o 


<NO  ® 
CNJM 


N ('«  O 
CM  K O 
. . ft. 

N m • 

- ^ ft 


l/)0 
^ — 
in  <t 


® CM  o 

N O «T 

® © o 


N m r-  rg 


— — o n 
n — • eg  ft 
OftO  • • 

• • s n 

— O—  cg(M 


n rg  n 

N N — 


OftO 
ri  N (M 


• -40ft  — — 

— ooo®«j  — ®fg 

• on— soo—  • 

OCVjCMCM  — — — — S 
ft 

• u*i«t<tsn— 

o — — nftKn  — n 


(M  o oj  s in  o 

(Njfg  ^ 


oon— soincvj  • 
OOJtMCM  — — — — 


UJUJ  UJUJ  lu  UJ  UJ  UJ  I*  MW  UJIW 

ocMN-ftCM'a^oa'rginino 
Non®in®oo®ocMvo 
<M®c— oinnc<rocM® 

4ft—n  — o®ecMft®K 

— — OOCMOOOOOCMin 

I I I I 1 I I I I I I I 

UJ  UJ  UJ  UJ  U UJ  UJ  UJ  UJ  U.'  UJ  UJ 
iT-ftnO^J  ftft®CMOCM<J 

• otn®  — Nnoin— 

n <j  ft  c ® fu  c ® o s - — n 

I 

ujcMCMnn4no<jcM  — — — 
np**^*»**»*-f 
cMn  — oo  — ooooocmo 
ft  I I I 1 I I I I I I I I 

• UJ  UJ  UJ  Ui  UJ  UJ  UJ  lu  itj  UJ  u;  UJ 
OSCMfvjrgrgcMON  a:^“  — i 

• ®o<jN— oo<J'no<pn 
nfMONnoor*)CMft<j®n 


no  — ^ — CJS  <3  CM  — — CM 


ooo 

IP  • 
O N O 



OOON  M MOrg 


O O ® CM  O C 
4 N K.  CM  rg  I. 
ft  ® ®c  O' 


— n fMn 
® CM  o o 


CM  n o in 
CM  4 ® o 


► ® ft  n o 
) r)  CM  o — 


o in  — 

® 4 U) 


N 'J  n o ( 
CM  CM  (M  CM  • 


K •-  K O ® 

K eg  O — — 

<J  • 4 O C J 

• 4 • • • 

o ^ o o in 


— ftO 

— <M  • 

— o<j 

• • (j 


ocM  4 in 

ft  O O CM 


u)  n — 

CM  CM  CM  CM 


o n — u> 

CM  eg  CM- 


«“•  ® O » 

4 o O'  • 

• • • • o 

® r*  o o o 

^ ^ 0.1 


ft  o — — oo  — ooooorgn 

o in  I I I I I I 1 I I I I I 

o • UJ  UJ  UI  IJ  UJ  Ui  UJ  Ui  UJ  UJ  UJ  UJ 

oiocM  • nin<f4®NO  — ®CM  — nn 

•®  — ^n  — <toocnftfM  — 
oo— otnn<tcM  — in®N 


■ uiin ® < 
’ o UI  at  f 


n o rw  If)  rj 
o n o o N 
4 • o o 

• CM  • • • 

0 — 000 


*r>  N n i 
«t  ® 
0—4' 


o c>  o n 

fg  — O <f 


® CM  ® rg 
in  O — CM 


oo  < 

Q 

— CMn«yoONO 


n o n — 
n CM  CM  CM 


®NO  — 

CM  n 4 o o N o 


NS 
N U 
O 

N < 

• N 

n 4 

o 

ft 


CM  cr>  <r  N 
N o N <# 

• • • • 

on  — CTiNOCM  • 

rg  eg  CM  --  — at 


(M 


nc 

I 

Qj«ot  — on  — ftnn  — CM  — 
o 


CM  aj  n o 

oft®® 


N ft  r 
CM  O CM  ( 


N n — ® f 

CM  eg  CM  — • 


• — O 
I I 
< UJ  UI 
5®  O 

N 4 

f o n 


n 

s s i 

CM  ® Q C 
CM  4 
O ft  < I 

• ♦ N 

®n<  - 

o 

— fgn^ooN®ft  — fMn^ooN®ft<cvj 


— o — — — ootjCMn 

I I I I 4 I I I I I 
UJ  Ul  iLiUiUMWUlliJlUUJ 
n «t  g)  o m (»  rg  o rv  ft 
5fNN— N»^Oft(\IO 
^0®®C»ft^NC''4 
• ••••••••• 

ON®®  — <Mn  — on 


u 

< 

Ci 


UUCMJUUUUUU 


RE*D  I CO  1 .RM  1 , SR.<  1 . EM  1 .SEMI  CCA5C550 

1001  ECR“*T ( 4 ( 7X ,Dl 3 .6 » ) COAS05&0 

1002  FORMATdH  ///IH  'Ml  = (RMl  ♦ SRMIAWL.  EMI  *•  SEM  1 * WL  ) • / 1 H , C0AS0570 


L-ISTINII_:r_CCSTFCt. -CAPO  ✓ INPgT.PATA  £APC_Ng 


C^O  — ^J^O<Mr)^ON<OC^O— Cgn-tfU^ON 
l^)«0'£>’O^D'£>^f>OO^O^“^'^*NNNNN 


OOOOOOOOOOOOOOOOOOO 
co<^o— (Viri'fiT'sCKcroo  — 
U^^O'0'00»0'C'0•00'£)'0^‘NSN^*^•^• 
OOOOOOOOOOOOOOOOOOO 

<<<<<<<<<<<<<<<<<<< 

ooooooouooooocouoco 

UUOUUUUUUUUUOUUUUOU 


X 

* 

n 


UJ 

*-* 

-«• 

• 

• 

M 

1 

1 

w 

w 

*« 

m 

u 

T. 

1 

1 

UJ 

<— 

— * 

• 

•— 

» 

w 

w 

X 

(D 

u 

w 

w 

• 

X 

X 

J 

o 

• 

«>% 

» 

n 

••ri 

1 

1 

Z 

UJ 

►-- 

UJ 

• 

w 

W 

in 

• 

(D 

u 

♦ 

1 

1 

II 

J 

z 

» 

» 

UJ 

< 

< 

• 

z 

w 

J 

a 

N 

« 

N 

# 

in 

•m 

N 

^ « 

• 

— » 

-J-t 

» A. 

o 

O 

•>4  «-» 

3T  z 

X X 

V 

H 

1 

►- 

1 o 

• a 

in  X 

• 

•-• 

•-•  UJ 

cr  m 

» o 

a 

O 

w 

O 

w 1 

UJ  ♦ 

o • 

< 

o 

o • 

cr  — 

• ro 

< 

1 

1 < 

D* 

D — 

o 

^ UJ 

uq: 

^UJ 

u 

J 

-J 

w 

IL  ^ 

UJ  • 

u. 

» 

w X 

w X 

• • 

w 

o < 

< o < 

o u 

- J 

<• 

O'  • 

w C> 

• 

a 

z a 

M 

K 

• H 

♦ • 

1- 

*f  2 

UJ  1 

II 

< 

-*  U 

tr  u 

o 

11  • Uf 

"•  II 

• UJ 

wD 

OO 

--  I 

u 

D 

1 - 

1 II 

UJ  M 

JUJ/"  ^ — Ct^  2*^ 

or  i^nr  X ^ o ^ o a >-  uj  ft 

• - -rf\ju»-n-«DcraD 


h-  j ^ 7 i\  w i'’  u o » o *-  n 

uj  4f  ^ a u-  o < o u.  o o u iiJ  u uj  -r 

uuj<0'-'UujO‘-uuju.Q:ujiLti  u> 

- N 

ON®  o N m 

— N N N 


152 


USTInG  0«"  CONTWQu  capo  / INPyT  PfTA  gAPp,  NQ 


— (^Jro<^u^'0^‘®c^o  — Nn«jin'C’Nc0C‘O  — 


OOOOOOOOOOOOOOOOOO 

0000000000-^—r-*--«^«^ 

OOOOOOOOOOOOOOOOOO 

UUUUUUUOU^JUUUUUUUU 

OOOOOUOCJOOOOOOOUUO 

UUUOUIJUUU'JUUUUUUUU 


OOOOOOOOO 

cD(7o  — cvjf'xtino 

^^CMN<'JP>i<vjCVJ<N< 
OOOO  <5  0000 
UUUOUUUUU 
<<<<<<<<< 
uouooaoutJ 

OUUUUUU'->U 


oo  o 

N ® O 

<\JCsl  <NJ 

o o o 
uuu 
<<  < 
aoo 
uu  u 


o o o 

O ^ f\J 

r>nr» 
ooo 
uuu 
< «<  < 
o uo 
ovu 


ooo 
r>  u5 
n ^5  D 

ooo 
uuu 
< < < 
ooo 
uuu 


ooo 
N o 
ro  r>  r, 
ooo 
uuo 
< < < 
ooo 
ouu 


ooo 
a o ^ 
n<t  <f 
ooo 
uou 

< <f  < 

ooo 

uuu 


ooo 

f\j  r-|  C 

^ <5- 

o o o 

uuu 

< % < 

U Q O 

uuu 


ooo 
(P'ON 
<r  <f 

ooo 
uuu 
< < 

oo  u 

ooo 


o o oo  o o 

o O O --  o 

< «J  U")  IfJ  tD 

o o o o o o 

O O «J  o 

«(  41  «<  < < < 

oooouo 

ouuouu 


o 

z 

o z 
O *.^UJ 
o(Nicr 
0X0 
UHU 
< < 
OOi/i 
UUUJ 
M H 
o UJO 
uj  z a 

® 3 

•wo 
X _>u 
w o 
ocr  M 
o)  rn  3 
II  n w 
z m < 
o a 
— o 


OJ 

q; 

o 


O'  o 

O N 

forg 

• • 


® ® 

• N (M<> 

® O rg  nj 

• • • • 


o • • • 

• «t  4f  — * 
•O  A}®t»- 
O ^ f>J 


0 

1 

UJ  • • • • 

o O*  ® 

CO  ® U5  lT> 

w -•  10  r) 


V 

*4 

n ro  o N gj  A# 

4f  in  CT;  w lO  • 

ro ^ ^ ^ 

O' 

P P \ »4  • P P 

fg  fg  p p p \ 

a 

n 

cu  04  OJ  P --  >0 

p p c>  u>  iz  1 

O <J  rg  o — 

o 

• 

f-4  04  -A  U > O ® OJ 

o o «t  g>  eg  O'  UJ 

1 -4  ^ m o 

u.^ 

r> 

<t  ntj^  p CO  rj  lO 

• — OJ  P — O'  o 

UJ  — ^ r)p 

o 

p 

p p p eg  p • • 

vC  • • • • • • 

p p p p 

• 

• X 1 

c» 

ro  r)  cvj  r^.  ^ 

fg  to  C/  N «0  P ® 

* ro 

— 

— 

v> 

UJUJ  < 

o 

o cn  p p p • 

— 

*— 

UJ 

0 

crOA^ 

® 

— -<4  0.  m 0.  Af  oj 

OO  < o — 

1 p —0  — 04 

— 

— 

u 

Hi 

oz  ^ 

• 

o N P Z)  o eg 

• -A  O — ® "•  1 

UJ  P O O <f  o 

1 

1 

(X 

u*--  * 

n 

< r»  O'  vo  N 

® •pAjfgPuJio  •O'i’PP 

•A 

M 

o 

o 

• -J 

p 

p p p ^ p p p 

fg  10  • • • • O 

p 

p 

c 

p 

A^ 

S4’ 

z 

o 

UJlL  UJ  < 

N 

n P rg  • N vO  < 

• • GU  N ® <t  p 

4f  1 - — 

u 

m 

Cf  CJ>  UJ 

N 

p p p o • • • 

O ^4  • p p p N 

P |l|  P P P p 

1 

I 

oiL 

o 

N 04  <#  p ^ 

• ^ <T  O <0<Z  • 

o O ® ^ N 

— 

— 

UJ  z o 

U X w • 

p 

O O -•  ® « O p 

O p P fg  o CO  fg 

1 O O N P to 

w o 
< z 
o < 
ua> 

uz 

lO 


u.o<  r • 
oza  I 
a-LL  u-  ® 

X UJ  A-*  c^ 
iu(i.  cr  o • 
oo  — »o 
Z UJ  * • 

•-  W I X - s 


U)<fO  — — CD(?  niOUilO 
«> 


n ro  ro  o ® o *«• 


c^  c-  ^ 
UlMin 
iO  <?  o 

• * • 

fO  r)  n 


N n ^ 

•CVJ  C?*  ^ 

o • • • 

• CO  10 


lO  k UJ  CO  o uuo  n 

• • • • UJ  o • 

® N 'O  o • o • 

• • • • • ^ 

•<fN^OOD  •O- 

O'  ooo  • o I 
>o  o u)  e cj  I UJ  I 
• • • • I UJ  ® ^ 

® N W5  <t  UJ  10  lO  ' 


>^C  rn  ^ 

I --  e «o 

jujinri  r; 

)N  • • • 

) ^ 


w 

X ex  w UJ  o <0 

p p P o • • • 

— gj 

p • p 

o or  a 

UJ  • u.  < J - 

Zl 

® — -g-  eg  ,0  o o 

^iC^vOvCOO  *P 

p • 

z>  — 

zoor 

o • UJ  f3  a u.  a Q 

p 

ON®  p o ® < 

p pp  — o — ® • 

« O lO  O lO 

— U.O 

7 X (X  O Z • 

p 

® «J  o O O P 

ON®Wpoj  PO 

p p 

At  P P 

w ® 

-UJOW  >.  o — 

p 

p • p • p • p 

♦ - • • • *04  » 

o — 

p p p 

< O D 

O u cr  a w tj  — 

<t 

p p r>  o ® N ® 

o • ® N p »n  1 UJ 

1 1 

— — — 

o •< 

UJ  Z < < QC  * vO 

p 

* p p p p p p 

lOlO  • p * *UJOUJUJ 

• p p 

U A- 

>-araz<-A^ 

p 

O eg  ® • p A#  fg 

p #o®P^iO  • 

O O 

-o  P 

O' 

V 


GO 

I 

-I 

< 


••  o — a M u • 

W UJU.  JO  1-0 
U>  < • 

< — xujxa 

awujo:«<uj<f  • 

U.  O Q z W 
uj<zcra  — ®£D  • 
cra*-ooo#  •ni 
u.  u.  u.  < • 

► OUiU.  X < • 

-Zao</>o*^uj«t 
- *-•  UJ  UJ  a A-'  ® '**  •>  o 

< •lO<ffU<7' 

'<orcruur  o*N<t- 
j o o < •/  / t-  z • n • • • 

j u u a UJ  UJ  woo  • p)  r<  n 

J OOtO  — *-N  o • • • 

'H  |I.JOOOUO<IOaJ'0® 
: <ujuiZ  — 7 w®on 

► N N 


lO  ® O O 'O  lO 

O «t  — • N N UO 

• • • o • * • 

D o ruo  ® w lo 


OWNO^yO  •r)0-*C7'f'J® 
♦ — owconcg  *>0  — r><tP) 
in 
^ < 


-■•  rg  o o o N 
O N lo  • vO  o rg 
w o o fO  s 

• • • o • • • 

n r>  n w ® 


jui  w UJ  rr  u.  u.  3.  a y w fO 

fiJiixa  i-<  non 

no  O O M H H u — o o 

“)<j<JO  — cg®4f 

<nw  »u-u.  <ff)UO 


o o o -o 
• o too 
o • ^ ^ 

o o • • 

W ® 

\ • . »\ 

CD  »0  CM  O O 
fg  O a» 

< ♦ lO  o 

w o • • • 

< on 

o 

-•  (\i  n ^ 


ON  OlOCNJ  I 

• ••••*  i UJ  • 

OOOJ'O'JOUJlOO 

• — NN<rCgo  •! 
o • r>  oj  CMO  • eg  UJ 

rt  •CT'NOlO  •O^ 

• 0^  ••  ••0|N 

o^mn®-«  I lij  • 

• 4fgin<na)ujo<j 

oo'00'og?o  • • 

O**«*4*0JO 

I o o N N ® in  • I 
.rg  • • • *N.oUj 
I • cj  o cr  ® o 1 N 
o o o rg  UJ  ® 
• {X)org®«finN 


a 


-J 

J » 

» • 


10  >0  N 

I ro  o ® 

UJ  P5  <♦  n 

r>  • • • 

«t  t-« 

«o  • • - 

• ^ O'  o 

a>  ai  — <» 

• (M 10  r> 

•-•  • * • 

I 

UJ  • • • N 

CO  lO  u>  z 
lO  oj  N U5  U)  O' 
rj  <g  u)  n o U 


tx 

w 

< 

o 

u 

u. 


a 

UJ 

or 

o 

u 

u. 


o a 

^rgrOAj®oN 


w o 

< • 

o -- 

U II 
K Z 
> w q: 
cr  < o 
K n 

Z U UJ 

UJU.  fr 


o 

N 


O 


_l 

— < 

^ • 


CJ 

N 


U 

I 

-J 

< 


« 04  « 

A"A  04  -A 


O -'•a 


I (1  • Ul 
I 1-^0 
.)  ►-•Z 

: > • o s,  A-* 

) nr  II  — ® w 

J Z • P u.  o 
' Ul  < o — u 


J 

^ » 
< - ^ 

+ •!- 

J 

U 


I o 

•^UJ 

^ I 

o • 

\ < 

->  UJ 


W X 

O J 
wCL 

3 

- --  u 

.,  •UJ-O 
1 I (I 

« -z^a  z 

w o W-.  w UJ  O' 

«r  eg  u w C5  or  o 
u w z i»  o w o 
<nu  oouuJZ 

UJf5-*UWU.ttUJ 


o M ® o eg  ® 

-4  — og  04  04 


u 

< 

<4 


UUUUUVJUUUU 


153 


— rsj  fo  <t  iTJ  »£»  N o>  (T'  o — N n <»  h>  cu  V *0  CO  (;  o oj  r<  «j  tf > -i)  ^ cn  o t > ► •«  ^J  n < tn  n 't)  ^7*  r>  ^ rv<  ^ O ^ oo 

— -4  ^ *-i  fj  ( j cj  rj  cvj  f'J  c\/  r\;  Ci  c\j  r'.  ri  fi  r)  r»»  r)  r>  r)  D <j  .j  <j  <f  <7  < ^ «7  «/  <7  lO  in  lO  m in  tf)  IT.  lO 


OO'^OOOOOOOoOoOOOOOO 
o — J"-  C'  ^ cr  c^  o — ro  <J  in  ‘f5  ^•  ^ 

o o c o c c b C'  C-  C — -<■  — 

O O O O O Ci  o o <j  o 0 o o o o o o c 

crcrc'CGr.'C'rtaan''Yfrc:frfrQcrcr 

U.ULU.U-li-l4.U-U.U-j.UU  li  U.U.U  IlU-U. 


o o — r . r";  n n .j 

■-■-■  Tv." tr rjT<3rj  Cw  o*  ro*r;(  ’ •'T^rnTyr/ 

r,  O < * O O O C.»  O O O t j O f ><*»<  ) o 

ct  cr  f/  n :r  cr  f/  cr  <r  t rr  f'  ir  a r 

< < < < <<<<■<<*'<•■<<< 

Luu.  u.axu.u.xiLu.  uii-uuu 


m «n  ^ j)  c?'  o c 

r"»  ^ ; TV  rj  ry  <j  vT" 
C>  O O O O o Tj  { 

cr  rr  cr  or  cr  rr  n:  i 
*t  <1  *r  < *t  * 
L.  U.  U.  U ii.  U-  IL  I 
J -J  _)  ^ ,J  -J  ~J  . 


• 

• 

Z M 

• 

r> 

Lu  y 

w 

lu  u» 

0 

UJ 

Q 2 

u 

• 

"W 

T in^  ' 

0 UJ 

-1 

H 

3 1 in 

UJ 

w 

o 

CD 

• 

a -J^- 

< 

X G.  m 

-J 

o-' 

»-  C5  -J 

• 

— 

Z 3. 

w 

tiJ  u 

in 

> 

>23 

-J 

— 0 

Ui 

>-  in  X 

> 

in 

_j 

< 

»-  U'  U-’ 

U‘ 

-J 

in  _i  0) 

Oj 

• 

a z *f 

< 

.~- 

-DU 

_J 

U. 

w 

T 

w 

UJ  < u 

in 

u 

XI- 

J 

X 1“  X 

- or 

II  '.u  7 i- 
y I-  — J 

C7  U U) 

z < *m 
X «* 


<1 

- 

— 

1-  (\j  a 3 > 

Kl 

■— 

w 

UJ  < • 

<tl 

0 

X 

to  • a 0 K 

a 

<f 

-J 

UJ  U 3 ^ U 

1 

W 

u 

s U U.  <?  ii. 

UJ 

D 

— • • iT*  X *t  X 

Dl 

3 

< 

• nr  3 0 . 

ai 

»- 

3 

C 0 O'  C 3 — X 

y ►-  ^ 

o • - 

oz-  - 

O — u 

OU--  -J 
OC  V H 

cror  . « 

< < K 

U.  L — ^ 
JJwO  ^ 
It  y..U‘  — 

C.J  w'J-' 

uiz  z r / o 


N I-  ti’  — • • 'vJ 

m 

‘J  L’  il>  UJ  -J  < •-•  <M 
• — h-  C LU  O l/> 

vozax-.  c 

• 2 — < V >•  I- 

— 'iS  ^ <f  < 

>-  I • < 7 O 

w 2 i-  X in  X o 

— < _j  — • 

-2  ZUJSZ^ 

V -•  n O <t  *-  o n. 

V a 

cj  .j  *r«  X o uj 


o >-  UJ 

- *u  j 

— o u.  >- 

fg  • X ►* 

^ r H 


m — n .jvij  .1 -T 'i»  V •—  • cr  o 2 n < 0 

• » *-u;  Mi*  I r un  f w •—  ju  < 

X 3 in  c:  <•  .n  O’.  2 h-  -* »-  w »--  y 0 a o n 

k-r)7<>'f  *r<U  V.C. 

err  in^  ►••  _j  »t  in  ►-  t-  n j r > / .j  j j c . 

OC  > “)  » • — lij  L I *r  f/  M (7  . J . J J H <7  • 

II  3 •-^  c < ininuKJU(,u  ^ < 2 II  I 


I t- 
D X y 

O W 

O J 

^ —o 

c in  rr: 

• M w 5 

II  :T  II  > 
7 Him 

t - n -J 


i H 

O XX 

o u — 

O . J 


1 un  o ^ c f ^ — CO  "V  <t  IT.  c ^ 

• <r o‘->tnmwgv  - 
) O o O O o C’  C-*  O r ) o r>  o o o 
'(yf/crticrcrorfrfrnj’afrfr 

.U.U.U-L.U-U.  UL.U.U.U.U.U.U. 

_l 


<n 

UJ 

u. 

X 

r- 

X 

UJ 

UJ 

K 

a n 

X 

Q 

tn 

< z 

UJ 

•“* 

-J 

*4  • 

u. 

u < 

0 

in 

< 

1 in 

0 

N 

z 

*-  UJ 

Uj  — 

in 

»-• 

X 

UJ 

CP  f 

Z 

li. 

a < 

CIj 

l J 

— 

►- 

0 z 

z 

7 

UJ 

> 

\- 

V-  — 

> 

X 

r\j 

in 

4- 

7 

1 

►- 

II  0 

3 

UJ  u 

X 

X 

0 • Z 

7 

> z 

1 

< 

3 3 

»-  r\i 

< — 

X 

7. 

< w 

Z II 

X 

X « 

< 

X 

> »- 

3X 

u 

U) 

7 

w 

u z 

1-  a 

.»■< 

3 |~ 

X 

» 

• 3 

UJ  mil  z 

X 

0 u 

•»- 

z 

0 

lO  z T ul 

« 

V Z 

• 

— 

Ul  — 

UJ  1-  X 

2 

7 

►-  II 

z ^- 

in 

tlU- 

— 

X 

y J 

• z 

>- 

uy  0 

• 

a 7 

c n^-  a 

u;  7 

0 

< UJ 

M a UJ  h ^ 

UJ  U 

O' 

y Z 

z 0 1-  0 

X — 

» 

• 

K .J  >• 

in 

> 

tJH 

1-  in  (5  a 

-JZ  • 

UJ 

Z -jJ 

UJUJ3  0 

•u  UJ  u 

n 

— m 

i/ix  Li  at  y 

3 > Z 

» 

< h < 

z •>  — 

in 

TJ 

ui2 

< 0 1 

» 

IT.  • 

• U CT  M X 

3 TJ 

— Ct 

• 

(T  C>  0 

UJ  in 

• 

• ocsiD  0 Cl  UJ 

UJ  > — 

>-  H 

ll 

Mh-»-  X 

X < 

n 

» < 3<?  ‘U  ►‘•K 

K X 7 

— 

X u a um  m 0 

D 

w 

• — 

0 

- >■  U. 

2 i/»- 

UJ 

-J  az  n<A  < c 

— > 1- 

3 

azoiz  — rrin  • 

< •— 

h* 

c — 

— J aa)Z  — 

y a y> 

z 

1-  < 5:  0 

a 0 o 

— o n ^ u • n « »-  • II  T-) 

mrri  — occinr'in  — — into 

Mw5  Li 

r II  > I UJ  -«  II  > X .J  X X J 

uiin  K 7 u:  m *i  n H in»- 

n.j  — ioo  — si  I o — * 

O h-  I — I jX^  — »-  + .JXO  — l- 

_J  w --  w I w •-«  J kC  — • 

1 1 *-  < M u.  II  u n c-  «i  < n o u 11  <1  » 11  n u. 


H il'khX  r>*>^uO<c. 

- 2 o UJ  c.  2 < z n 2 J 

O •_JlT<rC>“  Z •-  X2 

• CO.  a *“ 

o HOZOiciaJuzn 

- zzi^  a h-zujzaj 

O — c^xj<cc 

• H K-xcrrui  ♦ 

O h*  UJ  u;  O ♦-  m UJ  10  I-  V O 

w »-  X n Z I ► - O J 

m C X 2 K lu  in  ►-  n 2 jJ  3 

« _)  < iiHo  rr  z in 

X a '4  < < IV.  m n I-  'j; 

< 2 ^-..,-4 

cocTxn  V 7f-z 

I _j  It  3 cr  J •?  in  tij  — 

— _i  — »-»-an'  — xnw«^a-»3  z 

II  < -■•  u;  7 o o h-  7 3 u.  y »■  - X 


'71 

z in  n UJ  r.  / 

0 u.  a u. 

— L'3U 

— » — — 

. 0 X u w — — -4  0 X U 't  — 0 

- — U < “ 

c*-*  — u — a u 7 

cv 

71 

u - ^ 0 

c 

a 

a 

Ui 

1- 

— 7 

0 0 

in 

a 0 

f\J  0 — 0 CM 

n<t  0 

iTsC  7 

z 

u u u v-i  u u i J u l.  J u V J 


J54 


OC7‘O--<\inci/';O^^C00O-“M''>C  iTJOn-CClT'O  — (g'OctO 
^^CO'000'OJ^DO'£^“^•SN^'^•^-NN^•X)CCC*CDDa5COcc^L«COO'^;'C^^>‘0■C^C^C^C^^OO 


cj  O'  o — oj  n in  o r ■ cy  a o — rvj  r,  <;  m .f  ''  'O  <>■  o ‘'J  f*  ■-:  Ln  ^ n 'd  c o — vj  m n n m r>  o 

. lO  i/:>  c-  o o </  c h-  N ^ r rv-  f n n tn  x c r x x x>  x x x a c*  o-  a i>  (T  cr  ir  o*  o o 

c S rr  ? C£  rr  cr  cr  cr  rr  fr  a a £/  ry  O'  cr  ti-  cr  ^'cr  ti-'' E?  ty  ff  -cy-fnc.  t-' 

<<<<f<<<<  <<<<<<<<<!«{<<■<<<<<  •?•!<<<<<<<'<  <«f 

L.liaUUU.ll-LLU-li.U.li.Xll-U-li.U.U  U.  aii-iLXXlLU  UU  HXIlU.  U L LLU-U-U.  XU-XU.  x 


a H 2 • 

a 

< I 1 — 1- 

X 

Q 

»r  U z < 

7 

u»  1 z y - 

II 

» 

X u or  o o 

U) 

O 

♦--  z UJ  X u. 

u 

• 

UJ 

^ o 

j.  « ir  uj  i-  z 

»- 

—I  •• 

o y ui  ►-  — 

z 

<’  _J 

;■)  ^ 1 rr 

o < 

o >-  u 1 z u a 

• 

• 

o > 

or  a . — ' z 

z to 

Ul 

• u 

X — fy 

'j’  .j 

y 

o 

H > Z U -fa  UJ 

o 

< 

w if) 

UJ-*  111  y (T 

uc. 

z 

• u 

u »•  z > 

(T  y 

y • 

z o \.j  z y D 

y 

o 

c -* 

►-  uj  uj  o cr  z 

y to 

tl 

• N 

z H Z U Ul 

< 

o 

^ to 

z HZ  > or 

? U' 

o 

_l  — 

<J  D U Z UJ  O 

L 

u 

— 

< w 

u a < H 

y X 

1 

o 

y 

01  I/)  <t  o < 

»- 

y 

UJ 

U^  \ 

<1 

z ar  u’u 

n 

y 

K| 

Z i/)  ^ C I- 

H 

X 

2 

to^  y 

<1 

• 

o 

*-• 

CD  < *“ 

o< 

• a tn  o z 

2.  y 

- 

— . < > • w 

r/ 

1 

ctt  - 

C Ui 

to 

V.  2 

a 

1 •uT.y 

1 

HI 

Ui  C C 1/1 

y y 

XI 

< 

«t 

M OwO.  - 

• c 

oiy  CcOLJ  2 

< 

X 

> 

y 5 

y 

o Z tA/  < 

OC  lO 

ai 

^ — z in  < 

U UJ 

UJ 

Ul 

J 

y > 

• 

z<a'.cr^- 

• CJ 

ZI 

u y — — 

— r 

C 

c z 

II  II 

u 

• *T  O to  r > 

o 

y •-  to  u:  _j  X 

Oy 

•~» 

— o 

> 

«— 

z 

^ ^ w ►-  CJ 

1 O 

1 

z o z j o 

Z 

to 

to  — 

•— 

O O o Z *-  CJ 

• y 

Ki 

XI  Z < “"  < C ' 

--  a:  *- 

S to 

to 

* 

• • « o X 

o 

1 

? - j r,  X 

L N 

-«-• 

^ ►- 

n 

y y 

Nl 

a so 

c • o 

Of 

UJ  • • ^ n 

^ X to 

z 

? * 

2 

c «-o  y 

• lO  o 

a\ 

fr  UJ  o < o y J • 

cj  — 

•— 

K-.  ^ 

D 

• • -J  _J  X 

<1 

u o • • to  «r  3. 

2 >' 

z ^ ^ 

C) 

X 2 

• 

H y < «t  O 

• X 

ul 

z C'H  *y  > 0 

y < to 

-.>  1 

C ro 

y 

O J O 

OUJ  CM 

• M • u 1 fr  / < M y >'  > >"  z'  ^ i I 

un.if,  n n r 1 o «r  vz  lo  «»  < n • • * * r 

T J *1  *1  uj  • X — ^ ^ j — 

► V- V y vr  <■  1 ‘ Nj ->  i.;  X «j  > r i 

i/>  u »'  uhAk/)  y >^  7 «o  cn 


'f  r — c' 

Ul  V) M < 

.j  X — 

< UJ  (.■)  M 

u 


\'.  X «j  > r 

/ w I 

■«  u ' 
1^1^. 
^ * ^,l  N.  Ki 
I/)  im/j 


I O _J  • / ♦ 4 CM 

^ n »sN- 

5 / U X N •-• 

-ww  — (T 

- y >-  ► o i — — II 

/ .111  h ^ 

J *>-  ►-  U ’1  U‘- 

o: 

oo  cv 

^ tf) 


.1  Ml  -rr  <\ 

< < *«  ft  tif  • 

> > I uj  n o 

o u f»*  a r 'll 

^ If  . -i  • 


. y.  J I- 

< u II' 

. . - u • * tr 

^ If 

0.  o 


--  O I 

w m UJ  ^ 

y - *-  • H 

c»i-  j J 

•ICO  Cl 

X . i n z’ 

O C IJ  / 

uicr  a 

1 j / n.j  r> 

_j . f — -J  I-  n 
< < U lu  *t  Ui 

u o — ar  u rr  ur 


155 


